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ABSTRACT 

Ali,  Galal  Abdalla.   Ph.D.,  Purdue  University,  August 
1972.   A  Laboratory  Investigation  of  the  Application  of 
Transfer  Functions  to  Flexible  Pavements .   Co-Major 
Professors :   William  H.  Goetz  and  Milton  F.  Harr. 

Transfer  function  theory  was  applied  to  examine  the 
behavior  of  flexible  pavements.   Two  convolution  techniques, 
one  numerical  and  the  other  analytical,  were  employed.   The 
numerical  convolution  was  used  to  derive  pavement  response 
functions  from  impulse  testing  of  three-layer,  flexible  model 
pavements.   The  analytical  convolution  formed  the  basis  for 
calculating  deflections  resulting  from  static  and  repeated 
loads . 

It  was  hypothesized  that  more  significant  parameters 
than  commonly  used  could  be  obtained  under  controlled  labora- 
tory conditions.   By  using  non-linear  regression,  the  response 
functions  were  approximated  by  a  mathematical  model  to  include 
these  parameters.   The  derived  model  was  used  as  input  in 
the  analytical  convolutions.   The  adequacy  of  the  developed 
models  was  verified  by  comparing  predicted  and  measured 
deflections . 

Curve-fitting  of  the  dynamic  peak  deflection  data,  by 
the  non-linear  regression,  resulted  in  a  mathematical  function 
for  the  deflection  basin. 
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A  silty  sand  subgrade ,  a  crushed  aggregate  base,  and  an 
asphalt  concrete  surface  constituted  the  components  of  the 
three-layer  systems.   Each  layer  material  was  characterized 
by  a  series  of  conventional  laboratory  tests.   Model  pave- 
ments, of  two  different  surface  course  thicknesses,  were 
tested  statically  and  dynamically  at  three  different  stress 
levels.   Test  temperatures  were  50,  75  and  100°F.   Deflec- 
tions were  measured  at  five  spatial  locations. 

The  following  are  the  primary  findings  of  this  research: 

1.  The  profile  of  peak  deflections  of  a  flexible  high- 
way pavement  can  be  described  by  the  equation 

t    ^  -Dx2 

y(x)  =  y0e 

where  y(x)  is  the  deflection  at  a  distance  x  from  the  load 

center,  y0  is  the  central  deflection  and  D  is  a  constant 

reflecting  the  attenuation  of  the  deflection  profile  with  x. 

2.  The  central  deflection  y0  increases  with  increases 
in  temperature,  and  decreases  with  increasing  surface  course 
thickness.   Increases  in  temperature  (50°F  -  100°F)  or  surface 
course  thickness  (1  inch  and  2  inches)  decrease  the  value  of 
the  parameter  D. 

3.  The  time-dependent  behavior  of  a  flexible  pavement 
can  be  represented  by  a  response  function  G(t)  which  is  a 
function  of  time  t.   It  is  possible  to  obtain  this  function 
from  a  single  impulse  test  on  the  pavement.   The  response 
function  is  independent  of  the  magnitude  of  the  impulse  load. 

4.  The  response  function  of  a  flexible  highway  pavement 
is  of  the  form 
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G(t)  =  ae"et  sin  yt 
where  a,  B  and  y  are  descriptive  parameters  in  the  response 
function . 

5.  Parameters  in  the  response  function  are  believed  to 
be  pavement  properties  that  provide  for  better  understanding 
of  its  behavior  than  those  currently  used. 

6.  Temperature,  surface  course  thickness  and 
spatial  location  have  their  respective  influences  on  the 
response  function.   Increases  in  temperature  increase  the 
value  of  the  a  parameter  in  the  response  function,  while 
increasing  surface  course  thickness  or  the  distance  from 
the  load  center  decreases  the  value  of  a.   The  6  and  y 
parameters  do  not  seem  to  be  affected  appreciably  by  the 
above  factors. 

7.  Transfer  function  theory  is  capable  of  predicting 
static  or  repeated  load  deflections  of  flexible  pavements. 
The  favorable  agreement  between  predicted  and  measured 
values  of  the  deflections  in  this  study  validates  the 
hypothesis  that  the  parameters  in  the  response  function 
are  material  descriptors  which  are  independent  of  the  type 
of  load  input. 


INTRODUCTION 

The  pursuit  of  a  rational  method  for  designing  and 
evaluating  pavements  requires  considerable  knowledge  of 
their  response  behavior  to  realistic  loadings.   The  ability 
of  a  pavement  structure  to  withstand  such  imposed  loads  can 
be  expressed  in  terms  of  primary  and  ultimate  response  modes 
[1,  2] .    In  the  primary  response,  the  magnitudes  of  developed 
stresses,  strains  or  deflections  can  be  used,  whereas  failure 
must  be  considered  in  the  ultimate  response.   Hence,  an 
essential  part  in  pavement  design  is  the  development  of  a 
procedure  to  analyze  a  pavement  system  for  representative 
input  stimuli.   Another  necessary  requirement  is  to  define 
a  failure  criterion  in  terms  of  the  ultimate  responses,  for 
example  by  limiting  deflection,  so  that  the  designer  may 
assign  allowable  values  to  these  quantities.   When  deflection 
is  considered  as  the  criterion,  it  must  be  supplemented  by 
other  measurements,  such  as  the  slope  of  the  deflection 
profile.   Knowledge  of  the  analytic  form  of  the  deflection 
profile  is  a  prerequisite  to  determining  its  slope. 

Although  the  question  as  to  what  causes  failure  in 
asphaltic  concrete  pavements  is  a  complex  one,  several 
failure  mechanisms  have  been  postulated  [3,  4,  5,  6], 


*   Numbers  in  brackets  refer  to  entries  in  the  Bibliography 


Surface  cracking  and  rutting  may  be  caused  by  the  effect  of 
direct  wheel  loads,  shear  displacements  and  further  compac- 
tion of  the  base  structure,  or  due  to  irrecoverable  deforma- 
tions in  the  subgrade  material  resulting  in  creep  or  plastic 
flow.   The  elastic  or  recoverable  portion  of  the  deformations 
in  the  bottom  layers  under  the  action  of  repeated  loads  may 
result  in  fatigue  failure,  whereas  insufficient  subgrade 
shear  strength  can  lead  to  shear  type  failure.   Brittle 
surface  cracking  may  be  caused  by  non- load-associated 
stresses  such  as  those  effected  by  changes  in  temperature. 

Current  mathematical  models  to  predict  pavement  response 
and  to  evaluate  its  behavior  generally  evolve  from  elastic 
and  viscoelastic  theories.   These  models  incorporate  geo- 
metric aspects  and  material  properties  within  a  system  of 
equations  which,  in  turn,  are  solved  so  as  to  satisfy 
selected  boundary  conditions.   However,  any  model,  required 
to  be  predictive,  must  be  inclusive  of  realistic  in-situ 
material  properties.   Also,  pavement  materials,  and  particu- 
larly bituminous  concrete,  are  time-  and  temperature-depend- 
ent.  At  the  present  time  no  single  theory  of  pavement 
behavior  that  considers  time-  and  temperature-dependent 
characteristics  of  pavement  systems  is  used  in  designing 
or  evaluating  pavements.   Thus,  a  need  exists  to  obtain 
parameters  or  functions  capable  of  completely  describing  the 
time-  and  temperature-dependent  material  properties  of  pave- 
ments subjected  to  static  or  dynamic  loads. 


A  new  approach  undertaken  in  the  present  study  to 
evaluate  pavement  behavior  deals  with  the  primary  response 
mode.   The  basis  of  the  investigation  is  the  transfer 
function  theory  that  has  been  applied  successfully  in  the 
synthesis  and  design  of  various  systems  in  other  fields  of 
engineering . 

This  investigation  extends  the  transfer  function  con- 
cepts to  flexible  pavements.   More  significant  material 
properties,  or  indicators,  than  commonly  hypothesized  are 
defined  and  determined  under  controlled  laboratory  condi- 
tions.  The  adequacy  of  the  developed  mathematical-experi- 
mental models  is  evaluated  by  comparing  predicted  responses 
due  to  static  and  dynamic  loads  with  experimentally  measured 
values.   The  effect  of  surface  course  thickness,  temperature 
and  spatial  location  on  the  dynamic  parameters  and  the  pre- 
dictions are  examined. 

In  addition,  it  is  hypothesized  that  the  deflection 
basin  of  a  pavement  follows  a  normal  distribution  curve. 
A  mathematical  form  is  obtained  to  fit  the  deflection  basin. 

The  procedure  to  achieve  the  above  objectives  as  a 
pavement  evaluation  tool  is  depicted  in  the  flow  chart  of 
Figure  1.   A  known  impulse  load  is  applied  to  a  pavement 
system.   The  resulting  output  deflections  are  measured  as 
functions  of  time  at  several  locations  on  the  pavement 
surface.   The  peak  deflections  are  curve-fitted  using  a 
computer  program  NONLINR.   A  mathematical  function  is  thus 
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FIGURE   I      PROCEDURE   FOR    PAVEMENT   EVALUATION 


obtained  to  describe  the  deflection  basin.   The  parameters 
in  this  function  are  studied. 

The  time-dependent  input  and  outputs  resulting  from  the 
impulse  tests  are  analyzed  within  a  formulated  theory  to 
give  response  functions.   A  computer  subroutine  CONVLI  is 
employed.   The  response  functions  are  verified  using  sub- 
routine CONVLE. 

A  mathematical  model  for  the  response  functions  is 
evolved  through  NONLINR.   This  model  is  convoluted  analyti- 
cally with  step  and  repeated  load  functions.   Computer  pro- 
grams STALOD  and  REPLOD  facilitate  the  respective  calcula- 
tions.  Predicted  and  measured  responses  are  compared. 


BACKGROUND  AND  STATE  OF  THE  ART 

In  this  section  a  literature  review  is  presented  which 
develops  the  necessary  background  and  state  of  the  art  for 
this  investigation.   Five  aspects  of  material  characteriza- 
tion and  pavement  studies  pertaining  to  this  study  are  pre- 
sented : 

1.  Application  of  elastic  theory. 

2.  Developments  from  viscoelastic  theory. 

3.  Use  of  dynamic  models. 

4.  Transfer  function  approach. 

5.  Laboratory  studies  of  flexible  pavements. 

Application  of  Elastic  Theory 
One  of  the  earliest  applications  of  the  theory  of 
elasticity  to  an  infinite  layer  medium  was  the  investigation 
reported  by  Boussinesq  [7].   Later,  Burmister  [8,  9]  pre- 
sented the  stresses  and  displacements  in  two-  and  three- 
layered  systems.   The  layer  theory  was  further  extended  by 
Fox  [10],  Acum  and  Fox  [11],  Mehta  and  Veletsos  [12],  Jones 
[13]  and  others.   Generally,  all  these  theories  have  in 
common  the  following  assumptions:   (a)   Each  of  the  system 
layers  is  a  continuous,  homogeneous,  isotropic,  linearly 
elastic  medium  of  infinite  horizontal  extent;   (b)   each  layer 


is  continuously  in  contact  with,  and  supported  by,  the  under- 
lying layer;   (c)   the  surface  loading  can  be  represented  by 
a  vertical  pressure  uniformly  distributed  over  a  circular 
area;   (d)   inertia  forces  are  negligible;  and   (e)   displace- 
ments are  small.   That  no  pavement  layer  can  satisfy  these 
restrictions,  except  for  extreme  cases,  is  generally  accepted. 
Furthermore  the  above  theories  do  not  consider  dynamic  load- 
ing, since  the  stress  function  technique,  which  underlies 
these  methods,  is  applicable  only  to  static  loads.   Avramesco 
[14]  outlined  a  mathematical  solution  for  computing  displace- 
ments in  elastic  layered  structures.   The  method  accounts  for 
the  dynamic  phenomena  in  pavements  and,  according  to  Avramesco, 
can  be  extended  to  the  viscoelastic  case.   However,  only 
qualitative  description  of  the  dynamic  behavior  was  reported 
without  numerical  results. 

How  closely  the  material  properties  are  represented  in 
any  theory  is  reflected  in  its  predictive  capability.   The 
representation  of  material  properties  by  two  elastic  constants 
(as  is  required  by  elastic  theory)  is  simply  a  mathematical 
expedient.   From  laboratory  tests  on  specimens  of  bituminous 
mixtures,  Busching,  Goetz  and  Harr  [15]  determined  six  inde- 
pendent coefficients  necessary  to  describe  the  stress-strain 
relationships  considering  the  material  as  homogeneous  and 
transversely  anisotropic  when  subjected  to  small  strains. 

Despite  the  several  shortcomings  that  elastic  layer 
theory  has,  this  model  can  be  used  as  a  first  approximation. 
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However,  it  remains  to  be  determined  how  the  performance  of 
pavements  in  the  real  world  is  related  to  the  prediction  of 
stresses  in  idealized  elastic  pavement  systems. 

Developments  from  Viscoelastic  Theory 
For  several  years  the  component  materials  of  a  pavement 
system  have  been  recognized  to  have  time  dependent  stress- 
strain  properties  due  to  consolidation  and  creep  [16,  17]. 
Thus,  to  overcome  the  limitations  imposed  by  the  assumptions 
of  time- independent  elastic  materials,  the  viscoelastic 
approach  was  introduced  into  material  characterization  and 
pavement  analysis. 

Viscoelastic  model  representations  of  soils  have  been 
discussed  by  Schiffman  [18].   Lee  [19]  and  Lee  and  Rogers 
[20]  presented  more  general  methods  of  viscoelastic  stress 
analysis.   Secor  and  Monismith  [21]  studied  the  viscoelastic 
properties  of  bituminous  concrete  using  a  four  element  model. 
Papazian  [22]  represented  asphaltic  concrete  by  a  mechanical 
model  consisting  of  five  Kelvin  units  in  series  with  one 
Maxwell  unit.   He  pointed  out  that  an  infinite  number  of 
Kelvin  units  would  be  needed  for  accurate  representation  of 
the  material. 

The  development  of  the  correspondence  principle  for 
isotropic  media  in  1955  by  Lee  [23]  ,  and  shortly  thereafter 
its  extension  by  Biot  [24]  to  include  anisotropic  media,  made 
possible  the  analysis  of  pavements  using  time  dependent 
material  properties.   The  viscoelastic  solutions  are  based 


on  the  premise  that  the  time  dependent  boundary  value  prob- 
lem can  be  reduced  to  an  elastic  model  if  all  time  varying 
quantities  are  replaced  by  their  equivalent  Laplace  trans- 
forms and  the  elastic  constants  by  operator  forms  of  the 
stress-strain  relations.   Perloff  and  Moavenzadeh  [25] 
determined  vertical  deflections  for  a  static  load  and  a  point 
load  moving  at  a  constant  velocity  on  the  surface  of  a  homo- 
geneous, isotropic,  linear  viscoelastic  half-space. 

The  solutions  to  two-layer  systems,  where  either  one  or 
both  components  were  considered  viscoelastic,  was  the  subject 
of  the  papers  published  by  Pister  and  his  co-workers  [26,  27, 
28,  29].   The  problem  of  a  viscoelastic  asphaltic  concrete 
beam  on  a  set  of  independent  springs  approximating  a  Winkler- 
type  foundation  was  investigated  for  both  a  stationary  load 
[26]  and  a  concentrated  moving  force  [27].   The  case  when  the 
top  layer  was  assumed  to  be  elastic  and  the  subgrade  to  be 
viscoelastic  was  solved  by  Pister  and  Williams  [28]  and 
Hoskin  and  Lee  [30] .   Both  layers  were  considered  to  have 
time  dependent  (viscoelastic)  properties  in  the  analysis 
presented  by  Pister  [29]  and  Huang  [31] . 

In  response  to  the  demand  of  research  workers  for 
improvements  on  pavement  models  to  include  realistic  material 
representation,  more  elaborate  models  were  used  to  solve  the 
three-  and  four-layer  viscoelastic  problems  [2,  31,  32]. 
Barksdale  and  Leonards  [33]  presented  a  viscoelastic  analysis 
for  three-  and  four-layer  pavement  systems  subjected  to 
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repeated  loads.   An  attempt  was  made  to  compare  the  theoreti- 
cal predictions  of  the  resilient  and  permanent  deformations 
due  to  repeated  loads  with  those  determined  at  the  AASHO 
Road  Test.   The  measured  permanent  surface  deflection  for 
the  AASHO  Road  Test  Section  was  more  than  three  times  the 
predicted  value.   The  discrepancy  was  mainly  due  to  the 
assumption  of  an  elastic  base  course  made  in  the  visco- 
elastic  analysis,  while  the  results  of  the  AASHO  Road  Test 
indicated  that  59  percent  of  the  permanent  deformation 
occurred  in  the  base  and  subbase.   Thus,  in  predicting 
flexible  pavement  performance,  all  layered  materials  must 
be  considered  to  have  time  dependent  characteristics. 

Based  on  the  fact  that  pavement  materials  are  visco- 
elastic  and  temperature-dependent,  Ku  [34]  proposed  stress- 
strain  relations  that  accounted  for  temperature  variations. 
Although  the  temperature- dependent  relations  were  verified 
experimentally  [35,  36],  the  application  of  the  stress-strain- 
temperature  constitutive  relations  to  the  solution  of  layered 
systems  has  yet  to  be  developed.   Recently,  Moavenzadeh  [37] 
extended  the  primary  response  model  for  a  three-layer  visco- 
elastic  system  [2]  to  include  failure  in  flexible  pavements. 
The  developed  mathematical  formulation  accounts  for  environ- 
mental conditions,  such  as  temperature  and  moisture,  but  no 
numerical  results  were  provided. 

When  different  temperatures  are  being  considered,  the 
equivalence  of  spring  and  dashpot  models  are  lost.   To  obtain 
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the  temperature-dependent  constitutive  relations  for  visco- 
elastic  materials,  either  more  simplifying  approximations 
must  be  made  or  data  used  directly  [38] . 

Each  of  the  above  investigations,  within  the  framework 
of  viscoelastic  theory,  attempted  to  take  into  consideration 
the  time  dependence  of  pavement  materials  and  consequently 
provided  an  insight  into  their  behavior;  however,  the  corre- 
spondence principle  requires  that  the  associated  elastic 
problem  must  be  amendable  to  solution.   The  elastic  approach, 
in  turn,  imposed  certain  mechanistic  material  properties 
(constitutive  relations)  into  the  analysis  inasmuch  as  the 
equations  obtained  from  equilibrium  conditions,  strain  defini 
tions,  and  compatibility  relations  are  less  than  the  number 
of  unknowns  contained  within  these  relationships. 

In  both  elastic  and  viscoelastic  theories,  mathematical 
expediency  has  led  to  idealization  of  the  constitutive  rela- 
tions.  The  combination  of  springs  and  dashpots  selected  to 
characterize  the  stress-strain  properties  of  pavement  layer 
materials  depends  on  the  individual  investigator (s) . 
Basically,  the  criterion  is  to  obtain  a  better  curve  fit 
for  the  response  sought.   Usually,  some  difficulty  is 
involved  in  taking  the  inverse  Laplace  transform  to  obtain 
the  final  solution. 

Among  the  viscoelastic  models  commonly  used  are  the 
Maxwell  [39],  Voigt-Kelvin  [25],  van  der  Poel  (or  Standard 
Solid)  [40],  four-element  [21,  41]  or  combinations  of  these 
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models  [22,  31,  33].   The  model  constants  are  usually  evalu- 
ated by  running  a  creep  test.   The  response  to  a  unit  step 
load  is  then  approximated  by  the  Dirichlet  series 

n     -tE. 
i  =  l  x 

where  C.  and  E.  are  constants,  t  is  time  and  n  is  the  number 
of  terms  desired.   The  exponents  E.  are  given  specific  values 
as  suggested  by  Shapery  [42]  and  the  coefficients  C-  are 
determined  by  the  linear  least  squares  method  of  curve  fit. 
It  is  noticed  here  that  the  material  parameters  C-  are  pre- 
determined by  the  experimenter  so  as  to  avoid  the  non-linear 
set  of  equations  that  result  if  the  error  sum  of  squares  is 
minimized  with  respect  to  each  parameter  [42] . 

One  of  the  objectives  of  the  present  investigation  is 
to  determine  the  material  properties  defined  by  the  response 
to  a  unit  impulse  load  using  a  non-linear  least  squares 
curve  fitting  technique.   A  computer  subroutine  [43]  is 
employed  for  this  purpose. 

Use  of  Dynamic  Models 
Recently,  research  workers  have  recognized  the  fact  that 
pavements  should  be  regarded  as  dynamic  systems.   The  elastic 
and  viscoelastic  theories  do  predict,  within  a  limited  degree 
of  accuracy,  static  and  quasi-static  stresses  and  displace- 
ments, but  little  is  known  about  the  dynamic  behavior  of 
pavements.   Lattes,  Lions  and  Bonitzer  [44]  attempted  to 
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include  inertia  effects  but  the  method  was  in  the  develop- 
ment stage  and  no  numerical  results  were  published. 

In  dynamic  models,  the  mass,  dashpot  and  spring  function 
represent,  respectively,  inertia,  damping  effect  and  restor- 
ing force,  of  a  mechanical  system.   The  two  configurations 
mostly  used  to  represent  a  road  structure  are  shown  in 
Figures  2-a  and  2-b,  respectively,  for  a  single  degree  of 
freedom  (three-parameter)  model  and  a  seven-parameter  model. 
Using  sustained  vibrations  and  assuming  the  mechanical  system 
of  Figure  2-a  to  represent  a  road  system,  Heukelom  [45] 
analyzed  dynamic  deflections  of  pavements.   Harr  [46]  used 
a  three-parameter  model  to  study  the  effect  of  vehicle  speed 
on  dynamic  pavement  deflections  due  to  a  rectangular  pulse 
loading  function.   In  conjunction  with  their  study  to  calcu- 
late the  transient  responses  of  a  pavement  system,  Szendrei 
and  Freeme  [47]  experimentally  determined  the  shape  of  the 
impulse  to  which  a  road  is  subjected  under  traffic  loading. 
The  two  pulse  shapes  which  were  observed  are  illustrated  in 
Figure  3.   Simple  as  well  as  complicated  models  [47,  48,  49, 
50]  have  been  proposed  to  investigate  the  dynamic  response 
of  pavement  systems.   The  inadequacy  and  shortcomings  of  such 
models  have  been  discussed  in  the  literature  [45,  48,  49,  50]. 
The  fact  that  these  models  do  not  possess  sufficient  predic- 
tive capability  to  satisfy  either  existing  or  contemplated 
needs  is  not  surprising  since  the  representation  of  pavement 
behavior  by  a  spring-mass-dashpot  system  is  only  a  simple 
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b.     Seven-  Parameter   Model 


FIGURE  2      MECHANICAL     MODELS    COMMONLY     USEO 
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Time 
a.    Rectangular  Pulse    Near    the    Tire    Center 
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Time 
b.   Half -Sine    Pulse    Near    the    Tire    Edge 


FIGURE   3 


TYPICAL    PULSE    SHAPES    PRODUCED    BY 
A   MOVING   WHEEL    LOAD  (47). 
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idealized  model  in  which  distributed  effects  are  approximated 
by  lumped  parameters.   This  is  attributed  to  lack  of  a  true 
model  that  would  satisfactorily  serve  the  purpose. 

Transfer  Function  Approach 
Based  on  transfer  functions,  several  solutions  to 
mechanical  and  electrical  engineering  problems  are  found  in 
the  literature  [51,  52,  53,  54].   The  transfer  function 
represents  the  operator  which,  acting  on  the  time -dependent 
input,  yields  the  output  of  a  dynamic  system.   When  obtain- 
ing the  transfer  function  analytically,  the  system  param- 
eters are  defined  by  assumed  differential  equations.   For 
example,  the  equation  of  motion  for  the  system  shown  in 
Figure  2-a  is  given  by 


ray(t)  +  cy(t)  +  ky(t)  =  F(t) 


(1) 


where    m  =  mass, 

c  =  viscosity  coefficient, 
k  =  spring  constant, 
y(t)  =  output  displacement 

y  =  velocity, 
y(t)  =  acceleration, 
F(t)  =  input  forcing  function,  and 
t  =  time. 
The  transfer  function  G(s)  for  this  three-parameter  system 
is  obtained  from  the  ratio  of  the  Laplace  transform  of  the 
output  y(s)  to  that  of  the  input  F(s)  as 
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G(s)  -  ^Sl  (2D 

F(s) 


where  s  is  a  complex  variable. 

Neglecting  the  initial  conditions,  G(s)  becomes 


G(s) 


ms   +  cs  +  k 


A 
(s-Pj) (s-p2) 


(3) 


where     A  -  —  .  and  (4a) 

m  '  ^  J 


Pi»P»  =  '  fa  ±  /ts  -   i  (4b) 


Similarly,  the  differential  equations  of  motion  for  the 
seven-parameter  model  (Figure  2-b)  are  written  as  [55] 

miYiCt)  +  ciyi(t)  +  (k1+k3)y1(t)    k3y2(t)  =  F(t)   (5a) 

m2y2(t)  +  c2y2(t)  +  (k2+k3)y2(t)  -  k3yx(t)  -  0      (5b) 

and  the  transfer  function  becomes 

G(s)  =  i±tel 

F(s) 


_  A(s-ri)  (s-r2) 

(s-pi)  (s-p2)  (s-p3)  (s-pO 

where  the  constant  A,  the  zeros  r.  and  the  poles  p.  are 
functions  of  mi,  m2 ,  Ci,  c2,  ki ,  k2  and  k3. 


(6) 


In  general,  the  mathematical  model  from  any  mass-spring- 
dashpot  mechanical  system  can  be  written  in  the  form  [52] 

(s-rjfs-rj  •••  (s-r  ) 


m 


G(S)  =  7 T7 "i ? =r  (7) 

(s-Pl)(s-p2)  •••  (s-pn)  l/J 

However,  the  assumed  mechanical  model  specifies  the  number 
of  material  parameters  obtained  in  the  mathematical  repre- 
sentation . 

The  procedure  used  by  Swami  [56]  to  obtain  the  transfer 
function  for  asphaltic  concrete  presents  a  method  where  no 
mechanical  model  is  assumed.   The  technique  is  based  on 
direct  frequency  tests  which  require  steady  state  conditions 
between  the  input  and  output  when  both  vary  sinusoidally 
with  time.   The  amplitude  and  phase  relationship  between 
input  and  output  sine  waves  are  measured  and  the  process  is 
repeated  throughout  the  frequency  range  of  interest.   The 
ratio  of  the  magnitude  of  the  output  (displacement)  to  that 
of  the  corresponding  input  (force)  is  then  plotted  (in  deci- 
bels) against  the  log  of  frequency  to  obtain  the  frequency 
spectrum.   The  asymptotic  approximation  of  the  frequency 
spectrum  gave  the  following  transfer  function  of  the  system: 

r(*-\    -     A(s-r],)  (s-r2)  (s-r3)  f 

G(S)  "  (s-Pa)(s-p2)(s-p3)^-pJ  W 

The  study  [56,  57]  concluded  that  an  asphaltic  concrete  at 
a  constant  temperature  could  be  uniquely  represented  by  a 
transfer  function;  and  that  the  coefficients  in  the  transfer 
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function  were  believed  to  be  better  indicators  of  material 
performance  than  those  previously  used. 

The  frequency  response  method  of  determining  the  trans- 
fer function  may  be  impractical  in  cases  where  short  testing 
time  is  necessary.   Some  established  pavement  systems,  such 
as  Interstate  Highways,  may  not  tolerate  off-schedule  opera- 
tion for  steady-state  sinusoidal  measurements.   In  the 
present  investigation  the  impulse  response  of  a  pavement 
system  is  obtained  in  the  time  domain  from  a  single  impulse 
test  with  a  duration  of  less  than  one  second. 

Laboratory  Studies  of  Flexible  Pavements 
Literature  was  reviewed  to  provide  background  in  select- 
ing the  pertinent  factors  and  assigning  typical  values  to 
the  variables  used  in  this  study.   Results  of  the  laboratory 
studies  of  others  are  pertinent  in  this  regard. 

Subbaraju  [58]  conducted  model  studies  of  stresses  in 
the  upper  layers  of  pavements.   A  2\   -  inch  thick  asphaltic 
concrete  slab  was  placed  on  a  subgrade  soil  confined  in  a 
wooden  box  with  inside  dimensions  23|  by  25^  by  23  inches 
deep.   Static  loads  of  650  pounds  were  applied  through  a 
3$  -  inch  diameter  bearing  plate.   Assuming  the  stresses  for 
the  loading  used  would  not  extend  more  than  ten  inches  beyond 
the  center  loading  point,  no  edge  effect  was  considered.   The 
assumption  was  verified  by  his  test  results. 
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Monismith  and  Secor  [41]  attempted  to  validate  visco- 
elastic  theory  by  analyzing  a  thin  viscoelastic  plate  on  a 
Winkler  foundation.  Theoretical  deflections  were  compared 
with  those  measured  from  a  40  by  40  by  l\  inches  deep 
asphaltic  concrete  slab  resting  on  a  bed  of  closely  spaced 
springs.  In  some  instances,  discrepancies,  which  were 
larger  at  higher  test  temperatures,  were  attributed  to  the 
assumptions  made  in  the  analysis. 

To  validate  the  discrete-element  solutions  for  pavement 
slabs,  Agarwol  and  Hudson  [59]  tested  a  simple  two-layer 
model.   The  system  used  in  the  investigation  consisted  of 
a  small  aluminum  slab  seated  on  a  saturated  clay  subgrade 
in  a  2  by  2  by  2-foot  wooden  box.   Load  magnitudes  of  100 
and  200  pounds  were  transmitted  statically  through  a  \  -   inch 
diameter  rod. 

Model  studies  where  all  layer  materials  were  ideally 
simulating  [60]  or  actually  representative  [61,  62]  have 
been  conducted.   Surface  course  horizontal  dimensions  rang- 
ing from  9  by  9  inches  [59]  to  10  by  10  feet  [62]  have  been 
used. 

In  comparison  with  static  tests,  dynamic  tests  have 
been  relatively  few  in  number.   Brown  and  Pell  [63]  tested 
a  two-layer  system  subjected  to  a  pulse  load  having  an 
amplitude  of  approximately  7.6  kips  and  durations  of  0.1  and 
2  seconds.   The  study  confirmed  that  the  assumptions  of 
elastic  theory  are  not  generally  valid  for  pavement  materials 
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[63].   Waterhouse  [61]  tested  a  three-layer  test  pit  to 
failure  under  a  repeated  load  of  5  tons  magnitude  and  fre- 
quency of  8  cycles  per  minute.   The  sides  of  the  box  con- 
taining the  6-foot  square  by  5-foot  deep  specimen  were 
believed  to  have  only  a  slight  effect  upon  the  results  [61]. 

Drennon  and  Kenis  [62]  considered  pavements  as  linear 
systems  and  applied  Duhamel's  superposition  integral  to 
predict  repetitive  load  displacements  from  static  loads. 
The  maximum  load  intensity  applied  to  the  full-scale  asphaltic 
concrete  pavement  model  was  80  psi  at  a  frequency  of  one  cycle 
per  second.   It  was  concluded  from  the  experimental  results 
that  the  total  observed  deflections  increased  linearly  with 
load  and  that  the  dynamic  components  of  both  deflection  and 
strain  were  basically  linear.   Some  discrepancy  existed 
between  measured  and  predicted  displacements  partly  due  to 
using  only  the  creep  curve  rather  than  both  creep  and  recovery 
curves  as  input  to  the  superposition  integral  [62] . 

A  promising  field  test  to  evaluate  flexible  pavements 
by  impulsive  loads  was  carried  out  by  Isada  [49,  50].   Prom 
the  peak  displacement  and  the  magnitude  of  the  impulse,  a 
road  stiffness  was  determined  but  no  analysis  of  the  displace- 
ment time  data  was  performed. 

The  above  review  indicates  that  neither  the  classical 
theories  (elastic  or  viscoelastic)  nor  the  assumed  dynamic 
models  can  adequately  characterize  the  time-dependent 
response  of  a  pavement  system.   It  also  shows  that  the 
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transfer  function  theory  appears  to  be  a  sufficient  frame- 
work within  which  the  time-dependent  behavior  of  pavements 
can  be  investigated. 
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THEORETICAL  ANALYSIS 

The  basic  concepts  underlying  the  transfer  function 
theory  is  presented  in  this  section.   Three  areas  are  of 
major  interest:   general  assumptions,  experimental  determina 
tion  of  the  impulse  response  and  prediction  theory. 

Assumptions 

Transfer  function  theory  is  based  on  the  assumption 
that  pavements  behave  as  deterministic  systems.   A  relation- 
ship which  exists  between  known  input  and  output  conditions 
is  sought.   This  is  in  contradiction  to  classical  theories 
which  assume  this  relationship  to  be  known  initially.   The 
system  functions  are  schematically  shown  in  Figure  4.   Com- 
plex dynamic  situations  including  aircraft  flight  data  have 
been  synthesized  utilizing  the  assumption  that  a  determina- 
ble function  relates  the  output  to  the  input  [51,  52,  55]. 

Furthermore,  it  is  assumed  that  pavements  will  behave 
as  linear  systems.   Although  at  high  magnitudes  of  stresses 
and  temperatures  non-linearity  may  be  observed,  the  investi- 
gations reported  by  Heukelom  [45],  Drennon  and  Kenis  [62] 
and  others  [47,  48,  61]  have  shown  that  the  assumption  holds 
for  small  deflections  normally  experienced  by  pavement 
systems.   The  AASHO  Road  Test  results,  where  single  as  well 
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as  tandem  axle  loads  were  employed,  indicated  that  creep 
speed  deflections  were  linear  in  most  cases  [3,  pp.  102-103]. 

Experimental  Determination  of  the  Impulse  Response 
With  reference  to  Figure  4,  a  dynamic  system  can  be 
represented  by  a  transfer  function  G(s)  defined  as  the  ratio 
of  the  Laplace  transform  of  the  output  0(t)  to  the  Laplace 
transform  of  the  input  I(t).   By  this  definition, 

G(s)  =  5M 
I(s) 

°r  0(s)  =  G(s)  T(s)  (9) 

where  0(s)  =  L[0(t)]  (10a) 

T(s)  =  L[I(t)]  (10b) 

s  =  Laplace  transform  parameter, 

t  =  time  ,  and 

L  =  Laplace  transformation 
The  inverse  transform  of  G(s)  is 

G(t)  =  L-1[G(S)]  (11) 

where  L"1  is  the  inverse  Laplace  transformation. 
The  function  G(t)  is  determined  solely  by  the  system  param- 
eters and  is  called  the  impulse  response  or  the  response 
function  of  the  system.   The  operational  form  of  Equation  9 
can  be  expressed  in  the  time  domain  using  Borel's  theorem 
on  the  convolution  of  two  functions  [54]  to  obtain  the  output 
as 
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0(t) 


G(T)   I(t-T)dT 


(12) 


0 


G(t)  is  obtained  implicitly  by  the  finite  difference  method 
[64]  as 


0(tk)  -   ^  G(tj)  Ktk.j  +  i)  At 
(tk) J  KF   At 


(13) 


where   G(tn)  =  the  vector  quantity  of  the  response  function 

at  time  t,, 
0(t,)  =  the  value  of  the  output  at  t,  , 
1  (t,  )  =  the  value  of  the  input  at  t,  , 
K  =  conversion  factor, 
F   =  the  peak  value  of  the  input,  and 
At  =  time  interval. 
The  denominator  in  Equation  13  may  be  too  small  for 
Equation  13  to  converge.   If  this  is  the  case,  as  was  found 
to  be  generally  true  in  the  present  investigation,  K  is 
chosen  such  that  the  response  function  is  within  the  conver- 
sion region.   In  this  study,  a  denominator  value  of  9  or  more 
would  cause  Equation  13  to  converge.   The  K  factor  was  taken 
to  be  10  for  faster  convergence. 

The  experimental  basis  discussed  here  to  obtain  the 
response  function  is  the  impulse  test.   An  impulse  load 
input  of  known  peak  value  is  applied  for  a  very  short  time 
duration  to  the  system  whose  response  function  is  sought. 
Both  the  impulse  input  force  and  the  output  deflections  are 
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recorded  as  functions  of  time.   Figure  5  shows  schematic 
traces  of  the  input  and  the  output  necessary  to  determine 
such  a  response  function  in  conjunction  with  Equation  13. 
The  computation  results  in  vector  quantities  of  the  response 
function.   These  values,  shich  are  spaced  At  apart  in  the 
time  domain,  can  then  be  approximated  by  a  mathematical 
function  using  a  non-linear  least  square  curve  fit  method. 
The  closer  the  approximation  to  the  actual  response  function 
curve,  the  better  becomes  the  approximation  of  the  response 
function  sought.   The  function  used  in  this  investigation 
is 

G(t)  =  ae"3t  sin  yt  (14) 

where  G(t)  =  response  function, 
t  =  time ,  and 
a,  6  and  y  are  constants. 
As  a  check  on  the  response  functions  prior  to  curve 
fitting,  the  finite  difference  solution  of  Equation  12  was 
obtained  explicitly  [64]  as 

k 
0(tk)  =  J      GCtj)  I(tk.j+1)  At  (15) 

Thus,  the  output  calculated  from  Equation  15  should  match 
the  measured  output  used  to  obtain  the  response  function. 

The  modified  version  of  the  computer  program  [64]  used 
to  calculate  the  response  functions  and  to  predict  the  output 
is  shown  in  Appendix  A.   Subroutines  CONVLI  and  CONVLE  were 
employed  to  solve,  respectively,  Equation  13  and  Equation  15. 
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c.  Response  Function 
FIGURE  5       EXPERIMENTAL    DETERMINATION    OF 
SYSTEM    RESPONSE     FUNCTION 
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Prediction  Theory 
If  the  system  represented  by  Equation  9  (or  Equation  14) 
is  subjected  to  an  excitation  force  F(t)  and  the  correspond- 
ing output  time  function  is  y(t),  the  transformed  response 
is 

y"(s)  =  G(s)  F(s)  (16) 

By  the  convolution  theorem, 

y(t)  =  L-^GCs)  F(s)] 

=  G(t)*F(t) 

t 

G(t)  F(t-x)di  (17) 

0 

Equation  17  illustrates  that  one  can  calculate  the  system 
response  to  any  excitation  function  F(t)  if  the  response 
function  G(t)  is  known.   When  the  input  function  has  a 
mathematical  representation,  the  convolution  Equation  17 
can  be  solved  analytically.   The  solutions  for  a  step  load 
input  and  a  repeated  load  input  are  derived  below  for  a 
pavement  system  whose  response  function  is  G(t). 

Solution  for  Step  Loading 
Consider  the  step  function  input  F(t)  defined  as 

F(t)  =  0      for  t  <  0 


F(t)  =  FQ     for  t  >  0 


(18) 


where  F   is  the  magnitude  of  the  static  compressive  force. 
However,  due  to  the  introduction  of  the  conversion  factor  K 
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in  obtaining  the  response  functions,  Equation  18  has  to  be 
modified  as 


F(t)  =  0 
F(t)  =  KF( 
F(t)  -  Fo 


for  t  <  0 

for  0  <  t  <    t 

for  t  >  t 
n 


n 


(19) 


where  t   is  the  time  at  which  the  input,  from  which  the 

n  r 

response  function  is  obtained,  first  becomes  non-zero. 
Substituting  Equations  14  and  19  into  Equation  17, 


t 
r   n 


y(t)  =  kf, 


ae 


Bt       . 


6t 


sin  yt   di    +   Fn  ae  sin   yx    dx 


■   Fo   * 


KY 

2     ^    .,2 


CK-1) 


-et 


n 


y  +  r     zp  +  Y2 


n 


sin    (Ytn+<f) 


■et 


/62    +    Y2 


sin    (Yt+4>) 


(20) 


For    t   <   t    , 

n ' 


y(t)    =   KFQ  ae"BT   sin 


YT     ClT 


■   KF~   a 
o 


JL 


gt 


L62  +  y2     /% 


sin    (Yt  +  <j>) 


2        +        y  2 


where   <|>   =    tan_1(Y/&) 

At  t  =  0 ,  Equation  20a  yields: 


(20a) 
(21) 


y(0)  =  0 


(22) 


As  t  becomes  large,  the  steady-state  response  is  obtained 
from  Equation  20  as 
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y(t)  =  FQ  a 


32  +  Y2' 


(JLU 


3t 


n 


T  +  Y^ 


sin  (Ytn+4») 


(23) 


The  computer  program  STALOD  to  solve  Equation  23  is  given 
in  Appendix  B. 

Solution  for  Repeated  Loading 
With  the  modification  discussed  in  the  previous  section 
taken  into  account,  a  haversine  repeated  load  configuration 
shown  in  Figure  6  is  defined  as 


F(t)  =  0 


for  t  <  0 


KF, 


F(t)  =  -j^-  (1  -  cos  wt)      for  0  <  t  <  t 


n 


F(t)  =  -y-  (1  -  cos  cot)      for  t  >  t 


n 


(24) 


where  w  =  2tt/p,  in  which  p  =  period  of  loading. 


Substituting  Equations  14  and  24  into  Equation  17, 

t 

Br 


y(t)  = 


KFn  '  n 


ote 


sin  yt  (1-cos  a)(t-T))dx 


F   rt 
ro 


ae 


3t 


sin  yi    (1-cos  a)(t-T))dx 


n 


Ky_ 


(K-D   "3tn 


Y' 


/|F  +  Y: 


sin  (Ytn+<|>) 


6t 


/32  +  Y2 

F0a 

— t—  cos  wt- 


sin  (Yt+<t>) 


KT (Y-fQ +    Cy+uQ 

U2  + 


(Y-oo)2    62  +  (y+w)2J 


Fit) 


N 


period^ 
—  t=(2N-l)-| 


t°Np 


a.  Input:    F(  t)  «  -§>  ( l-cos^t) 
2  p 


y(t) 


32 


Time,t 


r\r\<?[r\wr\r\ 

I  Tir 


Time,  t 


b.  Output 


FIGURE   6      HAVERSINE     INPUT  AND    OUTPUT 
FUNCTIONS 
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-     (K-l) 


■Bt 


n 


sin 


L/V  +  (y-w)2 

-Bt, 


(Y-o))tn    +    $A 


11 


"V  +  (Y+w)2 

e-gt 
^B2  +  (y-o))2 


sm 


sin 


sin 


F0a 


/B2   +  (y+w) 2 

sin   oat    •(  K 

Bt 
(K-l) 


(Y  +  co)tn    +    <j>2 


(Y-aj)t    +    4>1 
(Y  +  oj)t    +    <pz 


B  B 

_B2  +  Cy-w)2        32  +  (Y+u))  2 


_/B2  +  Cy-w)  2 
-Bt, 


cos 


(Y-u)tn   +    <j>j 


n 


/B2  +    (y  +  w) 2 

e-gt 
/B2  +  (y-^)2 

e-gt 
/B2  +  Cy+oj)  2 


cos 


COS 


COS 


(Y  +  oj)tn    +    <J>2 


(Y-oj)t    +    4>  a 
(Y  +  aj)t    +    (J)2 


where      <J>    =    tan"1     (y/B)  , 


<f>x    =    tan"1 
<J>2    =    tan" x 


Y-oj 

I    B    J 
Y  +  oj 

[    B    J 

,  and 


(25) 


(26) 


Total  Response .   The  total  response  yT(t)  of  the  system 
is  evaluated  when  a  load  cycle  is  at  its  peak  position 
(Figure  6) ,  that  is  when 
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t  =  (2N-1)  | 


where  t  =  time, 

N  =  number  of  load  applications,  and 
p  =  period  of  loading. 


Since    sin    wt    =    sin 
and        cos    cot    =    cos 


^    (2N-1)    E 
P  2 


0, 


2^ 


(2N-1)    | 


1, 


Equation    26    reduces    to 

v    ft1    =    !°!   J     „ky  K    f     y    -    m  y    +    co 

XT^J  2         B2  +  y2         2    [3i+  (Y-w)2         62  +   (Y+a))2 


(K-l)         "6tn      .       ,    ^ 

sm    (ytn    +    <j>) 


/62  +  Y2 


K-l 


3t, 


_/B2  +  (y-w)2 


sin 


(Y-w)tn    +    tfrj 


n 


/B2  +  (y+w)2 

„-3t 


sm 


(Y  +  w)tn    +    <{>2 


/p  +  y2 


sin    (yt    +    (J>) 


-6t 


_/B2   +  (y-oj)  2 
"St 


sm 


(Y-w)t    +    <i>1 


/B2  +   (y+w) 2 


sin 


(Y  +  co)t    +    <f>2 


(27) 


(28) 
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Permanent  or  Accumulative  Response.   The  permanent 

nse  y  ft)  is  measure 
p 

(Figure  6) ,  that  is  when 


response  y  (t)  is  measured  at  the  end  of  a  load  cycle 


t  =  Np 


(2tt 
Since      sin   cot      =    sin      — 

1  P 

and     cos  cot   =  cos   — 

I  P 

Equation  26  yields 

y  (t)  .  !£  I   *y  4  -  5 

xpL  ;     2   ja2  + Y2    2 


Np 
Np 


=  0 


1, 


Y  -  co     +   _Y_+_ co_ 
82  +  CY"u) 2    82  +  (y+u) 2 


(K-l)   "6tn   .   ,  .     x. 
sin  (Ytn  +  <j>) 


/B2  +  Y2 
'K-l 


6t 


n 


_/32  +  (y-w)  2 


sin 


(Y-w)tn  +  <|)  j 


-Bt 


n 


/B2  +  (y+w) 2 


sin 


(Y+co)tn  +  <})2 


/B2  +  Y2 


sin  (Yt  +  cj>) 


Bt 


sin 


_/B2  +  (y-w)2 
-Bt 


/B2  +  (y+w) 2 


sin 


(Y-oj)t  +  $1 
(Y+co)t  +  <|) 


(29) 


(30) 


Appendix  C  presents  the  computer  program  REPLOD  to  solve 
Equations  28  and  30. 
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EXPERIMENTAL  INVESTIGATION 

The  experimental  phase  of  the  study  was  undertaken  to 
investigate  the  theory  of  pavement  behavior  developed  in  the 
preceding  section.   The  scope  of  the  investigation,  materials 
selected,  design  and  preparation  of  the  pavement  models, 
instrumentation  and  test  procedures  are  described. 

Scope 
The  study  was  limited  to  controlled  laboratory  testing 
of  three-layer  flexible  model  pavements  constructed  of 
representative  paving  materials.   A  silty  sand  subgrade,  a 
crushed  aggregate  base,  and  an  asphalt  concrete  surface  were 
the  components  of  the  three-layer  systems.   Surface  course 
thicknesses  of  1  and  2  inches  were  used.   The  3^  -  inch  base 
course  thickness  was  estimated  by  an  approximate  design 
procedure  briefly  outlined  in  the  next  subsection.   The  use 
of  an  asphalt- treated  base  and  a  3-inch  surface  course  thick- 
ness originally  planned  to  be  considered  were  deleted  because 
the  first  case  amounts  to  increasing  the  surface  course 
thickness  over  the  untreated  base  [3,  60],  while  in  the 
latter  case,  only  the  response  of  the  top  layers  would  have 
been  measured. 
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Static  loads  and  two  dynamic  loads,  namely,  impulse  and 
cyclic  loadings  were  thought  to  simulate  the  loading  condi- 
tions which  occur  in  a  prototype  pavement.   Contact  pressures 
of  14.93,  29.86  and  59.72  psi  were  applied  over  a  circular 
loaded  area  of  4-inch  diameter  at  the  center  of  the  pavements 
These  stress  magnitudes  were  considered  representative  of 
those  encountered  in  highway  pavements  by  vehicle  wheels. 
The  loading  plate  used  was  chosen  to  be  small  relative  to 
the  pavement  surface,  but  yet  to  allow  a  sufficient  portion 
of  the  subgrade  to  be  stressed.   The  tests  were  carried  out 
under  isothermal  conditions  in  a  constant  temperature  room 
at  50,  75  and  100°F. 

The  model  pavement  with  2-inch  surfacing  was  also  sub- 
jected to  corner  loadings  to  study  the  edge  effects  that 
might  result  as  a  consequence  of  the  sides  of  the  box. 
Loading  modes  were  impulse  and  cyclic.   The  test  temperature 
was  50°F  for  the  corner  loadings  since  any  edge  effects 
would  be  more  pronounced,  if  any  were  to  be  significant,  at 
these  test  conditions.   Only  the  largest  stress  magnitude, 
namely,  59.72  psi   was  applied  over  the  4-inch  diameter  plate 
for  the  edge  loading  condition. 

Approximate  Design  of  the  Model  Pavements 
Influence  curves  developed  by  Nijboer  [65],  from  the 
three-layer  elastic  analysis  presented  by  Jones  [13],  in 
conjunction  with  the  equivalent  layer  theory  [66]  ,  were 
used  for  the  purpose  of  estimating  the  thickness  of  the 
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base  course.   A  brief  outline  of  the  procedure  follows: 

1.  From  data  given  by  Yoder  [6],  the  elastic  moduli 
E3  =  2,500  psi  and  E2  =  5,000  psi,  were  obtained  for  the 
subgrade  and  the  base  course,  respectively.   The  correspond- 
ing value  for  the  surface  mixture  was  Ex  =  40,000  psi 

[15,  57] . 

2.  With  a  limiting  deflection  of  0.025  in.  [67,  68], 
an  applied  pressure  of  80  psi ,  a  radius  of  contact  area  of 

2  in.,  the  influence  curve  [65,  Graph  IV]  yields  a  thickness 
of  2.7  in.  of  asphaltic  concrete  over  the  subgrade. 

3.  A  one-inch  surface  was  thought  reasonable. 

4.  The  remaining  1.7  in.  multiplied  by  the  approximate 

surface  course/base  course  equivalency  of  y%^-  =    2  gave, 

/  b2 

upon  rounding  to  the  nearest  half-inch,  3^  in.  of  base  course 

Materials  and  Preparation  of  Model  Pavements 
The  model  pavements  were  contained  in  a  wooden  box 
reinforced  with  steel  angles.   Since  the  interest  was  in 
testing  the  pavement  models  inside  a  constant  temperature 
room,  the  maximum  possible  size  of  the  box  adopted  had  the 
internal  dimensions  32£  by  32 |  by  23£  inches  deep.   To 
minimize  moisture  losses,  the  box  was  made  of  exterior  type 
plywood  and  the  interior  sides  and  bottom  of  the  box  were 
coated  with  spar  varnish  diluted  with  linseed  oil. 

The  next  three  sections  present  the  results  of  tests 
performed  on  each  of  the  materials,  followed  by  a  descrip- 
tion of  placing  the  respective  material  in  the  box. 
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Subgrade  Soil 

A  silty  sand  was  used  for  the  subgrade.   Table  1  sum- 
marizes the  physical  properties  and  Figure  7  shows  the  grain 
size  distribution  curve  of  the  sand.   Compaction  test 
results  are  given  in  Figure  8. 

The  soil  was  stored  in  air-tight  bags  at  its  natural 
moisture  content  of  19.1  percent.   A  sufficient  amount  of 
soil  to  make  a  compacted  lift  of  about  one  inch  was  then 
spread  out  in  large  trays .   The  clods  were  broken  by  hand 
and  the  soil  was  mixed  frequently  to  insure  uniform  drying. 
Periodically,  the  moisture  content  was  determined  using  a 
Speedy  Moisture  Tester,  shown  in  Figure  9,  which  was  pre- 
calibrated  for  this  particular  soil  (Figure  10) .   When  the 
moisture  content  was  about  0.5  percent  above  the  optimum, 
the  soil  was  placed  in  the  box  in  a  predetermined  amount 
necessary  per  lift.   The  soil  was  then  compacted  in  the 
box  using  an  air  hammer,  weighing  14  pounds,  with  a  3- inch 
diameter  base.   The  air  supply  was  set  at  75  psi.   Compac- 
tion at  the  edges  and  leveling  were  achieved  by  a  hand  tamp 
weighing  16  pounds  and  with  a  base  plate  8  inches  square. 
Figure  11  shows  the  compaction  process.   Marked  lines  on  the 
interior  of  the  box  sides  indicated  when  the  required  com- 
paction was  obtained  as  calculated  from  Equation  31: 

Ywet  =  IT  (31) 
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TABLE  1 
PHYSICAL  PROPERTIES  OF  THE  SUBGRADE  SOIL 


Classification,  Unified  System  ---------  SM 

Liquid  Limit,1  %    ---------------  -  29.3 

Plastic  Limit,2  %      ---------------  23.3 

Plasticity  Index,2  %    -------------  -  6.0 

Optimum  Moisture  Content,3  °^    ----------  13.7 

Optimum  Dry  Density,3  lbs. /ft3  ---------  114. 4 

Specific  Gravity1*   ---------------  2.73 

Soaked  CBR  (California  Bearing  Ratio),5  I      -    -    -  10. 0 


1  AASHO  T  89 

2  AASHO  T  9  0 

3  AASHO  T  9  9 
*  AASHO  T  100 
5  AASHO  T  193 
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FIGURE   8       MOISTURE- DENSITY    RELATIONSHIP    OF    THE 
SUBGRAOE   SOIL 
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FIGURE   9      THE    SPEEDY    MOISTURE    TESTER 
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FIGURE  10     CALIBRATION    OF    THE    SPEEDY    MOISTURE 
TESTER    FOR    THE    SUB6RA0E    SOIL 
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FIGURE    II      SUBGRADE     SOIL     UNDER      COMPACTION 
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where  Ywet  =  wet  density  of  soil  in  lbs. /ft3, 
W  =  weight  of  soil  layer  in  pounds, 
h  =  compacted  thickness  of  layer  in  inches. 
The  constant  C  was  calculated  from  the  geometry  of  the  box 
to  be  1 . 64. 

Before  placing  each  subsequent  layer,  the  smooth  surface 
of  the  compacted  soil  in  the  box  was  scarified  to  minimize 
compaction  planes.   When  a  compacted  thickness  of  about  4 
inches  was  obtained,  in-situ  density  was  checked  by  the  sand- 
cone  method  (AASHO  T  191) .   Samples  were  taken  at  each  in- 
place  density  test  for  water  content  determination. 

The  subgrade  compaction  was  completed  in  20  layers  to 
give  a  total  compacted  thickness  of  1 6 |  inches.   The  average 
dry  density  of  the  compacted  lifts  was  115.5  lbs.  per  cu.  ft. 
which  was  1011  of  the  standard  AASHO  T  99.   The  moisture 
content  averaged  13.91. 

Base  Course 

The  material  selected  for  the  base  course  was  a  crushed 
limestone  aggregate  conforming  with  Indiana  specification 
Size  No.  53  (Figure  12).   Tests  performed  according  to 
AASHO  T  99  and  AASHO  T  193  gave,  respectively,  an  optimum 
density  of  126.4  lbs. /ft3  at  an  optimum  moisture  content  of 
6.2%,  and  a  soaked  CBR  of  51.3  percent. 

The  aggregate  was  compacted  to  the  total  of  3^  inches 
thick,  in  three  lifts,  with  the  same  air  hammer  used  for  the 
subgrade.   However,  for  the  base  course,  the  circular  base 
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plate  of  the  hammer  was  replaced  with  a    3\  -   inch  square  one 
Leveling  of  the  surface  was  accomplished  by  a  vibrator. 
A  6-inch  by  3-inch  base  plate  attached  to  the  vibrator 
served  as  a  leveling  foot.   Both  the  air  hammer  and  the 
vibrator  were  operated  at  a  line  pressure  of  75  psi.   The 
base  course  was  compacted  to  1001  of  the  maximum  dry  den- 
sity (AASHO  T  99)  as  controlled  by  Equation  31. 

Surface  Course 

The  bituminous  mixture  used  for  the  surface  course  was 
prepared  from  crushed  limestone  and  sand  aggregate  blended 
to  conform  with  Indiana  type  B  surface  mixture  gradation. 
Sieve  analysis  and  origin  of  the  aggregate  are  given  in 
Table  2.   Figure  13  shows  the  gradation  range  for  Indiana 
Type  B  Surface  Mixtures  and  the  gradation  chosen  for  this 
investigation. 

A  60-70  penetration  asphalt  cement  having  the  proper- 
ties shown  in  Table  3  was  used  for  the  surface  course.   The 
weight  of  asphalt  cement  mixed  with  the  aggregate  was  6 
percent  by  weight  of  aggregate  in  accordance  with  Indiana 
specifications.*   Specimens  were  made  and  tested  according 
to  the  Hveem  method  of  mix  design  [69] .   The  Hveem  test 
results  are  summarized  below: 


*   In  the  specifications,  the  range  of  asphalt  content  for 
Type  B  Surface  Mixture  is  5.5  -  6.6%  by  weight  of  mix. 
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TABLE  2 
SIEVE  ANALYSIS  OF  AGGREGATES  FOR 
THE  SURFACE  COURSE  MIXTURE 


Sieve 
Passing 

Size 
Retained 

Fraction 
% 

Material 

1/2 

in . 

3/8 

in. 

14 

Crushed  Delphi  Limestone 

3/8 

in . 

No. 

4 

36 

it       ii        it 

No. 

4 

No. 

8 

11 

H       ti       ii 

No. 

8 

No. 

16 

12 

Lafayette  River  Sand 

No. 

16 

No. 

30 

12 

ii        it      it 

No. 

30 

No. 

50 

7 

ii        ii      ii 

No. 

50 

No. 

100 

4 

ii        n      ii 

No. 

100 

No. 

200 

2 

ii        it      ii 

No. 

200 

2 

Greencastle  Limestone  Filler 
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TABLE  3 
RESULTS  OF  TESTS  ON  ASPHALT  CEMENT 


Penetration,  100  grams,  5  sec,  170?1      -----   62 

Penetration,  200  grams,  60  sec,  32°F1  -----   16 

Loss  on  Heating,  50  grams,  5  hr.,  325°F,2  %    -    -    -         0.02 
Penetration  of  Residue,  '  %    of  Original  -  -  -  -   83 

Softening  Point,  Ring  and  Ball,3   °F  -----  -  123 

Specific  Gravity  at  77°F  4  ----------  -    1.011 

Flash  Point,  Cleveland  Open  Cup,5  °F  ------  581 

Ductility  at  77°F,  5  cm/min.,6  cm   -------  150 

Kinematic  Viscosity  at  275°F,7  cSt  -------  460 

Kinematic  Viscosity  at  325°F,7  cSt  -------  137 


+ 


1 

ASTM 

D 

5 

2 

ASTM 

D 

6 

3 

ASTM 

D 

36 

it 

ASTM 

D 

70 

5 

ASTM 

D 

92 

E 

ASTM 

D 

113 

7 

ASTM 

D 

2170 
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Specimen 

Unit  Wei| 
(lbs. /ft 

3) 

Hveem 
Stabili 

ty 

A 

ir  Voids 

(*0 

A60-1 

150.72 

36.8 

3.66 

A60-2 

151.01 

57.8 

3.85 

A60-3 

149.82 

49.7 

5.07 

Average 

150.52 

48.1 

4.19 

Prior  to  preparing  the  bituminous  mixture  to  be  used 
as  surface  course  for  the  pavement  models,  compaction  and 
leveling  experimentation  was  conducted  with  the  air  hammer 
and  vibrator.   Although  the  devices  proved  to  be  satisfactory 
in  placing  the  base  course,  it  was  thought  essential  to 
correlate,  by  trial,  the  time  of  compaction  and  leveling 
with  the  compactive  effort  required  to  obtain  the  design 
density  at  the  selected  asphalt  content. 

The  mixture  for  the  surface  course  was  prepared  by 
weighing  5000-gram  batches  of  total  aggregate,  heating  the 
aggregate  to  325°F  and  then  mixing  it  with  the  asphalt, 
heated  to  the  same  temperature,  in  a  mechanical  mixer  for 
two  minutes.   After  mixing  was  completed,  each  batch  was 
put  in  a  flat  pan  and  placed  in  a  forced-air  draft  oven  to 
cure  for  15  hours  at  140°F. 

To  minimize  the  temperature  loss  from  the  mixture  due 
to  air  drifts  from  the  air  hammer  during  compaction,  the 
constant  temperature  room  was  set  to  the  maximum  operating 
temperature  of  113°F.   Compaction  of  the  bituminous  mixture 
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was  started  after  the  mixture  was  reheated  to  325°F.   The 
hot  mixture  was  spread  uniformly  over  the  base  course  in 
one  layer  such  that  the  compacted  thickness  would  be  one 
inch.   The  sides  of  the  box  served  as  the  form  for  laying 
and  compacting  the  mixture.   The  bituminous  mixture  was 
compacted  and  leveled  by  the  heated  air  hammer  and  vibrator 
for  the  predetermined  time.   The  density  obtained  for  the 
one-inch  asphaltic  concrete  surfacing  was  149.8  pcf.   The 
in-place  void  content  corresponding  to  this  density  was 
calculated  to  be  4.6  percent. 

Upon  completion  of  the  tests  with  the  model  pavement 
having  a  one-inch  surfacing,  more  batches  of  bituminous 
mixture  to  increase  the  surfacing  thickness  to  two  inches 
were  prepared  in  the  above  manner.   A  wooden  form,  with 
32j  by  32j  by  3|  inches  internal  dimensions  connected  to 
the  top  of  the  box  was  used  during  spreading  and  compaction 
of  the  mixture  over  the  one-inch  asphalt  surface.   Figure  14 
shows  the  bituminous  mixture  under  compaction.   Leveling  of 
the  surface  by  the  vibrator  is  shown  in  Figure  15.   The  same 
compactive  effort  was  applied  to  give  the  same  density  as 
for  the  one- inch  surface. 

Instrumentation 
Instrumentation  was  carried  out  with  consideration  for 
recording  the  time-dependent  input  and  output  parameters. 
A  crucial  part  of  the  instrumentation  technique  lies  in 
devising  a  system  which  is  capable  of  applying  an  impulse 
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FIGURE    14       BITUMINOUS      MIXTURE     UNDER 

COMPACTION 
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FIGURE     15     LEVELING    OF    THE    SURFACE 
COURSE 
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force  of  a  desired  magnitude  and  duration,  and  of  detecting 
the  surface  motion  at  any  required  location  on  the  pavement 
surface . 

The  equipment  consisted  of  a  loading  frame  designed  and 
fabricated  of  steel  members  coupled  with  an  MTS  electro- 
hydraulic  actuator.   Details  of  the  loading  frame  and  the 
rigid  support  for  the  model  pavements  are  shown  in  Figure  16. 
A  mechanical  holder  was  made  of  methyl  methacrylate  and 
connected  to  the  top  plate  (item  9  )  of  the  loading  frame. 
The  purpose  was  to  hold  the  servoram  while  the  MTS  console 
was  being  programmed  in  tension  prior  to  applying  a  single 
compressive  square  wave.   The  generation  of  an  impulse  by 
this  technique  is  schematically  illustrated  in  Figure  17. 

Loads  were  applied  by  the  MTS  machine  and  measured  with 
a  1000-lb.  capacity  Strainsert  load  cell,  type  FL1U,  mounted 
between  the  hydraulic  actuator  and  the  loading  plate.   The 
output  of  the  load  cell  was  recorded  with  a  Brush  Mark  280 
two-channel  recorder. 

Five  Sanborn  Linearsyn  linear  variable  differential 
transformers  (LVDT's)  were  used  to  sense  the  deflections. 
For  the  center  loading,  one  LVDT ,  designated  as  no.  3,  was 
located  3.25  inches  from  the  center  of  the  loading  plate. 
The  other  four  were  spaced  2.5  inches  apart  along  the  same 
radial  line  as  LVDT  no.  3.   Figure  18  shows  the  locations 
of  the  LVDT's  during  corner  loading. 


57 


FIGURE    16       LOADING      FRAME     AND    MODEL    PAVEMENT    SUPPORT 
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Changes  in  output  voltage  from  the  LVDT ' s  were  ampli- 
fied and  recorded  by  an  8-channel  Sanborn  recorder.   Cali- 
bration of  the  LVDT ' s  was  effected  by  displacing  the  core 
of  the  LVDT  with  a  micrometer  and  relating  this  displacement 
to  the  magnitude  of  the  pen  deflection  of  the  recorder.   A 
typical  calibration  sheet  for  the  LVDT ' s  is  shown  in  Table  4. 
One  dial  gage  located  at  the  same  radial  distance  as  LVDT 
no.  2  from  the  center  of  the  load,  but  as  close  as  possible 
to  the  LVDT,  served  as  a  check  during  static  and  cyclic 
testing . 

The  LVDT ' s  were  fixed  to  the  pavement  surface  by  screw- 
ing one  end  of  a  brass  extension  rod  into  the  LVDT  core  and 
the  other  end  into  a    I  -   inch  by  \   -  inch  by  \   -  inch  thick 
brass  plate.   The  plate  was  in  turn  glued  to  the  appropriate 
location  on  the  pavement  surface  using  Eastman  9]0  Adhesive. 
The  LVDT ' s  were  supported  by  aluminum  channels  and  plexiglas 
beams  connected  to  the  top  angles  of  the  box.   Figure  19 
shows  the  general  view  of  the  test  set-up  inside  the  constant 
temperature  room  for  the  center  loading.   A  close  view  of 
the  LVDT ' s  is  shown  in  Figure  20. 

Test  Procedures 
Three  tests  are  described  in  this  section:   impulse, 
static  and  repeated  tests.   The  coding  shown  in  Figure  21 
was  used  to  make  data  identification  orderly. 
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TABLE  4 
TYPICAL  FORMAT  FOR  LVDT  CALIBRATION 


Temperature:  75°F 

Power  Amplifier  Setting:  X«l 

Preamplifier  Setting,  mv/cm:  20 

Pen  Displacement,  mm:  20 


Channel 


LVDT         LVDT  Core  Displacement, 
Designation  x  0.001  inch 


1  1  3.0 

2  2  3.2 

6  3  3.0 

7  4  3.2 

8  5  2.8 
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FIGURE   19       GENERAL  VIEW    OF   TEST    SETUP 
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FIGURE   20       CLOSE   VIEW    OF  TEST    SETUP 
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Series      I    L  I    I 


► 


Load    Intensity: 

14.93  psi   (I),    29.86  psi  (2),  or   59.72  (3) 

Input    Function: 

Impulse  (I),     Static  (S),  or    Repeated  (R) 

Test  Temperature: 

50°  F  (L),    75°F(M),  or    IOO°F  (H) 


Thickness   of    Surface   Course: 
I  inch  (I),  or  2  inches  (2) 


*  Applied   to  a   4-inch    diameter   plate,   these   pressures 
amount  to,  respectively,  187.5,  375.0,  and  750.0  lbs. 


FIGURE    21      COOING    OF    DATA 
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Impulse  Load  Tests 

With  the  function  generator  in  the  square-wave  mode  and 
by  suitable  manipulation  of  the  MTS  console,  impulse  load 
magnitudes  of  187.5,*  375.0  and  750  lbs.  were  applied  for  a 
duration  of  about  0.16  sec.   This  time  and  a  drop  height  of 
0.5  inch  were  determined  from  preliminary  tests  during  con- 
ditioning of  the  model  pavements. 

Typical  traces  of  input  load  and  output  deflections  are 
shown  in  Figure  22.   A  format  for  data  reduction  is  illus- 
trated in  Table  5. 

Static  Load  Tests 

For  static  load  tests,  the  frequency  of  the  MTS  machine 
was  set  at  one  cycle  per  second,  the  function  generator  in 
the  ramp  mode  and  the  control  mode  in  the  manual  trigger. 
By  these  adjustments,  a  static  compressive  force  of  the 
programmed  magnitude  was  applied.   The  load  was  held  on 
until  the  deflections  became  almost  constant  with  time. 
This  was  found  to  be  in  the  order  of  a  few  minutes. 

Figure  23  shows  representative  graphical  data  for  this 
type  of  test. 

Repeated  Load  Tests 
In  advance  of  starting  a  cyclic  load  test,  the  function 
generator  of  the  MTS  machine  was  programmed  in  the  compressive 
haversine  form  and  in  the  control  mode.   Once  testing  started, 


*   For  Series  1HI1,  the  peak  load  was  found  to  be  180.0  lbs 
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FIGURE  22      TYPICAL    TRACES    OF   AN     IMPULSE 
LOAD    TEST 


TABLE  5 

TYPICAL  DATA  REDUCTION  FOR  IMPULSE  TEST 

Series  :   1HI 3 
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Time 
(sec .  ) 

Distance  £ 

rom  Loac 

Center,*  inch 

les 

0 

3. 

25 

5. 

75 

8 

25 

Pen 

Dispt. 

Load 

Pen 
Dispt. 

Deflec. 

Pen 
Dispt. 

Deflec. 

Pen 
Dispt. 

Deflec. 

mm 

lbs. 
xlOO** 

mm 

in.xio"* 
x  1 .  5  0** 

mm 

in.xio-4 
x0.80** 

mm 

in.xlO"" 
xO.40** 

0. 

0 

0 

0 

0 

0 

0 

0 

0.02 

0 

0 

0 

0 

0 

0 

0 

0 

0.04 

1.0 

100.0 

0.6 

0.90 

1.0 

0.80 

1.0 

0.40 

0.06 

6.0 

600.0 

10.0 

15.00 

8.0 

6.40 

5.6 

2.24 

0.08 

7.5 

750.0 

11.4 

17.10 

12.0 

9.60 

8.6 

3.44 

0.10 

7.5 

750.0 

13.0 

19.50 

13.0 

10.40 

9.7 

3.88 

0.12 

7.5 

750.0 

13.6 

20.40 

13.4 

10.72 

10.2 

4.08 

0.14 

7.5 

750.0 

14.0 

21.00 

13.8 

11.04 

10.3 

4.12 

0.16 

7.5 

750.0 

14.1 

21.15 

14.0 

11.  20 

10.4 

4.16 

0.18 

7.  5 

750.0 

14.2 

21.30 

14.4 

11.52 

10.6 

4.24 

0.20 

7.5 

750.0 

14.4 

21.60 

14.6 

11.68 

10.7 

4.28 

0.  22 

7.5 

750.0 

14.5 

21.75 

14.8 

11.84 

10.8 

4.32 

0.  24 

0 

0 

8.0 

12.00 

6.0 

4.80 

8.0 

3.20 

0.26 

0 

0 

3.0 

4.50 

3.5 

2.80 

4.0 

1.60 

0.  28 

0 

0 

1.6 

2.40 

2.6 

2.08 

2.8 

1.12 

0.30 

0 

0 

1.4 

2.10 

2.  2 

1.  76 

2.0 

0.80 

0.32 

0 

0 

1.2 

1.80 

2.0 

1.60 

1.9 

0.76 

0.34 

0 

0 

1.1 

1.65 

1.9 

1.52 

1.8 

0.72 

0.36 

0 

0 

1.0 

1.50 

1.8 

1.44 

1.6 

0.64 

0.38 

0 

0 

0.8 

1.20 

1.7 

1.36 

1.5 

0.60 

0.40 

0 

0 

0.7 

1.05 

1.5 

1.20 

1.4 

0.  56 

0.42 

0 

0 

0.6 

0.90 

1.4 

1.12 

1.3 

0.52 

0.44 

0 

0 

0.6 

0.90 

1.3 

1.04 

1.2 

0.48 

0.46 

0 

0 

0.6 

0.90 

1.  2 

0.96 

1.1 

0.44 

0.48 

0 

0 

0.6 

0.90 

1.1 

0.88 

1.1 

0.44 

0.50 

0 

0 

0.6 

0.90 

1.1 

0.88 

1.1 

0.44 

** 


The  deflection  functions  recorded  at  distances  10.75  and 
13.25  inches  were  too  small  to  justify  this  type  of  data 
reduction. 

Calibration  Factors,  lbs. /mm  or  in.  x  10"Vmm 
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FIGURE  23      TYPICAL    TRACES    OF    A     STATIC    LOAD    TEST 
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the  hydraulic  actuator  moved  up  and  down  from  zero  load  to 
maximum  and  back  to  zero.   This  sequence  was  repeated 
automatically  for  the  programmed  number  of  250  cycles  at 
a  rate  of  15  cycles  per  minute. 

Typical  traces  of  repeated  load  input  and  deflection 
outputs  are  illustrated  in  Figure  24.   Formats  for  data 
reduction  are  shown  in  Tables  6  and  7,    respectively,  for 
static  and  repeated  load  tests. 
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TABLE  6 
TYPICAL  DATA  REDUCTION  FOR  STATIC  TEST 
Series:   1HS3 


X 

Preamp . 
Setting 

Calib . 
Factor 

Pen 
Dispt. 

Total 
Deflec. 

in. 

mv/cm 

in.  xlO"  Vnim 

mm 

in.xio-1* 

3.25 

50 

4.000 

6.4 

25.60 

5.75 

20 

1.600 

9.0 

14.40 

8.25 

10 

0.800 

7.0 

5.60 

10.75 

5 

0.375 

2.0 

0.75 

13.25 

2 

0.140 

0.0 

0.0 

TABLE  7 
TYPICAL  DATA  REDUCTION  FOR  REPEATED  TEST 
Series:   1HR3 


X 

Preamp . 
Setting 

Calib. 
Factor 

Pen 
Dispt. 

Total 
Deflec. 

in. 

mv/cm 

in.  xlO"  Vmm 

mm 

in.xKT1* 

3.25 

50 

4.000 

6.2 

24.80 

5.75 

20 

1.600 

8.2 

13.12 

8.25 

10 

0.800 

6.1 

4.88 

10.75 

5 

0.375 

1.9 

0.71 

13.25 

2 

0.140 

0.0 

0.0 
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RESULTS  AND  DISCUSSION 

The  theory  and  experimental  techniques  to  evaluate  the 
response  behavior  of  a  pavement  system  were  advanced  in 
previous  sections.   The  test  data  and  analysis  of  results 
are  presented  in  this  Section  and  Appendix  A  together  with 
their  implications  and  the  indicated  relationships.   The 
scope  of  the  investigation  included  two  flexible  model  pave- 
ments, differing  in  the  thickness  of  the  surface  course, 
tested  at  three  stress  levels  and  three  temperatures. 
Impulse,  static  and  repeated  loads  were  applied  at  the 
center  of  the  pavements.   Corner  loadings  were  also  con- 
ducted to  investigate  the  edge  effects  that  might  result 
as  a  consequence  of  the  sides  of  the  box.   For  all  series, 
surface  deflections  were  recorded  as  functions  of  time  at 
five  locations.   The  results  are  presented  and  discussed  in 
the  following  sequence: 

1.   Impulse  Load  Test  Results 

a.  Normal  Distribution  and  the  Deflection  Basin 

b.  Response  Function  of  Flexible  Pavements 

c.  Form  of  the  Response  Function 

d.  Check  for  the  Response  Function 

e.  Load- Independency  of  the  Response  Function 

f.  Effect  of  Temperature  on  the  Response  Function 
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g.   Effects  of  Surface  Course  Thickness  and  Spatial 
Location 

2.  Static  Load  Results 

3.  Repeated  Load  Results 

4.  Corner  Loading  Test  Results 

Impulse  Load  Test  Results 
The  peak  load  magnitudes  employed  in  most  of  the 
impulse  tests  were  187.5,  375,0  and  750.0  lbs.   These  corre- 
spond, respectively,  to  applied  pressures  of  14.93,  29.86 
and  59.72  psi.   The  time  duration  was  0.16  second.   This 
was  established  by  adjusting  the  drop  height  of  the  hydraulic 
actuator  and  the  frequency  of  the  MTS  function  generator. 
The  capability  of  the  testing  equipment  was  a  limiting  factor 
in  using  'time  duration  less  than  0.16  second. 

Normal  Distribution  and  the  Deflection  Basin 
The  peak  values  from  the  deflection  functions  are  pre- 
sented in  Table  8  and  plotted  in  Figures  25  through  30. 
These  data  indicate  that  in  all  cases  the  deflections  in- 
crease with  increase  in  temperature.   This  was  to  be  expected 
since  asphaltic  materials  are  thermoplastic  and  an  increase 
in  temperature  will  result  in  increase  of  their  flow  charac- 
teristics . 

As  the  thickness  of  a  pavement  structure  is  increased, 
the  resistance  to  imposed  loads  increases  and  the  system 
tends  to  deform  less.   However,  this  behavior  was  not  uniform 
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TABLE  8 
PEAK  DEFLECTIONS  FROM  IMPULSE  TESTS 


Series 

y(x), 

inch  (x 

ID"*) 

Distance  from 

Load  Center,  inches 

3.25 

5.25 

8.25 

10.75 

13.25 

1LI1 

1.50 

0.64 

0.16 

_  _ 

-  - 

1LI2 

2.48 

0.77 

0.32 

1LI3 

7.65 

3.04 

0.67 

1MI1 

3.23 

1.52 

0.38 

1MI2 

7.50 

4.00 

0.86 

1MI3 

15.00 

8.00 

2.00 

1HI1 

4.05 

2.24 

0.78 

0.03 

1HI1* 

4.28 

2.24 

0.86 

0.08 

1HI2 

11.70 

6.16 

2.11 

0.15 

1HI3 

19.20 

10.40 

3.44 

0.36 

1HI3* 

21.75 

11.84 

4.32 

0.57 

2LI1 

0.83 

0.48 

0.16 

2LI1* 

0.75 

0.42 

0.18 

2LI2 

1.43 

0.83 

0.22 

2LI2* 

1.43 

0.86 

0.29 

2LI3 

4.88 

3.36 

1.36 

0.69 

0.20 

2MI1 

1.91 

1.40 

0.74 

0.38 

0.22 

2MI2 

3.75 

2.76 

1.44 

0.  75 

0.32 

2MI2* 

3.64 

2.52 

1.09 

0.68 

0.25 

2MI3 

8.10 

5.04 

2.48 

1.28 

0.49 

2MI3* 

8.40 

5.12 

2.52 

1.43 

0.53 

2HI1 

3.60 

2.40 

1.04 

0.68 

0.22 

2HI1* 

3.30 

2.40 

1.12 

0.68 

0.25 

2HI2 

7.20 

5.08 

2.64 

1.43 

0.78 

2HI2** 

6.90 

4.64 

2.16 

1.20 

0.48 

2HI2** 

6.90 

4.80 

2.24 

1.28 

0.48 

2HI3 

12.00 

9.92 

4.48 

2.70 

1.05 

*   Represent  duplicate  series  of  tests. 
**   Represent  triplicate  series  of  tests 
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regarding  the  spatial  locations  x  where  measurements  were 
made.   At  the  two  locations  closest  to  the  loaded  position, 
decreased  deflections  were  observed  with  increase  in  thick- 
ness, while  at  the  other  three  locations,  namely  8.25  inches 
and  beyond,  deflections  increased  with  increasing  thickness. 
See  the  typical  comparisons  made  in  Figure  31.   This  phenome- 
non was  attributed  to  the  slab  action  of  the  pavement  system. 
When  a  system  with  a  thicker  surface  is  subjected  to  load, 
it  will  demonstrate  less  flexibility  than  would  be  the  case 
for  a  thinner  pavement.   Consequently,  the  rate  of  decrease 
in  deflection  with  increasing  x  is  less  than  that  for  a 
thinner  surface. 

The  function  selected  to  fit  the  peak  deflection  data 
of  Table  8  was  of  the  form 

y(x)  =  y0e"Dx2  (32) 

where  y(x)  =  the  deflection  at  a  distance  x  from  the  load 

center , 
y0  =  the  maximum  deflection  of  the  deflected  basin, 

and 
D  =  a  constant,  reflecting  the  attenuation  of  the 

deflected  basin  with  x. 
The  values  of  y0  and  D  determined  by  a  non-linear  regression 
method  are  given  in  Table  9.   The  non-linear  regression 
analysis  gave  squared  correlation  coefficients  (R2),  as 
reported  in  Table  9.   These  range  from  a  little  less  than 
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TABLE  9 
VALUES  OF  yQ  AND  D  IN  THE  EQUATION  y(x)  =  y0e"Dx 


Series 

y0 

in.  (x  10"1*) 

D 
1/in.2 

R2 
% 

Number 

of 
Points 

1LI1 

2.256 

0.0385 

99.95 

5 

1LI2 

4.071 

0.0475 

99.26 

5 

1LI3 

11.881 

0.0416 

99.97 

5 

Average 

0.0425 

1MI1 

4.716 

0.0353 

99.85 

5 

1MI2 

10.723 

0.0325 

99.32 

5 

1MI3 

21.236 

0.0317 

99.49 

5 

Average 

0.0332 

1HI1 

5.676 

0.0288 

99.66 

10 

1HI2 

16.124 

0.0299 

99.83 

5 

1HI3 

27.907 

0.0287 

99.07 

10 

Average 

0.0291 

2LI1 

1.066 

0.0274 

98.88 

10 

2LI2 

1.959 

0.0281 

99.07 

10 

2LI3 

6.106 

0.0199 

99.38 

5 

Average 

0.0251 

2MI1 

2.254 

0.0153 

99.54 

5 

2MI2 

4.437 

0.0168 

98.87 

10 

2MI3 

9.997 

0.0195 

99.44 

10 

Average 

0.0172 

2HI1 

4.163 

0.0176 

98.85 

10 

2HI2 

8.388 

0.0172 

99.11 

15 

2HI3 

14.560 

0.0149 

97.88 

5 

Average 

0.0166 
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0.98  to  close  to  unity.   This  is  an  indication  that  the 
model  chosen  to  fit  the  data  of  deflected  profiles  is 
adequate  to  the  task. 

It  is  observed  from  Figures  25-31  that  the  portion  of 
the  pavement  surface  under  the  loading  plate  was  assumed  to 
have  experienced  a  constant  displacement  at  the  indicated 
test  conditions.   This  was  a  consequence  of  the  use  of  a 
rigid  loading  plate.   The  choice  of  the  rigid  plate  itself 
was  made  to  obtain  a  distinct  trace  of  load  duration  during 
impulse  tests. 

The  effect  of  temperature  on  y0  is  shown  in  Figures  32 
and  33,  respectively,  for  the  model  pavements  having  one- 
inch  and  two-inch  surface  courses.   As  temperature  increases, 
the  less  stiff  and  more  responsive  the  material  becomes. 
Thus,  for  both  model  pavements,  y0  increases  with  increases 
in  temperature.   However,  the  plots  of  y0  versus  temperature 
do  exhibit  differences  between  the  two  cases.   For  the  one- 
inch  surface,  the  bottom  layers,  whose  response  constitutes 
most  of  the  total  deflection,  are  less  affected  by  tempera- 
ture changes  than  the  bituminous  surface  course  and  hence 
the  rate  of  increase  in  y0  decreases  with  increasing  tempera- 
ture.  On  the  other  hand,  as  the  thickness  of  the  asphaltic 
concrete  surface  is  increased,  the  influence  of  the  surface 
layer  becomes  more  pronounced.   Inasmuch  as  the  bituminous 
mixture  is  more  temperature  susceptible  than  the  other  layer 
materials,  increasing  temperature  increases  the  rate  of 
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increase  in  y0.   This  is  also  considered  as  evidence  in  this 
investigation  that  the  subgrade  and  the  base  course,  in 
addition  to  the  surface  layer,  responded  to  the  input  loads. 

Figure  34  depicts  the  effect  of  surface  course  thick- 
ness on  y0  for  the  three  load  magnitudes  and  the  three  test 
temperatures.   As  the  temperature  increases  the  slopes  of 
the  lines  increase,  supporting  the  above  discussion  about 
the  effect  of  increasing  surface  course  thickness  on  y0. 

The  factor  D  in  Equation  32  is  a  measure  of  how  the 
deflected  basin  attenuates  with  spatial  distance  x.   Smaller 
D  values  reflect  more  spreading  of  the  deflection  basin. 
Table  9  indicates  that  the  D  parameter  does  not  change 
appreciably,  even  though  the  magnitude  of  the  peak  load  is 
increased  at  any  one  thickness  or  temperature.   Average 
values  of  D  versus  temperature  are  plotted  in  Figure  35  for 
the  two  thicknesses.   Inasmuch  as  increases  in  temperature 
increase   pavement  response  for  a  given  thickness,  attenua- 
tion of  the  deflection  basin  with  x  decreases  and  hence  D 
decreases.   The  figure  also  indicates  that  the  slab  action 
is  reflected  in  the  parameter  D.   As  the  surface  course 
thickness  is  increased,  the  deflected  basin  spreads  farther 
and  thus  D  decreases.   It  appears  from  Figure  35  that  there 
is  no  interaction  between  the  effects  of  temperature  and 
surface  course  thickness  on  D. 
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Response  Function  of  Flexible  Pavements 
The  response  functions  were  calculated  from  the  input 
load  and  the  output  deflection  functions  using  Equation  13. 
Computations  were  made  by  the  computer  utilizing  the  sub- 
routine CONVLI  of  the  program  listed  in  Appendix  A.   The 
results  of  the  calculations  are  also  shown  in  Appendix  A. 
For  each  test  series,  the  table  of  the  input  and  outputs  is 
followed  by  the  response  function  table.   The  three  param- 
eters are  given  as  functions  of  time.   The  third  table  in 
each  series  is  the  deflection  function  predicted  by  Equation 
15  and  as  computed  by  subroutine  CONVLE. 

Typical  plots  of  the  input  and  outputs  are  shown  in 
Figures  36,  37  and  38,  and  the  corresponding  response 
functions  are  illustrated  in  Figures  39,  40  and  41.   It  is 
seen  that  the  response  functions  reach  a  first  peak  and  then 
oscillate  with  time.   The  peak,  in  general,  occurs  in  the 
time  band  of  0.04  to  0.08  seconds.   The  general  shape  of  the 
response  functions  and  the  corresponding  results  agree  with 
those  obtained  from  mechanical  models  [45,  49,  50].   However, 
when  mechanical  models  are  employed  to  describe  pavement 
behavior,  the  parameters  must  be  defined  initially,  usually 
in  the  form  of  differential  equations.   Furthermore,  no 
reliable  method  exists  today  to  determine  these  parameters 
"ab  initio."   This  study  is  characterized  by  the  fact  that 
the  response  function  of  a  pavement  system  is  determined 
directly  from  the  input-output  data  of  impulse  tests.   The 
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consistency  of  the  response  functions  obtained  in  this 

study  appears  to  validate  the  first  assumption  that  "a 

determinable  relationship  exists  between  known  input  and 
output  of  a  pavement  system." 

Form  of  the  Response  Function 
Some  typical  response  functions  are  represented  in 
Figures  39,  40  and  41  for  the  indicated  test  conditions. 
In  the  interest  of  obtaining  more  information  about  pavement 
response  functions  and  to  provide  a  rational  interpretation 
of  the  results,  the  response  function  curves  were  approxi- 
mated by  the  mathematical  model 

G(t)  =  ae"6t  sin  yt  (14) 

where  a,  g  and  y   are  parameters  to  be  determined  by  a  non- 
linear regression  analysis.   The  results  of  the  analysis  are 
shown  in  Tables  10,  11  and  12.   It  is  seen  from  these  tables 
that  the  squared  correlation  coefficient  (R2)  ranged  from 
66  to  93  percent. 

From  the  technique  followed  in  approximating  the 
response  function  curve,  it  is  obvious  that  the  number  of 
parameters  required  to  describe  the  response  function  depends 
on  the  model  selected.   As  the  number  of  parameters  in  the 
model  increases,  the  accuracy  of  curve- fitting  improves  and, 
consequently,  the  value  of  R2  increases. 
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TABLE  10 
DESCRIPTIVE  PARAMETERS  IN  THE  RESPONSE  FUNCTIONS  AT  50°F 

G(t)  =  ae"et  sin  yt 


Thick-   x       F         a         6         y  R: 

p 
ness  r 

in.    in.       lbs. 


lb.-sec.       1         1 


x  10-6      sec.      sec 


3.25     375.0  6.86  10.764  13.382  85.73 

750.0  10.40  10.551  15.848  87.48 

Average  8.63  10.658  14.615 

5.75     375.0  1.36  7.127  14.593  78.09 

750.0  3.60  8.784  17.285  86.79 

Average  2.48  7.956  15.939 

8.25     375.0  0.48  6.140  14.470  75.86 

750.0  0.72  8.195  15.716  86.32 

Average  0.60  7.168  15.093 

3.25     187.5  5.39  12.879  12.336  92.39 

375.0  3.57  9.847  12.709  87.60 

750.0  5.39  8.785  16.018  84.34 

Average  4.78  10.504  13.688 

5.75     187.5  2.41  11.302  14.222  83.88 

375.0  1.70  7.517  12.988  83.16 

750.0  3.51  7.540  16.406  84.82 

Average  2.54  8.786  14.539 

8.25     187.5  0.70  9.610  12.598  65.68 

375.0  0.51  7.409  12.996  73.44 

750.0  1.56  8.792  16.824  83.86 

Average  0.92  8.604  14.139 
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TABLE  11 
DESCRIPTIVE  PARAMETERS  IN  THE  RESPONSE  FUNCTIONS  AT  75°F 

G(t)  =  ae"et  sin  yt 


Thick- 

X 

F 
P 

lbs. 

a 

e 

Y 

R2 

ness 
in. 

in. 

in. 

1 
sec 

1 
sec. 

Ib.-sec. 
x  10~6 

% 

1 

3.25 

187.5 

21.28 

12.604 

13.875 

88.  26 

375.0 

21.89 

11.537 

14.989 

86.00 

750.0 

24.61 

12.348 

13.752 

87.44 

Average 

22.59 

12.163 

14.205 

5.75 

187.5 

10.31 

13.079 

12.791 

85.77 

375.0 

13.73 

12.908 

13.496 

88.41 

750.0 

14.42 

13.084 

12.478 

89.07 

Average 

12.82 

13.024 

12.922 

8.25 

187.5 

2.08 

10.434 

11.703 

87.44 

375.0 

2.52 

11.177 

15.026 

87.98 

750.0 

3.51 

13.199 

13.934 

87.35 

Average 

2.70 

11.603 

13.554 

2 

3.25 

187.5 

9.49 

9.500 

13.868 

90.80 

375.0 

7.98 

8.  232 

12.914 

88.23 

750.0 

12.26 

11.000 

15.272 

88.52 

Average 

9.91 

9.577 

14.351 

5.75 

187.5 

6.9  7 

9.525 

14.191 

90.31 

375.0 

4.97 

7.164 

14.398 

81.43 

750.0 

9.17 

13.347 

14.020 

89.63 

Average 

7.04 

10.012 

14.203 

8.25 

187.5 

4.21 

10.523 

13.452 

92.32 

375.0 

3.53 

9.195 

14.479 

91.03 

750.0 

3.  73 

11.471 

15.260 

87.63 



Average 

3.82 

10.396 

14.397 

100 


TABLE  12 
DESCRIPTIVE  PARAMETERS  IN  THE  RESPONSE  FUNCTIONS  AT  100°F 

G(t)  =  ae"6t  sin  yt 


Thick- 

X 

F 

a 

0 

Y 

R2 

ness 
in. 

in. 

P 
lbs. 

in. 

1 
sec . 

1 
sec. 

% 

lb. -sec. 
x  lO"6 

1 

3.25 

180.0 

19.33 

8.953 

16.968 

84.62 

375.0 

26.03 

7.950 

15.353 

83.28 

750.0 

26.30 

8.633 

15.334 

84.30 

Average 

23.89 

8.512 

15.885 

5.75 

180.0 

12.35 

10.044 

15.383 

88.75 

375.0 

22.13 

12.506 

12.335 

82.95 

750.0 

14.93 

9.478 

14.911 

81.78 

Average 

16.47 

10.676 

14.210 

8.25 

180.0 

4.41 

10.523 

15.061 

88.57 

375.0 

6.15 

10.951 

14.452 

78.24 

750.0 

6.46 

10.659 

12.977 

85.20 

Average 

5.67 

10.711 

14.163 

2 

3.25 

187.5 

17.54 

11.273 

14.046 

83.11 

375.0 

23.25 

13.187 

13.870 

86.90 

750.0 

19.33 

12.960 

12.700 

80.87 

Average 

20.04 

12.473 

13.539 

5.75 

187.5 

10.89 

9.051 

14.274 

87.87 

375.0 

11.08 

8.779 

12.825 

89.93 

750.0 

15.66 

12.068 

13.983 

89.15 

Average 

12.  54 

9.966 

13.694 

8.25 

187.5 

4.56 

8.202 

13.798 

84.81 

375.0 

7.46 

11.846 

11.680 

92.34 

750.0 

6.68 

11.081 

13.015 

93.19 

Average 

6.23 

10.376 

12.831 
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Check  for  the  Response  Function 
The  response  functions  determined  from  the  input-output 
data  of  the  impulse  tests  were  verified  employing  Equation 
15.   Essentially,  this  equation  uses  the  impulse  load  input 
and  the  derived  response  function  for  a  particular  case. 
Computations  result  in  the  deflection  function  which  was 
previously  used  to  obtain  the  response  function.   Hence,  the 
predicted  deflection  function  should  match  the  corresponding 
measured  value.   This  checking  procedure  was  applied  to  all 
the  response  functions  as  shown  in  Appendix  A.   The  agreement 
was  observed  to  be  close  to  perfect.   Typical  plots  are 
illustrated  in  Figures  36,  37  and  38.   Table  13  shows  the 
predicted  and  measured  values  for  Series  1HI3  as  a  typical 
illustration. 

Load- Independency  of  the  Response  Function 
In  this  investigation,  the  response  functions  for  almost 
all  cases*  were  determined  at  three  different  peak  stress 
levels,  namely  14.93,  29.86  and  59.72  psi  designated  respec- 
tively, by  the  numbers  1,  2  and  3  in  the  data  coding  system. 
It  was  found  that  the  response  functions  derived  at  these 
stress  levels  for  any  test  condition,  that  is  at  the  same 
surface  course  thickness,  temperature  and  spatial  location, 
did  not  change  significantly.   Comparisons  of  representative 


*   The  response  functions  for  the  1LI1  Series  were  not  deter- 
mined because  the  deflection  functions  were  not  clearly 
recorded.   Only  the  peak  deflection  values  are  included 
in  the  analysis  of  test  results. 
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TABLE  13 
TYPICAL  PREDICTED  AND  MEASURED  DEFLECTION  FUNCTIONS 

Series :   1HI3 


Time 
sec. 

Deflection  Funct 

Lons ,  inch 

(x  10"') 

Distance 

from  Load  Center, 

inches 

3 

25 

5 

.75 

8. 

25 

Pred. 

Meas  . 

Pred. 

Meas  . 

Pred. 

Meas  . 

0 

0 

0 

0 

0 

0 

0 

0.02 

0 

0 

0 

0 

0 

0 

0.04 

0.90 

0.90 

0.80 

0.80 

0.40 

0.40 

0.06 

15.00 

15.00 

6.40 

6.40 

2.24 

2.24 

0.08 

17.10 

17.10 

9.60 

9.60 

3.44 

3.44 

0.10 

19.50 

19.50 

10.40 

10.40 

3.88 

3.88 

0.12 

20.40 

20.40 

10.72 

10.72 

4.08 

4.08 

0.14 

21.00 

21.00 

11.04 

11.04 

4.12 

4.12 

0.16 

21.15 

21.15 

11.  20 

11.20 

4.16 

4.16 

0.18 

21.30 

21.30 

11.52 

11.52 

4.24 

4.24 

0.  20 

21.60 

21.60 

11.68 

11.68 

4.28 

4.28 

0.22 

21.75 

21.75 

11.84 

11.84 

4.32 

4.32 

0.24 

12.00 

12.00 

4.80 

4.80 

3.20 

3.  20 

0.26 

4.50 

4.50 

2.80 

2.80 

1.60 

1.60 

0.  28 

2.40 

2.40 

2.08 

2.08 

1.12 

1.12 

0.30 

2.10 

2.10 

1.76 

1.76 

0.80 

0.80 

0.32 

1.80 

1.80 

1.60 

1.60 

0.76 

0.76 

0.34 

1.65 

1.65 

1.52 

1.52 

0.72 

0.  72 

0.36 

1.50 

1.50 

1.44 

1.44 

0.64 

0.64 

0.38 

1.20 

1.20 

1.36 

1.36 

0.60 

0.60 

0.40 

1.05 

1.05 

1.20 

1.20 

0.56 

'0.56 

0.42 

0.90 

0.90 

1.12 

1.12 

0.52 

0.52 

0.44 

0.90 

0.90 

1.04 

1.04 

0.48 

0.48 

0.46 

0.90 

0.90 

0.96 

0.96 

0.44 

0.44 

0.48 

0.90 

0.90 

0.88 

0.88 

0.44 

0.44 

0.50 

0.90 

0.90 

0.88 

0.88 

0.44 

0.44 
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plots  of  the  response  functions  illustrated  in  Figures  39, 
40,  and  41  support   this  observation.   Table  14  summarizes 
the  first  positive  peak,  the  time  of  its  occurrence,  the 
first  negative  peak  and  the  time  of  its  occurrence  for  each 
of  these  response  functions.   It  is  evident  from  this  table 
that  these  four  parameters  do  not  change  significantly  as 
the  peak  impulse  load  F   is  changed  at  any  spatial  location 
considered. 

Inasmuch  as  the  parameters  in  Equation  14  approximately 
represent  the  respective  response  function  curve,  the  values 
of  the  a,  3  and  y   parameters  shown  in  Tables  10,  11  and  12 
are  examined  to  further  demonstrate  the  load- independency 
of  the  response  functions.   It  is  seen  from  these  tables 
that  the  a,  0  and  y   parameters  for  each  case  do  not  vary 
appreciably  as  the  magnitude  of  load  changes.   One  exception 
with  a  distinct  variability  is  evident  in  Table  10,  and  that 
is  the  magnitude  of  a  for  the  one-inch  surface  when  x  is 
equal  to  8.25  inches  and  the  load  is  750  lbs.   It  is  possi- 
ble that  this  discrepancy  may  be  reduced  if  more  parameters 
are  included  in  the  curve  fitting  process. 

From  the  above  discussion  it  is  seen  to  be  feasible  to 
use  the  average  values  for  a,  3  and  y   in  order  to  study  the 
effect  on  these  parameters  of  the  other  test  variables. 
This  is  the  subject  of  the  next  two  sections. 
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TABLE  14 
COMPARISON  OF  RESPONSE  FUNCTIONS 
OBTAINED  FROM  DIFFERENT  LOAD  MAGNITUDES 


5.75 


Surface  Course  Thickness :   1  inch 
Test  Temperature:   100°F 


First  Positive  Peak  First  Negative  Peak 

x        F       Magnitude  Time    Magnitude  Time 

inches     lbs.   lkin>   x  10" 6  sec.  -,<in'      *10"6  sec. 

lb. -sec.  lb. -sec. 

3.25             180.0                  8.40  0.06                -3.32  0.20 

375.0                12.48  0.06                -5.25  0.22 

750.0                11.10  0.06                -4.37  0.22 


180.0 

5.19 

0.04 

-1.31 

0.22 

375.0 

6.46 

0.04 

-1.64 

0.24 

750.0 

6.01 

0.04 

-2.06 

0.22 

8.25  180.0  1.64  0.06  -0.38  0.22 
375.0  2.19  0.06  -0.76  0.22 
750.0        2.24      0.06       -0.54      0.24 


105 

Effect  of  Temperature  on  the  Response  Function 
The  parameters  a,  3  and  y   are,  respectively,  measures 
of  the  peak  of  the  response  function  G(t)  and  of  how  G(t) 
varies  with  time.   As  such,  a  can  be  regarded  as  represent- 
ing the  stiffness  characteristic  of  a  pavement  system, 
whereas  3  and  y   reflect  the  damping  characteristics  of  the 
system.   Pavement  materials  having  large  a   values  will 
provide  less  resistance  to  imposed  loads. 

Referring  to  the  curves  of  a  versus  temperature  shown 
in  Figure  42,  it  is  readily  apparent  that  temperature  plays 
an  important  role  in  the  response  function  of  pavement 
systems.   It  was  mentioned  elsewhere  that  flexible  pavement 
components,  and  asphaltic  concrete  in  particular,  are 
thermoplastic.   They  yield  increased  response  for  increases 
in  temperature,  other  test  variables  being  the  same,  and 
hence  yield  increased  a  values.   This  trend  is  indicated  by 
the  curves  of  Figure  42. 

Examination  of  the  parameters  in  Tables  10,  11  and  12 
reveals  that  changes  in  3  and  y   are  small  as  temperature 
changes.   This  is  possibly  due  to  the  fact  that  no  signifi- 
cant change  is  observed  in  the  shapes  of  the  response  func- 
tion curves.   Furthermore,  the  slight  variations  in  the 
values  of  these  parameters  are  not  critical  since  they 
appear,  respectively,  as  exponential  and  sinusoidal  functions 
in  the  response  function  model  (Equation  14) . 
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FIGURE    42       CURVES     OF     a    PARAMETER 

VERSUS     TEMPERATURE 
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Effects  of  Thickness  and  Spatial  Location 
From  the  graphical  results  shown  in  Figure  42,  it  is 
observed  that  at  x  equals  3.25  and  5.75  inches,  the  value 
of  a  decreases  with  increases  in  surface  course  thickness  at 
all  temperatures.   However,  at  x  equals  8.25  inches  a  in- 
creases with  increases  in  thickness.   This  trend  reflects 
the  slab-like  action  of  the  system  that  a  thicker  pavement 
responds  more  at  remote  locations  than  a  thinner  one. 
Further  examination  of  the  figure  shows  that  the  curves 
indicate  the  interaction  between  surface  course  thickness 
and  temperature.   For  the  one-inch  surface  the  rate  of  in- 
crease in  a  decreases  with  increases  in  temperature,  whereas 
this  rate  is  increased  for  the  two-inch  surface.   It  is 
recalled  that  similar  observations  were  made  previously 
regarding  the  effects  of  temperature  and  surface  course  thick- 
ness on  the  central  deflection  y0.   Obviously,  the  value  of  a 
decreases  with  increasing  values  of  x,  thus  depicting  the 
attenuation  of  the  magnitude  of  a  with  spatial  location. 

Referring  to  Tables  10,  11  and  12,  it  is  seen  that  the 
values  of  the  parameters  3  and  y   remain  approximately  constant 
with  changes  in  the  surface  course  and/or  spatial  location. 
This,  again,  is  probably  due  to  the  similarity  of  the  shapes 
of  the  response  function  curves . 

Static  Load  Results 
Static  load  deflections  y(t)  were  calculated  using 
Equation  23  of  "Solution  for  Step  Loading."   Computations 
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were  carried  out  employing  the  computer  program  STALOD 
listed  in  Appendix  B.   The  predicted  deflections  and  the 
corresponding  measured  values  for  all  test  series  are  sum- 
marized in  Table  15.   Typical  plots  are  illustrated  in 
Figures  43,  44  and  45.   Close  agreement  in  almost  all  cases 
is  indicated  in  this  table.   It  is  very  significant  at  this 
stage  to  recapitulate  what  has  been  done:  • 

1.  Pavement  parameters  were  determined  from  impulse 
load  tests. 

2.  The  derived  parameters  were  used  in  conjunction 
with  a  formulated  theory  to  predict  static  load 
deflections  for  various  test  conditions  of  tempera- 
ture, surface  course  thickness  and  spatial  location. 

3.  The  agreement  between  predicted  and  measured 
values  was  observed  to  be  tolerable  within  experi- 
mental error. 

Thus,  the  following  conclusions,  which  will  further  be 
supported  by  the  results  and  discussions  of  the  next  section, 
can  be  made: 

1.  The  pavement  parameters  obtained  in  this  investiga- 
tion are  indeed  descriptors  of  pavement  behavior. 

2.  Inasmuch  as  the  two  types  of  loads  are  different, 
one  being  dynamic  and  the  other  static,  these 
parameters  are  independent  of  the  type  of  input 
load. 
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TABLE  15 
PREDICTED  AND  MEASURED  STATIC  LOAD  DEFLECTIONS 


Thick- 
ness 

in . 

Temp. 

F 

0 

y(i 

t)  ,    inch 

(x    10- 

*) 

Distance 

from   Load   Center ,    inc 

hes 

°F 

lbs. 

3. 

25 

5.  75 

8. 

25 

Pred. 

Meas  . 

Pred. 

Meas  . 

Pred. 

Meas . 

1 

50 

375.0 

3.00 

4.50 

0.72 

1.20 

0.26 

0.32 

750.0 

9.95 

10.50 

3.92 

3.92 

0.  75 

0.80 

75 

187.5 

4.38 

4.16 

1.98 

2.08 

0.43 

0.48 

375.0 

9.81 

9.20 

5.48 

5.24 

1.  15 

1.12 

750.0 

20.36 

20.70 

10.90 

10.40 

2.83 

2.80 

100 

187.5 

5.05 

4.95 

2.97 

3.20 

1.03 

1.12 

375.0 

13.51 

12.15 

8.54 

8.16 

2.  77 

2.92 

750.0 

26.62 

25.60 

14.46 

14.40 

5.58 

5.60 

2 

50 

187.5 

1.02 

0.98 

0.53 

0.48 

0.16 

0.16 

375.0 

1.58 

1.80 

0.84 

0.88 

0.26 

0.  27 

750.0 

5.52 

5.55 

3.79 

3.44 

1.63 

1.44 

75 

187.5 

2.23 

2.06 

1.65 

1.60 

0.93 

0.84 

375.0 

3.95 

4.50 

2.59 

2.36 

1.  71 

1.  52 

750.0 

11.33 

11.25 

7.36 

7.04 

3.38 

3.36 

100 

187.5 

3.84 

4.35 

2.63 

2.72 

1.13 

1.20 

375.0 

9.33 

9.00 

5.19 

5.12 

2.88 

2.64 

750.0 

14.86 

15.00 

13.23 

12.00 

5.67 

5.60 
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Repeated  Load  Results 

Haversine  load  deflections  were  predicted  employing 
Equations  28  and  30,  respectively,  for  the  total  deformation 
yT(t)  and  the  permanent  deformation  yp(t).   Calculations 
were  effected  by  the  computer  program  REPLOD  of  Appendix  C. 
Table  16  shows  comparison  of  predicted  and  measured  values 
of  yT(t).   Typical  plots  are  shown  in  Figures  46,  47  and  48. 
The  agreement  is  observed  to  be  as  good  as  that  in  the  case 
of  static  loads.   The  results  indicate  that  repeated  load 
deflections,  as  well  as  static  load  deflections,  can  be 
predicted  using  the  theory  developed  in  this  investigation. 

Comparisons  of  the  calculated  deflections  in  Table  15 
for  the  static  loads  and  the  corresponding  values  in  Table 
16  for  the  haversine  loads  show  that  the  magnitudes  are 
essentially  the  same  in  both  cases.   This  is  not  surprising 
since  it  is  easily  verified  that  the  repeated  load  equation 
(Equation  28)  reduces  to  the  static  load  equation  (Equation 
23)  as  the  time  t  gets  large.   It  appears  that  the  time 
duration  of  4  seconds  used  in  haversine  load  tests  is  large 
enough  to  consider  Equation  23  as  the  limiting  case  of 
Equation  28. 

Regarding  the  measured  values  of  the  static  load  and 
repeated  load  deflections  presented,  respectively,  in  Tables 
15  and  16,  it  is  noticed  that  the  deflections  observed  from 
repeated  load  tests  are  always  less  than  those  resulting 
from  static  load  tests.   This  is  due  to  the  difference  in 
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TABLE  16 
PREDICTED  AND  MEASURED  HAVERSINE  LOAD  DEFLECTIONS 


Thick- 
ness 

in. 

Temp. 

F 

0 

yT 

ft) ,    inc 

-h    (x    10" 

") 

E 

istance 

from   Load   Center,    inc 

hes 

°F 

lbs  . 

3. 

25 

5. 

75 

8. 

25 

Pred. 

Meas  . 

Pred. 

Meas . 

Pred. 

Meas . 

1 

50 

375.0 

2.99 

3.00 

0.71 

0.80 

0.  26 

0.24 

750.0 

9.94 

9.15 

3.81 

3.36 

0.  75 

0.68 

75 

187.5 

4.37 

3.98 

1.98 

1.92 

0.43 

0.40 

375.0 

9.80 

9.08 

5.48 

5.04 

1.15 

0.96 

750.0 

20.34 

19.80 

10.88 

10.24 

2.82 

2.  76 

100 

187.5 

5.04 

4.32 

2.97 

2.52 

1.03 

0.96 

375.0 

13.50 

10.24 

8.  52 

6.08 

2.  77 

2.32 

750.0 

26.61 

24.80 

14.45 

13.12 

5.57 

4.88 

2 

50 

187.5 

1.02 

0.94 

0.53 

0.48 

0.16 

0.16 

375.0 

1.58 

1.73 

0.85 

0.84 

0.26 

0.24 

750.0 

5.52 

4.95 

3.79 

3.40 

1.63 

1.40 

75 

187.5 

2.22 

1.95 

1.65 

1.44 

0.93 

0.80 

375.0 

3.95 

3.90 

2.59 

2.  24 

1.  71 

1.44 

750.0 

11.32 

9.15 

7.35 

6.40 

3.38 

2.88 

100 

187.5 

3.84 

3.75 

2.63 

2.48 

1.13 

1.04 

375.0 

9.32 

8.25 

5.19 

5.04 

2.87 

2.48 

750.0 

14.84 

14.63 

13.22 

11.04 

5.66 

5.04 
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load  durations  which  were,  respectively,  4  seconds  and 
a  few  minutes . 

Both  the  calculated  and  measured  values  of  the  repeated 
load  deflections  did  not  vary  with  increasing  number  of  load 
applications.   From  the  theoretical  point  of  view,  the  terms 
that  would  reflect  this  change  are  exponential  (Equations  28 
and  30) .   Because  of  the  rapidly  decaying  feature  of  exponen- 
tial functions,  the  contribution  of  these  terms  to  the 
calculated  deflections  becomes  negligible  with  increasing 
number  of  load  applications.   Experimentally,  observed 
repeated  load  deflections  at  locations  other  than  directly 
under  the  load  result  from  the  first  few  cycles  and  remain 
constant  thereafter.   This  is  indicated  in  the  representa- 
tive output  traces  of  Figure  24.   The  observation  is  also 
substantiated  by  the  results  of  Drennon  and  Kenis  [62] 
obtained  from  laboratory  tests  of  a  full  scale,  three- 
layered  flexible  pavement. 

The  magnitudes  of  both  calculated  and  measured  values 
of  repeated  load  permanent  deformations  yp(t)  were  too  small 
(in  the  order  of  10" 5  inches  or  less)  for  any  meaningful 
interpretations.   Theoretically,  it  is  evident  from  Equation 
30  that  yp(t)  becomes  zero  in  the  limit  as  t  increases. 
From  the  practical  point  of  view,  field  observations  of 
flexible  pavement  behavior  indicate  that  rutting,  which  is 
mainly  a  measure  of  permanent  or  non-recoverable  deforma- 
tions, occurs  directly  under  the  wheel  loads. 
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Corner  Loading  Test  Results 
A  question  arose  whether  the  confinement  of  the  model 
pavements  in  a  box  would  have  any  effect  on  the  center  load- 
ing test  results.   Advanced  indications  of  what  might  be 
expected  were  three-fold: 

1.  Previous  laboratory  studies  dealing  with  pavements 
found  edge  effects  to  be  negligible  [58,  61]. 

2.  Coating  the  interior  sides  of  the  box  with  spar 
varnish  diluted  with  linseed  oil  would  prevent 
adhesion  of  the  layer  materials  to  the  box  sides 
and  thus  would  place  no  restriction  on  the  verti- 
cal pavement  response. 

3.  Deflections  measured  by  LVDT  no.  5  at  3  inches  from 
the  edge  during  center  loading  tests  were  small 
indicating  that  the  response  at  the  pavement  edges 
would  be  small  enough  to  minimize  any  concern  about 
the  effects  of  confinement. 

Nonetheless,  to  gain  more  confidence  in  the  over-all 
results,  a  series  of  corner  loading  tests  was  conducted  on 
the  model  pavement  with  a  two-inch  surface,  and  at  the  most 
severe  test  conditions  regarding  edge  effects.   The  two 
dynamic  loads,  namely  impulse  and  cyclic,  were  employed  at 
the  59.72  stress  level.   The  test  temperature  was  50°F. 
These  tests  also  had  as  another  objective  to  see  whether  the 
response  would  vary  much  at  different  locations  but  equi- 
distant from  the  center  of  the  load.   Thus,  any  non-uniformity 
in  the  preparation  of  the  model  pavements  would  be  detected. 
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The  test  results  of  the  corner  loadings  are  compared 
in  Table  17  with  those  of  the  center  loadings.   It  appears 
from  this  table  that  there  is  reasonable  agreement  between 
the  values  measured  at  the  inner  sides  of  the  plate  during 
corner  loading  and  the  corresponding  values  from  center 
loading.   The  corner  loading  results  indicate  that  the 
magnitudes  near  the  edge  are  slightly  larger  than  the  corre- 
sponding values  recorded  at  the  inner  sides  of  the  plate. 
This  is  probably  due  to  the  cantilever  action  that  might 
have  influenced  the  responses  near  the  pavement  edges. 
From  the  discussion  of  this  section,  it  can  be  safely  con- 
cluded that  the  edge  effects  were  negligible  and  that  the 
manner  in  which  the  model  pavements  were  prepared  is  satis- 
factory . 
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TABLE  17 
EFFECTS  OF  CONFINEMENT 


Surface  Course  Thickness :   2  inches 
Test  Temperature:   50°F 
Magnitude  of  Load:   750  lbs. 


inches 


Location* 


Measured  Deflections,  inch  O  10"1*) 

Mode  of  Loading 

Impulse  Haversine 

Center   Corner   Center   Corner 
Loading   Loading   Loading   Loading 


3.25      Edge  Side  4.88  4.73  4.95  5.03 

5.75      Inner  Side  3.36  3.28  3.40  3.44 

Edge  Side  --  3.36  --  3.52 

8.25      Inner  Side  1.36  1.28  1.40  1.44 

Edge  Side  --  1.44  --  1.52 


For  corner  loading  tests.   See  Figure  18. 
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CONCLUSIONS 

In  the  present  investigation  three-layer  flexible 
pavement  models  utilizing  one-inch  and  two-inch  surface 
course  thicknesses  were  tested.   Impulse,  static  and 
repeated  loads,  of  intensities  14.93,  29.86  and  59.72  psi, 
were  applied  at  the  center  of  the  pavements.   A  few  test 
series  were  also  conducted  at  the  corner.   Test  temperatures 
were  50°,  75°  and  100°F.   Based  on  the  results  and  within 
the  scope  of  this  study,  the  conclusions  are  enumerated  here 
It  should  be  recognized  that  these  conclusions  pertain  to 
the  materials  and  testing  procedures  used  in  this  study. 
Justifiable  extrapolation  of  the  results  should  be  made  only 
upon  further  testing. 

1.  Transfer  function  theory  is  capable  of  predicting 
static  or  repeated  load  deflections  of  flexible  pavements. 
The  favorable  agreement  between  predicted  and  measured 
values  of  the  deflections  in  this  study  validates  the 
hypothesis  that  the  parameters  in  the  response  function 
are  material  descriptors  which  are  independent  of  the  type 
of  load  input. 

2.  The  time-dependent  behavior  of  a  flexible  pavement 
can  be  represented  by  a  response  function  G(t)  which  is  a 
function  of  time  t.   It  is  possible  to  obtain  this  function 
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from  a  single  impulse  test  on  the  pavement.   The  response 
function  is  independent  of  the  magnitude  of  the  impulse  load 
justifying  the  assumption  of  linearity. 

3.  The  response  function  of  a  flexible  highway  pave- 
ment is  of  the  form 

G(t)  =  ae     sin  yt 
where  a,  $  and  y    are  descriptive  parameters  in  the  response 
function . 

4.  Parameters  in  the  response  function  are  believed 

to  be  pavement  properties  that  provide  for  better  understand- 
ing of  its  behavior  than  those  currently  used. 

5.  Temperature,  surface  course  thickness  and  spatial 
location  have  their  respective  influences  on  the  response 
function.   Increases  in  temperature  increase  the  value  of 
the  a  parameter  in  the  response  function,  while  increasing 
surface  course  thickness  or  the  distance  from  the  load 
center  decreases  the  value  of  a.   The  3  and  y   parameters 

do  not  seem  to  be  affected  appreciably  by  the  above  factors . 

6.  The  profile  of  peak  deflections  of  a  flexible 
highway  pavement  can  be  described  by  the  equation 

yU)  ■  y0e 

where  y(x)  is  the  deflection  at  a  distance  x  from  the  load 
center,  y0  is  the  central  deflection  and  D  is  a  constant 
reflecting  the  attenuation  of  the  deflection  profile  with  x. 

7.  The  central  deflection  y0  increases  with  increases 
in  temperature,  and  decreases  with  increasing  surface  course 
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thickness.   Increases  in  temperature  (50°F-  100°F)  or 
surface  course  thickness  (1  inch  and  2  inches)  decrease  the 
value  of  the  parameter  D. 

8.   Instrumentation  and  experimental  procedures  con- 
ducted were  satisfactory  for  the  test  conditions  of  this 
investigation.   The  edge  effects  due  to  the  confinement  of 
the  model  pavements  in  the  box  are  negligible  with  regard 
to  the  results  of  this  study. 
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SUGGESTIONS  FOR  FURTHER  RESEARCH 

The  results  of  this  project  have  been  based  on  labora- 
tory testing  of  flexible  model  pavements.   It  is  realized 
that  this  is  a  first  attempt  in  applying  transfer  function 
theory  to  highway  pavements  and  more  work  needs  to  be  done. 
The  following  research  proposals  are  recommended  as  exten- 
sions of  this  investigation: 

1.   In  this  study,  the  response  functions  were  approxi- 
mated by  the  mathematical  model 

G(t)  =  ae"et  sin  Yt  Model  (1) 

The  following  models  could  be  used: 

G(t)  =  ae"6t  sin  yt    -  6  sin  et       Model  (2) 
G(t)  =  ae"gt  sin  yt  +  6  cos  et       Model  (3) 

where  G(t)  is  the  response  function,  t  is  time  and  a,  0,  y» 
6  and  e  are  parameters. 

Preliminary  analysis  of  some  of  the  data  indicated  that 
using  Models  2  and  3  improved  the  squared  correlation  coef- 
ficient (R2),  as  shown  in  Table  18.   The  theory  of  this 
investigation  can  be  extended  after  replacing  Model  1  by 
either  Model  2  or  3. 
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TABLE  18 
COMPARISON  OF  R2  FOR  DIFFERENT  MODELS 


e  X 

Series 


R: 


inches    Model  (1)    Model  (2)    Model  (3) 


1HI1  3.25  84.62  86.98  87.73 

5.75  88.75  91.38  91.85 

8.25  88.57  91.14  92.39 

1HI2        3.25  83.23  90.26  86.89 

5.75  83.04  89.52  90.52 

8.25  78.24  87.10  86.94 


1HI3 

3.25 

84.30 

88.38 

88.74 

5.75 

81.78 

88.11 

87.53 

8.25 

84.92 

91.89 

91.81 

2HI1 

3.25 

83.11 

85.24 

89.42 

5.75 

87.90 

92.25 

91.48 

8.25 

84.81 

91.40 

88.  70 

2HI2 

3.25 

86.90 

88.14 

90.22 

5.75 

89.93 

95.07 

93.36 

8.25 

92.34 

93.75 

96.77 

2HI3 

3.25 

80.87 

83.23 

89.71 

5.75 

89.15 

90.03 

93.85 

8.25 

93.19 

94.78 

96.38 
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2.  The  results  of  this  investigation  indicate  feasibil 
ity  of  conducting  field  tests  on  highway  pavements.   The 
signatures  from  a  truck  travelling  over  a  bump  can  be  used 
to  derive  the  response  functions.   Application  of  transfer 
function  theory  to  evaluate  airport  pavements  has  just  been 
completed  [64] . 

The  test  program  should  cover  a  variety  of  construc- 
tion types  based  on  a  large  statistical  sample  to  permit 
quantitative  evaluation  of  the  effects  of  the  test  factors 
involved.   In  the  present  study,  the  effects  of  variables 
were  evaluated  qualitatively. 

3.  The  mathematical  function  obtained  in  this  study  to 
describe  pavement  deflection  profiles  might  be  correlated 
with  failure  in  flexible  pavements.   Available  data  or 
deflections  measured  from  planned  field  tests  may  be  used 
for  this  purpose.   By  relating  the  slope  of  the  deflected 
profile  to  the  nature  and  extent  of  pavement  failure,  a 
rational  method  of  evaluation  could  emerge. 

4.  The  pavements  in  this  investigation  were  tested  in 
their  entirety  and,  as  such,  isolation  of  the  material 
properties  for  each  layer  was  not  possible.   A  research 
program  that  would  allow  determination  of  the  component 
properties  would  include  testing  a  single  subgrade  layer 
and  successively  building  up  to  the  testing  of  a  complete 
pavement . 
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Limited  experimental  work  with  a  3-inch  diameter  plate 
in  this  study  and  the  results  of  other  investigators  [70] 
indicated  that  the  loaded  area  may  have  an  effect  on  pave- 
ment properties.   Variation  of  this  factor  may  be  incorporated 
in  the  laboratory  testing  program. 

5.  In  this  investigation,  the  measured  output  parameter 
was  deflection,  based  on  the  premise  that  surface  deflections 
are  indicative  of  a  flexible  pavement's  ability  to  withstand 
stresses  produced  within  the  pavement  by  imposed  loads  [71, 
72].   However,  this  places  no  restriction  on  measuring  other 
parameters,  such  as  stress  or  strain.   Inventory  tests  con- 
ducted in  this  study  using  strain  gages  showed  that  instru- 
ments with  better  resolution  would  be  needed  to  record 
strain  magnitudes  of  the  order  of  20  microns. 

6.  In  the  present  study,  pavement  response  was  analyzed 
without  resort  to  elastic  or  viscoelastic  theory.   However, 
individuals  interested  in  extending  the  application  of  such 
theory  can  use  the  experimental  technique  of  this  investiga- 
tion.  By  impulse  testing  material  specimens,  more  reliable 
properties  can  be  determined  and  used  in  conjunction  with 
these  theories. 

7.  The  parameters  determined  in  this  investigation  may 
be  correlated  and  compared  with  those  obtained  by  other 
researchers . 
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Program 

PROGRAM  MAIN  (INPUT  .OUTPUT  , PUNCH  , TAPE 5=  INPUT  ,TAPIi6  =  OUTPUT  , 
1 TAPE 7=PUNCH, PLOT) 

C 

r  ****************************************************** 

c 

c 

C  COMPUTER  PROGRAM  TO  CALCULATE  RESPONSE 

C        FUNCTIONS  AND  PREDICT  DEFLECTIONS  AS  A  CHECK  ON  THE 

C         RESPONSE  FUNCTIONS 

C 

C 

r  ****************************************************** 

C 

C  DESCRIPTION  OF  VARIABLES  USED 

C  I         IDENTIFICATION  FOR  DEFLECTION  RESPONSE  FUNCTIONS 

C  J         =DATA  SET  NUMBER. 

C  N         =NUMBER  OF  SEQUENCE  VALUES  IN  THE  DEFLECTION 

C  RESPONSE  FUNCTIONS. 

C  IPLOT     =COMMAND  FOR  PLOT  OF  VARIABLES. 

C  DT        = INCREMENT  OF  TIME  BETWEEN  EACH  SEQUENCE  VALUE 

C         IN  THE  DEFLECTION  RESPONSE  FUNCTIONS. 

C  DEFL(I ,K)=ARRAY  OF  SEQUENCE  VALUES  FOR  I -TH  DEFLECTION 

C  RESPONSE  FUNCTION. 

C  GT(I,K)   =ARRAY  OF  SEQUENCE  VALUES  FOR  I  - TH  TIME 

C  DEPENDENT  TRANSFER  FUNCTION. 

C  GF(I,K)   =DUMMY  VARIABLE  FOR  TEMPORARY  STORAGE  OF 

C  SEQUENCE  VALUES. 

C  IPS  (I)    =NUMBER  OF  SEQUENCE  VALUES  FOR  PHASE  SHIFT  OF 

I-TH  DEFLECTION  RESPONSE  FUNCTION. 

C  GXX       =PEAK  LOAD 


DIMENSION  Vl(40) , V2 (40) -V3 (40) ,V4(40) ,V5(40) ,V6(40) 

DIMENSION  XINP(40)  ,XOUT(40) ,XG(40)  ,Y(40)  ,Z(40)  ,W(6) 

DIMENSION  DEFL(6,4  0) ,GF(6,40)  ,GT(6,40)  ,DE(6,40) 

COMMON/AXIS/  X(50) 

DIMENSION  MONT(IO) 

DIMENSION  IPS(12) ,M912) 

DIMENSION  PEAK(IO) ,CONT(10) ,SPK(10) 

DIMENSION  NAME (10) 

INTEGER  CONT, BEGIN, END 

REAL  MOD 

READ (5, 100)  BEGIN,END,N,DT,IPLOT,INEG,IMOD,ISMOTH 

N1=N 

DO  1000  J=BEGIN,END 

WRITE(6,208)  J 

READ(5,106)  (NAME (L) ,L=1 ,9) 

106  FORMAT (4X,9A8) 
WRITE(6,107)  (NAME(L) ,L=1 ,9) 

107  F0RMAT(1H1,////1X,9A8) 

re  ad ( 5 , 1 0 1 )  1  DATA , I  FLAG , JDATA , JFLAG , KDATA , KFLAG , LDATA , 
1LFLAG.MDATA, 
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2MFLAG,NDATA,NFLAG, (M (I ) , I  =  1 , 6)  , (1  PS (I ) , I  =  1 , 6  J  , GXX 
GO  TO  (29,700) ,IFLAG 

29  DO  30  1=1,6 

30  READ(5,102)  (DEFL ( I ,K) , K=l ,N) 
WRITE(6,200) 

DO  261  K=1,N 

261  X(K)=FLOAT(K)*DT-DT 
DO  262  K=1,N 

262  WRITE(6,201)  X (K) , (DEFL (I ,K) , 1=1 , 4) 

500  DO  501  1=1,6 
DO  501  K=1,N 

501  GF(I,K)=0.0 

IF(IFLAG.EQ.l)  GF (2 ,NPHASE) = SPIKE 

DO  504  1=1,5 

IF(M(I) .EQ.0)GO  TO  504 

DO  502  K=1,N 

XINP(K)=DEFL(1,K) 

502  XOUT(K)=DEFL(I+l,K) 

CALL  CONVLI (XINP ,XOUT ,XG ,DT ,N , GXX J 

DO  503  K=1,N 

GT(I+1,K)=XG(K) 

503  GF(I+1,K)=XG(K) 

504  CONTINUE 

WRITE(6,107)  (NAME(L)  ,L  =  1  ,9) 
WRITE(6,203) 
1572  DO  505  K=1,N 

X(K)=FLOAT(K)*DT-DT 

WRITE (6, 2 02)  X  (K)  ,  (GF (I ,K)  , 1  =  2  ,4) 

505  CONTINUE 
C1120  PUNCH 

IF(JFLAG.EQ.  0)  GO  TO  525 

DO  530  K=1,N 

V2(K)=0.0 

V3(K)=0.0 

V4(K)=0.0 

V5(K)=0.0 

V6(K)=0.0 
530  CONTINUE 
C 

C      DETERMINE  SCALING  PARAMETERS  FOR  RESPONSE  FUNCTIONS 

C 

DO  506  K=1,N 

V2(K)=GF(2,K) 

V3(K)=GF(3,K) 

V4(K)=GF(4,K) 

V5(K)=GF(5,K) 

V6(K)=GF(6,K) 

506  CONTINUE 

CALL  SCALES (10 . 0 ,V2 ,N , 1 ,V3 ,N , 1 ,V4 ,N , 1 ,V5 ,N , 1 , V6 ,N , 1 ) 

CALL  SCALE(x,10.0,N,l) 
C 

C      PLOT  AXES  AND  ANNOTATION 

C 
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CALL  AXI S ( 0. 0, 0. 0.12HTIME, SECONDS, -1 2, 1 0 . 0 , 0 . 0 ,X (N  +  l )  , 
lX(N+2) ,0) 

CALL  AXIS(0.0,0.0,17HRESPONSE  FUNCTION , 1 7 , 1 0 . 0 , 9 0 . 0 , 
1V2(N+1) ,V2(N+2) ,-1) 

CALL  SYMBOL(3.5,10.0,0.2,18HRESPONSE  FUNCTIONS , 0. 0 , 18 j 

CALL  SYMBOL (8. 0,8. 2 0,0. 0  8,1 ,0.0,-1) 

CALL  SYMBOL(8.1,8.14,0.08,11H-RESPONSE  1,0.0,11) 

CALL  SYMBOL (8. 0,8.1 ,0.08,2,0.0,-1) 

CALL  SYMBOL(8.1,8.04,0.08,11H-RESPONSE  2,0.0,11) 

CALL  SYMBOL (8 . 0,8.0,0.08,3,0.0,-1) 

CALL  SYMBOL(8.1,7.94,0.08,11H-RESPONSE  3,0.0,11) 

CALL  SYMBOL (8.0,7.9,0.08,4,0.0,-1) 

CALL  SYMBOL(8.1,7.84,0.08,11H-RESPONSE  4,0.0,11) 

CALL  SYMBOL (8.0,7.8,0.08,5,0.0,-1) 

CALL  SYMBOL(8.1,7.74,0.08,11H-RESPONSE  5,0.0,11) 

CALL  SYMBOL (8.0,7.7,0.08,6,0.0,-1) 

CALL  SYMBOL(8.1,7.64,0.08,11H-RESPONSE  6,0.0,11) 
C 

C      PLOT  RESPONSE  FUNCTIONS 

C 

DO  517  1=1,5 

IF(M(I) .EQ.0)  GO  TO  517 

DO  509  K=1,N 
509  Z(K)=GF(I,K) 

Z(N+1)=V2(N+1) 

Z(N+2)=V2(N+2) 

CALL  LINE(X,Z,N,1,2,I) 
517  CONTINUE 

CALL  PLOT(18. 0,0. 0,-3) 
525  CONTINUE 
C 

700  IF(MDATA.EQ.O)  GO  TO  702 

701  READ(5,102)  (DEFL  (1 ,K)  ,K  =  1 ,N) 

702  IF(MFLAG.EQ.O)  GO  TO  704 
DO  703  1=2,4 

703  READ(5,102)  (GT (I  ,K)  ,K=1 ,N) 

C      CONVERT  STORED  AND/OR  INPUT  DATA  TO  SUBROUTINE  FORM 

704  DO  706  K=1,N 
GF(1,K)=DEFL(1,K) 

706  XINP(K)=DEFL(1,K) 
XINP(3)=10.*GXX 

707  DO  710  1=2,4 
DO  708  K=1,N 

708  XG(K)=GT(I,K) 

CALL  CONVLE(XINP,XOUT,XG,DT,N) 
DO  709  K=1,N 

709  GF(I ,K)=XOUT(K) 

710  CONTINUE 

C      WRITE  DEFLECTION  RESPONSES 
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WRITE(6,107)  (NAME(L) ,L=1,9) 

WRITE(6,207) 

DO  711  K=1,N 

711  X(K)=FLOAT(K)*DT-DT 
DO  712  K=1,N 

712  WRITE(6,201)  X(K) , (GF( I ,K) , 1=1 , 4) 
1001  N=N1 
1000  CONTINUE 

IF(IPLOT.NE.O)  CALL  PLOT (0 , 0 , 999) 

C 

C 

C 

C      FORMAT  STATEMENTS 

C 

r  ****************************************************** 


****************************************************** 


c 

100  FORMAT(313,F10. 5,411) 

101  FORMAT(18I1,6I2,F10.5) 

102  FORMAT (8F10. 5) 

103  FORMAT (12) 

108  FORMAT (8F1 0.8) 

200  FORMAT (///25X, 'RECORD  OF  RAW  DATA'  , ///8X  ,' TIME ', 2X , 

1' LOAD, LBS. »  ,3X,'D    EFLECTION    S  ,  IN . '  , // 7X  ,  ' (SEC) ' 
2,4X,'X  =  0'  ,6X,'X=3.25'  ,5X,'X=5.75'  , 5X , ' X  =  8 . 25  * /// ) 

201  FORMAT (( IX, Fll. 2 ,F8 . 1 , 3 (Fll . 6) ) 

202  F0RMAT(1X,F12.3,3(F14.8)) 

203  FORMAT (///,8X, 'RESPONSE    FUNCTIONS    FROM    IMPLICIT    CONVO ' 
l'LUTION'//8X, 'TIME* ,7X, 'R   ESPONSE      FUNCTI* 
2'0    N    S7/7X, '  (SEC)  '  ,7X, 'X=3.25*  ,8X, •X=5.75l  ,8X, »X=8.25' 
3///) 

204  F0RMAT(1X,F12.6,12X,5F12.6) 

206  F0RMAT(1X,2F12.6) 

207  FORMAT (///,9X, 'DEFLECTIONS  FROM  RESPONSE  FUNCTIONS  AN' 
I'D  LOAD'  ,///8X, 'TIME'  ,2X, 'LOAD, LBS.  '  ,3X, 'D  E  F  L  E  C  ' 
2'T  I  O  N  S,IN.  7/7X, ' (SEC)  '  ,4X,'X=0'  ,6X,  *X=3.25'  , 5X , 
3'X=5.75'  5X  'X=8.25'///) 

208  FORMAT (1H1,1X///////////////////////// IX ,' 

*  DATA  SET  NUMBER'  ,2X, 12) 

210  FORMAT(21(8E10.3/) ,7E10.3) 
214  FORMAT(2F10.8) 
C 

C 

STOP 
END 

C 

C     SUBROUTINE  CONVLI 
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C  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = =  =  =  =  =  =  =  =  = 

C  PURPOSE 

C         TO  COMPUTE  THE  VECTOR  OF  VALUES  OF 

C  XG  FROM  CONVOLUTION  EQUATION  XINP*XG=XOUT 

C  FOR  GIVEN  XINP(I)  AND  XOUT(I) 

C        DISCRIPTION  OF  PARAMETERS 

C        XINP- INPUT  VECTOR  OF  VALUES 

C         XOUT- OUTPUT  VECTOR  OF  VALUES  AFTER  CONVOLUTION 

C         XG-VECTOR  OF  VALUES  OF  CONVOLUTION  FUNCTION 

C         DT- INCREMENT  OF  THE  SUMMATION  POINTS 

C         N-NUMBER  OF  SUBINTERVAL  PLUS  1. 

C      SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  NEEDED 

C         NONE 

C      METHOD 

C         1.   SET  THE  VALUE  OF  'TEST'.  THEN  TEST  THE  VALUES 

C       OF  XINP(I)  TO 

C  SEE  WHETHER  XINP (I)  IS  GREATER  THAN  TEST  OR  NOT 

C         2.   IF  MAGNITUDE  OF  XINP(K)  IS  GREATER  THAN  'TEST' 

C      THEN  MAKE  UP  K  MORE  VALUES  OF  XOUT (I J  AND  K  MORE 

C         SUMMATION  EQUATIONS 

C         5.   CONSIDER  THE  SUMMATION  EQUATION   WHICH  IS 

C      TRANSFERED  FROM  INTRGRATION  EQUATION  OF  CONVOLUTION. 

C         4.   BY  FORWARD  SUBSTITUTION  METHOD  TO 

SOLVE  FOR  XG(I)  ,1  =  1,. . .  ,N. 
C      

SUBROUTINE  CONVLI (XINP , XOUT ,XG ,DT ,N ,GXXJ 
C      

DIMENSION  XINP(l) ,XOUT(l) ,XG(1) 

DIMENSION  Z(200) ,Y(200)  ,W(200) 

A=0. 

TEST=. 000001 

DO  150  1=1, N 

150  XG(I)=0.0 

DO  151  K=1,N 

IF(ABS(XINP(K)) .LT.TEST)  GO  TO  151 

L  =  K 

GO  TO  152 

151  CONTINUE 

152  IN=L-1 
IM=N+L-1 

IF(IN.EQ.O)  GO  TO  156 
XDOT=XOUT (N) /FLOAT (IN) 
DO  153  K=1,IN 

153  XOUT (N+K) =XOUT (N) -XDOT*FLOAT (K) 
156  CONTINUE 

DO  155  I=L,IM 

11=1-1 

DO  154  K=L,II 
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154  A=A+XG(K-L+1)*XINP(I-K+L) 
XINP(L)=10.*GXX 
XG(I-L+1)=(X0UT(I)-A*DT)/(XINP(L)*DT) 

155  A=0.0 
RETURN 
END 

C  llllllllllllllllllllllllllllll lllllllllll lllllllll llll 

C  SUBROUTINE  CONVLE 

C  

C  PURPOSE 

C  TO  COMPUTE  THE  VECTOR  OF  VALUES  OF  XOUT(I)  FOR  GIVEN 

C  VECTORS. . .XINP(I)  AND  XG(I)  BY  USING  CONVOLUTION 

C  (I.E.  CONVOLUTE  XINP (I)  AND  XG(I)  BY  SUMMATION  METHOD) 

C  DISCRIPTION  OF  PARAMETERS 

C  XINP- INPUT  VECTOR  OF  VALUES 

C  XOUT-OUTPUT  VECTOR  OF  VALUES  AFTER  CONVOLUTION 

C  XG-VECTOR  OF  VALUES  OF  CONVOLUTION  FUNCTION 

C  DT- INCREMENT  OF  THE  SUMMATION  POINTS 

C  N-NUMBER  OF  SUBINTERVAL  PLUS  1. 

C 

C  SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  NEEDED 

C  NONE 

C  METHOD 

C  TRANSFER  THE  INTEGRATION  FORMULA  OF  CONVOLUTION  INTO 

C  FINITE  SUMMATION  EQATIONS  WITH  SUBINTERVALS  OF  LENGTH  DT 

C  AND  NUMBER  N-l.  COMPUTE  XOUT(I)  BY  GIVEN  VALUES  OF  XINP 

C  AND  XG  DIRECTLY 


C 


SUBROUTINE  CONVLE (XINP ,XOUT ,XG ,DT ,N) 

DIMENSION  XINP(l)  ,XOUT(l) ,XG(1) 

DO  100  1  =  1, N 

XOUT(I)=0. 

DO  150  J-1,1 

XOUT (I ) =XOUT ( I ) +XG (J) *XINP (I  - J  +  l) *DT 
150  CONTINUE 
100  CONTINUE 

RETURN 

END 
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Reduced  Data 
The  reduction  of  the  impulse  test  data  was  implemented 
using  the  program  of  this  Appendix.   Each  test  series 
consists  of  three  tables.   These  are:   input  and  outputs 
(raw  data)  ,  response  functions  and  predicted  deflection 
functions,  in  that  order. 
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i 


• 


11.  T? 


RECORD    OF    PAW    nATA 


TIME       LOAD»l  PS.  P    F    F    I      ErTTOM    StIN. 

(SEC)  X  =  0  X=3.?^  X  =  5. 75  X  =  «,^ 


0.00  0.0  0.000000  0.000000  0.000000 

.0?  o.o  0.000000  o.ooonoo  0.000000 

.04  100,0  .000060  .000013  ,000004 

,0ft  ?00,0  .000150  .000016  .ooooo* 

,Ofl  375,0  ,000173  .000035  ,00001? 

,10  375,0  ,00019c,  ,000045  .000016 

.1?  375,0  .000218  .000061  .000020 

,14  375,o  ,000??5  ,000067  .000024 

.16  375.0  .000240  .000070  .00002« 

,lfl  375,o  ,000?44  .000074  .00003? 

20  375.0  .00024P  .000077  .00003? 

.??  0.0  .000165  .O0004P  ,0000?0 

.24  0.0  .00009R  ,00003?  .00001? 

.2*  0.0  .00006P  .000013  .00000« 

.?fl  0.0  .000049  .000010  .000004 

30  O.o  .000045  .00000*  0.000000 

3?  0.0  .000038  .000005  0.000000 

.34  0.0  .000030  .000003  0.000000 

.36  0.0  .000023  .000003  0.000000 

•3fl  0.0  .000015  .000003  0.00000^ 

.40  O.n  .000008  .00000"*  0.000000 

.4?  0.0  .000008  .00000?  0.000000 

.44  O.o  .000008  .000003  0.000000 

.46  O.o  .000008  .00000?  0.000000 

.4P  0.0  ,000008  .000003  0.000000 

.50  0.0  .000008  .00000^  0.000000 
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1LT? 

HESPONSF    FUNCTIONS    FROM    T^Pl  TCTT    CONVOLUTION 
TIMF  PFSPONSF       CU    MOTIONS 

(SEC)  X=3.?5  X=*.7^  X=8.25 


.000 

.00000080 

.00000017 

•OOOOOOO* 

.0?O 

.00000196 

,00000020 

.0000001 0 

.040 

.00000212 

.00000044 

.0000001^ 

•  060 

.0^000221 

.00000054 

.0000001P 

.OHO 

.0O000?30 

.00000070 

.Q00000?3 

.100 

.00000217 

.0000007? 

.000000?6 

.1?0 

.00000215 

,fl0on0069 

.000000?9 

.140 

.00000198 

.00000067 

.00000011 

.uo 

.00000183 

.O0OOO0A4 

,000000?8 

•  l«o 

.00000061 

.00000021 

.00000010 

•  ?00 

-.00000020 

-.00000003 

-.0000000? 

•  ??0 

-.00000041 

-.P00000?4 

-.OOOOOOOf, 

•  ?40 

-.00000041 

-.00000021 

-.00000009 

.260 

-. 0^000019 

-.00000018 

-.0000001 1 

.?80 

-.00000004 

-.00000010 

-.OOOOOOOR 

.300 

.00000008 

-.00000004 

-.00000004 

.320 

.00000018 

.00000003 

-.00000000 

•  340 

.00000025 

.0000001 0 

.00000003 

.360 

.00000019 

.0000001 1 

.00000004 

.380 

.00000015 

.00000010 

.00000003 

.400 

.00000009 

.00000006 

.0000000? 

.420 

.00000004 

.00000003 

.oooooooi 

.440 

.00000001 

,00000001 

-.00000000 

.460 

.00000001 

-.00000000 

-.oooooooi 

.4R0 

-.00000004 

-.00000002 

-.oooooooi 

.500 

-.00000007 

-.00000004 

-.00000001 
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1LT2 

DFFLFCTIONS  FRCM  RFSPONSF  FUNCTIONS  ANn  LOAP 

TIMF      LOAD»lPS.         DFFLErTION    S«IN. 
(SEC)  *  =  0  X=3.2*  Xr5.75  Xsfl.?^ 


•  00 

0.0 

0.000000 

0.000000 

0.000000 

.0? 

0,0 

0.000000 

0.000000 

0.000000 

.04 

100.0 

.000060 

.0000) 3 

.000004 

.06 

200.0 

.000150 

,00001ft 

.OOOOOP 

.OR 

375.0 

.000173 

.00003* 

.000012 

.10 

375,0 

.OOOlQc, 

,onC04P 

.00001ft 

.12 

375.0 

.00021R 

.000061 

.000070 

.14 

375.0 

.0002?^ 

.000067 

.000074 

.1ft 

375.0 

.000240 

.000070 

.0000?P 

.18 

375.0 

.000244 

.000074 

.000032 

.20 

375,n 

.00024P 

.000077 

.000032 

.22 

0,0 

.OOOlftS 

.000048 

.000070 

.24 

0.0 

.00009P 

.000032 

.000012 

.2ft 

0,0 

.000068 

.000013 

.OOOOOP 

.28 

0,0 

.000049 

.000010 

.000004 

.30 

0.0 

.000045 

.00000ft 

-.000000 

.32 

0.0 

.000038 

.000005 

.000000 

.34 

0.0 

.000030 

.000003 

0.000000 

.3ft 

O.n 

.000023 

.00000"^ 

-.000000 

.3R 

0.0 

.000015 

.000003 

-.000000 

.40 

O.o 

.000008 

,000003 

-.000000 

.4? 

0.0 

.000008 

.000003 

-.000000 

.44 

O.o 

.000008 

.000003 

.000000 

.46 

0.0 

.000008 

.000003 

0.000000 

.48 

0.0 

.000008 

.00000^ 

.oooooo 

.50 

0,0 

.000008 

.000003 

.000000 
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1LI3 


RECORD  OF  PAW  DATA 


TIME      LOADtLPS.         DEFLECTION    SfIN, 
(SEC)  X  =  0  X«3.25  X.5.75  X«8.25 


0.00 

0.0 

0.000000 

0,000000 

0.000000 

.0? 

0.0 

0.000000 

0.000000 

0,000000 

.04 

400,0 

.00007S 

.000032 

.000010 

.06 

720.0 

.000525 

.0001P2 

.000022 

•  OP 

750.0 

.000675 

.000256 

.000051 

.10 

750.0 

.000705 

.000296 

.000059 

•  1? 

750.0 

.000720 

.000300 

.OOOOM 

.14 

750.0 

.000735 

.000300 

.000062 

.16 

750.0 

.000750 

.000304 

.000064 

•  IP 

750.0 

.000758 

.000304 

.000066 

.26 

750.0 

.000765 

.000304 

.000067 

.22 

0.0 

.000450 

.000152 

•000048 

.24 

0.0 

.000165 

.000032 

•000019 

.26 

0.0 

.000105 

.000024 

.000010 

.28 

0.0 

.000090 

.000016 

.000006 

.30 

0.0 

.000075 

.000012 

.000003 

.32 

0.0 

.000060 

.000008 

,000002 

•  34 

0.0 

.000045 

.000008 

.000002 

.36 

0.0 

.000030 

.00000P 

.000002 

.38 

0.0 

.000023 

.000008 

.000002 

•  40 

0,0 

.000015 

.000008 

.00000? 

.42 

0,0 

.000015 

.000008 

.00000? 

.44 

0.0 

.000015 

.000008 

.0000ft? 

.46 

0.0 

.000015 

.000008 

.00000? 

.4R 

o.o 

.000015 

.000008 

.00000? 

.50 

O.o 

.000015 

.OOOOOP 

.00000? 
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1  LI? 


RESPONSE    FUNCTIONS    FROM    lvP|  TCIT    CONVOLUTION 
TIME  RFSPOMSF      FUNCTIONS 

(SEC)  X=3.25  X  =  =;.7«5  X  =  fl.?«^ 


.000 

.00000050 

.00000021 

.00000007 

.020 

.00000345 

.00000126 

.00000014 

•  040 

.00000412 

.00000156 

.0000003? 

.060 

.00000391 

.00000168 

.00000034 

.080 

.00000362 

.00000154 

.00000032 

.100 

.00000335 

.0000013S 

.OOOOOOOO 

.170 

.00000312 

.00000127 

.O00000?fl 

.140 

•0O000?86 

.000001 14 

.00000026 

.1*0 

.00000262 

.00000103 

,000000?4 

•  180 

•O0000031 

-.00000007 

.0000001 0 

•  200 

-.00000129 

-.00000074 

-.00000009 

•220 

-•000001  15 

-.00000057 

-.0000001 } 

.240 

-.00000075 

-.00000039 

-.OOOOOOOQ 

•  260 

-.00000041 

-.0*000023 

-.00000007 

.280 

-.00000013 

-.00000009 

-.00000004 

.300 

.00000009 

.00000004 

-.oooooooi 

.320 

.00000027 

.00000015 

.0000000? 

.340 

.00000046 

.00000024 

.00000004 

.360 

.00000039 

.O00000?l 

.00000005 

.380 

.00000022 

.00000012 

.00000003 

.400 

.00000009 

.oooooooo 

.00000002 

.420 

.oooooooo 

.oooooooo 

.oooooooi 

.440 

-.00000004 

-.00000002 

•oooooooo 

•  460 

-.00000005 

-.00000003 

-.OOOOOOOO 

•  480 

-.0000000* 

-.oooooooo 

-.oooooooi 

.500 

-.ooooooio 

-.00000005 

-.oooooooi 
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ilt? 


DEFLECTIONS  FROM  RESPONSF  FUNCTIONS  AND  LOAD 

TIME       LOAD»L«S.  D    F    F    I     ErTTON    S.IN. 

(SEC)  X=0  X=3.?S  Xr5.75  XsP.PS 


LOO 
.0? 
.04 
.06 
.08 
.10 
.1? 
.14 
.16 
.18 
.20 
.22 
.24 
.26 
.28 
.30 
.32 
.34 
.36 
.38 
.40 
.4? 
.44 
.46 
.48 
•  SO 


0 

0 

40C 

720 

750 

750 

750 

750 

750 

750 

750 

0 

0 

0 

0 

0 

0 

0 
0 
0 
0 
0 
0 
0 
0 
0 


0.000000 
0.000000 
.00007s 
.0005?5 
•00067s 
.00070S 
.0007?0 
.000735 
.000750 
.000758 
.000765 
,000450 
.000165 
.000105 
.000090 
.00007s 
.000060 
.000045 
.000030 
.000023 
.000015 
.000015 
.000015 
.00001S 
.000015 
.000015 


0.000000 
0.000000 

.000032 
.00019? 
.0002^6 
.000296 
.000300 
.000300 
.000304 
.000304 
.000304 
.00015? 
.00003? 
.000024 
.000016 
.000012 
.000008 
.000008 
.000008 
.000008 
.000008 
.000008 
.000008 
.000008 
.00000* 
.000008 


0.000000 
0.000000 
.000010 
.0000?? 
.000051 
.000OS9 
.0000* 1 
.00006? 
.000064 
.0000*6 
.000067 
.000048 

.000019 

.ooooio 

.00000* 
.000003 
.00000? 
.00000? 

.ooooo? 

.00000? 
.00000? 
.00000? 
.00000? 
.00000? 

.ooooo? 

.00000? 
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1MI1 


RECORD  OF  PAW  DATA 


TIME      LOAD»LPS.         DFFLEcTION    S.IN. 
(SEC)  X«0  Xs3,?5  Xe5.75  X«8,?5 


.00 

0.0 

0.000000 

0.000000 

0.000000 

.02 

0,0 

0.000000 

0.000000 

0.000000 

.04 

100.0 

.000045 

.000040 

,000008 

.06 

150.0 

.00025Q 

.000104 

,000019 

.08 

187.  5 

.000270 

.000120 

,000022 

.10 

187.5 

.000285 

.000128 

,000029 

.1? 

187.5 

. 000293 

.000136 

.000032 

.14 

187.? 

.000300 

.000140 

.000034 

.16 

187.5 

.000308 

.000144 

.000035 

.18 

187.5 

.000315 

.000148 

,000037 

•  20 

187.  « 

.000323 

.000152 

,000038 

.22 

0.0 

.000195 

,000080 

,000029 

.24 

O.p 

.000135 

.000056 

,000019 

.26 

0.0 

.000008 

.000048 

,000016 

.28 

0.0 

.000045 

.000040 

,000013 

.30 

0.0 

.000038 

.000032 

,000010 

.32 

0.0 

.000030 

.000024 

,000006 

.34 

0.0 

.000023 

.000016 

,000003 

.36 

0.0 

.000015 

.000008 

,000003 

.38 

0.0 

.000008 

.000008 

.000003 

.40 

0.0 

.000008 

.000008 

,000003 

.42 

O.n 

,000008 

.000008 

,000003 

•  44 

0.0 

.000008 

.000008 

,000003 

.4ft 

0.0 

.000008 

.000008 

.000003 

.48 

0,0 

.000008 

.000008 

.00000"* 

.50 

0.0 

.000008 

.000008 

.000003 
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1  MI  1 


RESPONSE  FUNCTIONS  FPOM  IMP|  TCTT  CONVOLUTION 
TIME        RESPONSF   FUNCTIONS 
(SEC)         X«3.25         X=*.7*  X*8.2* 


1.000  .00000120  .00000107  .000000P1 

.020  ,00000681  .00000269  .00000049 

.040  .00000654  .000002*8  .000000*3 

.060  .00000628  ,000002«1  .000000*6 

.080  .00000586  .00000274  ,000000*8 

.100  .00000545  ,00000257  .000000*6 

.120  .00000511  .00000242  .000000*2 

.140  .00000478  .00000228  .000000*1 

.160  .00000451  .0000021*  .000000*8 

•190  .00000076  .00000012  .00000030 

•200  -.00000031  -.00000030  .0000000* 

.220  -.00000064  -.00000020  .00000001 

•240  -.00000136  -.00000011  -.00000000 

•260  -.00000084  -,00000004  -.00000002 

•280  -.00000042  .OOOnOOOl  -.00000006 

.300  -.00000004  .00000004  -.00000007 

.320  .00000023  .00000005  -.00000000 

•340  .00000048  .00000026  .0000000ft 

.360  .00000051  .00000025  .00000008 

.380  .00000043  .00000019  .OOOOOOOfl 

.400  .00000032  .0000001*  .00000007 

.420  .00000015  .00000013  .00000007 

.440  .00000005  .00000011  .0000000ft 

.460  .00000000  .00000010  .00000005 

.480  -.00000011  -.OOOnOOOl  -.00000001 

.500  -.00000019  -.00000012  -.0000000* 
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1  MT  1 


DEFLECTIONS  FROM  RESPONSE  FUNCTIONS  ANO  I0»n 

TIME      LOADtLPS.         DEFLErTION    S.IN. 
(SEC)  X*0  X«3.25  X=5.75  X=P.?P 


•  00 

0.0 

0.000000 

0.000000 

0.000000 

•  02 

0.0 

0,000000 

o.oooooo 

0.000000 

.04 

100.0 

•00004P 

.000040 

.OOOOOP 

.06 

150.0 

.000250 

,000104 

.000019 

•  0« 

187. « 

.000270 

.000120 

.000022 

.10 

187.5 

.00028P 

.000128 

.000029 

.1? 

187,  * 

.000293 

.000136 

.000032 

.14 

187.5 

.000300 

.000140 

.000034 

•  1ft 

1P7.C 

.000308 

.000144 

.00003P 

.IP 

187.  P 

•00031P 

.00014P 

.000037 

.20 

187. c 

.000323 

.000152 

.00003P 

.22 

0.0 

.000195 

,0000*0 

.000029 

.24 

0.0 

.000135 

.000056 

•00001P 

•  2ft 

0.0 

.000098 

.OO004P 

•000016 

.28 

0.0 

.000045 

.000040 

•  oooon 

.30 

0.0 

.000038 

.000032 

•000010 

.3? 

0.0 

.000030 

.000024 

•000006 

.34 

0.0 

.000023 

.000016 

.000003 

.36 

0.0 

.000015 

.000008 

.000003 

.38 

o.o 

.000008 

•OOOOOP 

.00000"* 

.40 

o.n 

.OOOOOP 

.000008 

.000003 

.4? 

O.o 

.000008 

•OOOOOP 

.000003 

•  44 

O.o 

.000008 

.000008 

.000003 

.46 

O.o 

.000008 

•000008 

.000003 

.4P 

0.0 

.OOOOOP 

.OOOOOP 

.000003 

.50 

0.0 

•OOOOOP 

.OOOOOP 

. OOOOOP 

]MI? 


RECORD    OF    PAW    DATA 


TIME      LOAD*LPS.         DEFLECTION    S»IN, 
(SEC)  X  =  0  Xe3.25  Xs5.75  X«R.?* 
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i«00 

•  02 

•  04 
.06 

•  on 

.10 

.1? 

•  14 
.16 

•  IP 
.20 
.22 
.24 

•  26 
.28 
.30 
.32 
.34 
.36 
.38 
.40 
.42 
.44 
.46 
.48 
.50 


0 

0 

200 

300 

375 

375 

375 

375 

375 

375 

375 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


.0 

•  0 

■  0 
*0 
.0 

■  0 

•  0 

•  n 

•  0 
«n 
.0 

•  0 
.0 

•  0 
.0 

•  0 

•  0 

.n 

.0 
.0 
.0 

•  0 
.0 
.0 
.0 
.0 


0.000000 
0.000000 
.000120 
.000608 
.000638 
.000653 
.000668 
.000600 
.000713 
.000735 
.000750 
.000413 
.000293 
.000128 
.000098 
.000060 
.000045 
.000030 
.000023 
♦000015 
.000015 
.000015 
.000015 
.000015 
.000015 
.000015 


0,000000 
0.000000 
.000064 
.000320 
.000336 
.000352 
.000368 
.000376 
.000384 
.O003R2 
.000400 
.000256 
.000176 
.000112 
.000080 
.000056 
.000040 
.000024 
.000016 
.000016 
.000016 
.000016 
.000016 
.000016 
.000016 
.000016 


0.000000 

0.000000 

.000003 

,0000*6 

.0000^2 
.000077 

.0000*0 
•0000P2 
.0000*1 
•0000*5 
.0000** 
.000058 
.00002* 
.000016 
.000013 
.000010 
.000006 
.00000^ 
.000003 
.000003 
.000003 
.000003 
.000001 
.000003 
.000001 
.000003 
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1MI? 

RESPONSE    FUNCTIONS    FROM    IMPLICIT    CONVOLUTION 
TIME  RESPONSF      FUNCTIONS 

(SEC)  Xb3.25  X*c.75  X«8.?5 


tOOO 

.00000160 

.00000085 

.00000004 

.020 

,00000798 

,00000420 

.OOOOOOPR 

•  040 

•00000771 

•00000406 

•000000*9 

•  0*0 

.0^000713 

,000o03«* 

.OOOOOOP* 

•  080 

•00000661 

•000003*9 

•00000OA2 

•  100 

•00000623 

•0O000342 

•0000007* 

•  120 

•00000591 

•00000318 

.00000070 

•  140 

•00000560 

•0000029* 

•000000*5 

•  160 

•00000524 

•00000277 

.000000*0 

•  180 

•0O00003T 

•00000065 

.00000017 

•  200 

-•00000057 

-.00000010 

-.00000019 

•  220 

-•00000196 

-•00000055 

-.00000023 

•  240 

-•00000148 

-.0000005=1 

-.0000001* 

•  260 

-.00000116 

-.00000044 

-.0000001 0 

•  280 

-.00000062 

-.00000027 

-.00000007 

.300 

-.00000016 

-.00000013 

-.00000000 

•  320 

.00000034 

.00000007 

.00000004 

.340 

.00000073 

.00000035 

•OOOOOOOR 

.360 

.00000070 

.00000038 

.00000010 

.380 

.00000057 

.00000034 

.00000007 

.400 

.00000032 

.00000025 

•00000004 

.420 

.00000013 

.00000016 

.00000002 

•  440 

.oooooooo 

.00000010 

.00000001 

.460 

-.00000006 

.00000006 

.00000000 

•  480 

-.00000017 

-.00000006 

-.00000002 

.500 

-.00000023 

-.00000016 

-.00000003 
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1MI? 

DEFLECTIONS  FROM  RESPONSF  FUNCTION*:  ANn  LOAD 

TIME      LOADtLRS.         DEFlErTTON    S»IN, 
(SEC)  X=0  X*3.25  **5.75  X«P.PS 


.00 

0.0 

0.000000 

.0? 

0,0 

0,000000 

•  04 

200. ft 

.000120 

•  06 

300,ft 

.00060ft 

•  08 

375,0 

.000638 

.10 

375,0 

.000653 

.12 

375,0 

.000668 

•  14 

375. ft 

•000690 

•  16 

375, ft 

•000713 

.18 

375,0 

.000735 

.20 

375,0 

.000750 

•  2? 

0.0 

.000413 

.24 

0,0 

.000293 

•  26 

0.0 

.00012* 

.28 

0,0 

.000098 

.30 

0,0 

.000060 

.32 

0,0 

.  00004* 

.34 

0*0 

.000030 

•  36 

0*0 

.000023 

.38 

0,0 

.000015 

.40 

O.ft 

.000015 

.42 

0.0 

.000015 

•  44 

0.0 

,000015 

.46 

0.0 

.000015 

.48 

0.0 

.000015 

.50 

0.0 

.000015 

0  .  ft  ft  C  0  0  0 
O.ftftOOOO 
,ftft0064 
•ftft0320 
•ftft0336 
,ftft035? 
,ftftOT68 
•ftft037ft 
,ftft0?84 
•ftft0392 
,ft00400 
,ftft0256 
•ftft0176 
.ftft0112 
.ftftOOPO 
.00005ft 
,000040 
,ftft0024 
.ftftOOlft 
.ftftOOlft 
.ftft0016 
.ftft0016 
,ft00016 
,ftft00)6 
.ftftOOlft 
,ftft0016 


0,000000 
O.OOOftOO 
,000003 
.OOOftftft 
,000072 
.000077 
.OOOftPO 
.0000*2 
.0000*3 
.OOOOP* 
,0000*ft 
.00005* 
,0000?ft 
•00001ft 

•  oooon 

•OOOftlO 
.00000ft 
.00000* 
.000003 
.000003 
.000003 
.000003 
.000003 
.000003 
.000003 
.000003 
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1MI3 


RECORD    OF    RAW    DATA 


TIME      L0AD»LPS.         DEFLErTTON    S.IN, 
(SEC)  X=0  X«3.?5  X=5.75  X=P.?* 


o.ppoooo  0.000000 

0.000000  o.oooooo 

,oooi?P  ,00000* 

,ooo*?4  ,0001*4 

.00065*  .00017? 

.0006PP  .0001^0 

,ooO??0  .0001*4 

.00075?  ,o001»P 

,o007*P  ,0001<3? 

.0007P4  .0001Q* 

.000*00  ,000?00 

.00057*  ,00013* 

,00041ft  .0000*0 

.00025*  .000044 

.oooi"?  .uooo-xs 

,0001?P  .00003? 

.00009*  .000074 

.0000*4  .00001* 

.00004P  .oooonp 

,fl0003?  .00000* 

.00003?  .OOOOOP 

,00003?  .OOOOOP 

.ftO003?  .OOOOOP 

.00003?  .OOOOOP 

.00003?  .OOOOOP 

.^0003?  .OOOOOP 


.00 

0.0 

0.000000 

.0? 

0.0 

0.000000 

.04 

400.0 

,000?40 

.06 

600.0 

.001200 

•  OB 

750.0 

.0012*0 

•  10 

750.0 

.001320 

.1? 

750.0 

.0013P0 

•  14 

750. 0 

.001425 

•  16 

750.0 

.001455 

.IP 

750.0 

.0014P5 

.20 

750.0 

.001500 

.2? 

0.0 

.001110 

.24 

0.0 

,000570 

.26 

O.O 

,000360 

.28 

0.0 

.000270 

•  30 

0.0 

.000210 

.3? 

O.o 

.000150 

.34 

0.0 

.000120 

.36 

0.0 

.000090 

.3fl 

0.0 

.000060 

.40 

0.0 

.00004* 

.4? 

0.0 

.000030 

•  44 

0.0 

.000030 

.4* 

0.0 

.000030 

•  4P 

0.0 

.000030 

.50 

0.0 

.000030 
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1MI3 

RESPONSE  FUNCTIONS  FROM  IMPLICIT  CONVOLUTION 
TIME        RESPONSF   FUNCTIONS 
(SEC)        X-3.25         X«5.75         X»8,25 


1.000  .00000160  .00000085  .00000005 

•0?0  .00000787  .00000400  .00000109 

.040  .00000761  .0000039ft  .00000105 

.060  .0^000724  .00000378  .00000100 

.080  .00000691  .00000361  .00000093 

.100  .000006&1  .00000346  .000000*6 

•120  .00000606  .00000321  .000000«0 

•140  .00000564  .00000299  .00000074 

•160  .00000517  .00000280  .00000070 

•180  .00000220  .00000111  .00000020 

.200  -.00000089  .00000030  -.00000022 

•  220  -.00000150  ~if>O0OO04l  -.00000021 

•240  -.00000124  -.00000044  -.00000014 

.260  -.00000082  -.00000046  -.00000006 

.280  -.000000*8  -.00000028  -.00000002 

.300  -.00000002  -.00000014  .00000001 

.320  .00000036  .00000007  .00000003 

•340  .00000065  .00000024  .00000010 

.360  .00000071  .00000033  .00000011 

•3fl0  .00000045  ,00000033  .00000003 

.400  .00000025  .00000025  .00000005 

.420  .00000010  .00000018  .00000003 

,440  .00000000  .00000012  .00000002 

.460  -.00000005  .00000008  .00000001 

.480  -.00000015  -.00000005  -.00000001 

.500  -.00000020  -.00000015  -.00000004 
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1MT1 

DEFLECTIONS  FROM  RESPONSF  FUNCTION*!  ANn  LOAO 

TIME      LOADtLPS.         PFFLErTTPM    S.IN. 
(SEC)  X=0  X=3.?5  x=S.75  X=8.?n 


l.OO 
.0? 

.06 
.Oft 
.08 
.10 
.1? 
.14 
.16 

•  IP. 
.20 
.2? 

•  ?4 

•  2ft 
.28 

•  30 
.32 
.34 

•  3ft 
.38 
.40 
.4? 
.44 

•  4ft 

•  48 
.50 


0 

0 

400 

600 

750 

750 

750 

750 

750 

750 

750 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


.0 
.0 

•  o 

•  n 

•  o 

.0 
.0 

.0 
.0 

•  0 
.0 

•  0 

.0 

•  o 

*0 
.0 
.0 
.0 
.0 
.0 
.0 
.n 

•  n 
»0 
,o 

.0 


0.000000 
0.000000 
.000240 
.001200 
.0012ft0 
.001320 
.001380 
.0014?S 
.001455 
.001485 
.001500 
.001110 
.000570 
.000360 
.000270 
.000210 
.000150 
.000120 
.000030 
.000060 
.00004*; 
.000030 
.000030 
.000030 
.000030 
.000030 


0.00 
O.ro 

•  no 
,  n  o 
,oo 
.or 

•  no 
,  on 
,oo 
.on 

,PO 

.no 
,no 

,00 
.00 
.00 
.00 
.  00 
.00 
.00 
.00 
.00 
.00 
.00 
.00 
.00 


0000 
0000 
0128 
0624 
0656 
0*88 
07?0 
07^2 
07ft8 

0784 
0800 
0^76 
041ft 
025ft 
01Q2 
0128 
009ft 
0064 
0048 
0032 
0032 
0032 
0032 
0032 
0032 
0032 


0.000000 

o.ooooon 

.000008 
.0001*4 
.000172 
.000180 
.000184 
.000188 
.00019? 

.00019ft 
.000200 
.00013ft 
.0000*0 
,000044 
.00003ft 
.000032 
•0000?4 
.00001ft 
.000008 
.000008 
,000008 
.000008 
.000008 
.000008 
.000008 
.000008 


1HT1 


PFCOPp  OF  RAW  DATA 


TIME   LOAD.LPS.    D  E  F  L  E  r  T  I  0  N  S.IN. 
(SEC)     X=0       X=3.25      X=5.75      X=P,2* 
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0.00 

.o? 

.04 
.Oft 
.08 
.10 

.1? 

.14 
.16 

•  IB 
.20 
.2? 
.24 
.2ft 
.28 
.30 
.3? 
.34 
.3ft 
.3P 

•  40 
.4? 

•  44 
.4ft 
.49 
.50 


0.0 

0.0 

100.0 

170.0 

180. 0 

180.0 

180. 0 

1P0.0 

180.0 

180.0 

180.0 

0.0 

C.n 

0.0 

0.0 

0.0 

0.0 

O.o 

O.o 

0.0 

0.0 

0  ,0 
0.0 
0.0 
0,0 

O.o 


0.000000 

o. oooooo 

.000022 
,00028q 

.000315 
.000360 
.000367 
.000375 
.000390 
.000308 
.00040* 
.000180 
.000060 
.000030 
.000022 
.00001* 
.000007 
.000007 
.000007 
.000007 
.000007 
.000007 
.000007 
.000007 
.000007 
.000007 


0.00 
O.OO 
.00 
.00 
.00 
.00 
.00 
,00 
.00 

.oo 

.00 
.00 
.  00 
.00 
.00 
,00 
.00 
.00 
.00 
.00 
.00 

.00 

.00 
.00 
,00 
.00 


0000 

oooo 

0016 
0124 
0200 
0208 

0212 
0212 
021ft 
0220 
0?24 

oiftfl 

0064 
0040 
0024 
0020 
0016 
0012 
OOOP 
OOOP 
0008 
0008 
OOOP 
OOOP 
OOOP 
0008 


0.000000 

o. oooooo 

,000006 
.OOOOM 
.0000*4 
,000070 
,00007? 
,000074 
.00007* 
.000077 
.000078 

,0000c1 
,0000?6 
.000016 
.000010 
.000008 
.000006 
.000005 
.000003 
.000003 

.oooooo 

,000003 
,000003 

. oooooo 

,000003 
.000003 
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1HI1 


RESPONSE  FUNCTIONS  FROM  IMP|  KIT  CONVOLUTION 
TIMF        RFSPONSF   FUNCTIONS 
(SEC)         X=3.25  X«c;#75  X  =  8.?e; 


.000 
.070 
.040 
.0^0 
.0*0 
.100 

.170 
.140 
.160 
.180 

•  ?00 

•  ??0 
.240 

•  ?60 

•  ?80 
.300 

•  3?0 
.340 
.360 
.380 
.400 
.420 
.440 
.460 
.480 
.500 


.0000006? 
.00000786 

.00000795 
.00000840 
.00000777 
.00000720 
.00000689 
.00000641 
.00000598 
.00000081 
•0O00033? 
.00000304 
.00000210 
.00000132 
.00000067 
.00000009 
.00000073 
.00000126 
.00000105 
.00000061 
.00000025 
.00000001 
,00000012 
.00000018 
.00000026 
.00000026 


.00000044 
.00000340 
.00000519 
.00000490 

.0000045? 
.00000407 
.00000377 
.000003*0 
.00000376 
.0000014? 
.00000178 
.00000131 
•00000114 
.00000068 
.00000031 
.00000001 
.00000073 
.000000*3 
.0000006? 
.00000043 
.00000026 
.0000001? 
.00000004 
.00000000 
.00000011 
•.00000019 


.00000018 
.00000141 
.000001*3 
.000001*4 
.000001^2 
.0000014? 
.0000013? 
.00000173 
•0000011* 
•  000000-sO 

-.oooooo-n 

-.000000^8 

-.000000^6 

-.00000072 

-.00000010 

.00000000 

.00000008 

.00000019 

.00000020 

.0000001* 

.00000010 

,00000005 

.0000000? 

.00000001 

-.00000003 

-.00000007 
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1HI1 


DFFLFCTTONS  FROM  RFSPONSF  FUNCTIONS  ANn  LOAD 

TIME      LOADERS.         DFFLErTTOM    S»IN, 
(SFC)  XsO  X=3.?5  X=5.75  ***.?* 


0.00 

0.0 

0.000000 

0. 000000 

0.000000 

.0? 

0.0 

0.000000 

0.000000 

0.000000 

.04 

100. o 

.oooo?? 

.00001* 

.000006 

.06 

170. o 

,000?B5 

,0001?4 

.ooooci 

•  OR 

lftO.o 

.00031? 

,ono?oo 

.0000*4 

.10 

1R0.0 

.000360 

,000?OP 

.000070 

.1? 

IflO.O 

.000367 

.noO?1? 

.00007? 

.14 

iflo.n 

.00037^ 

,noo?l? 

.000074 

.16 

1R0.O 

.000390 

,0O0?1* 

•00007R 

•  ]fl 

IflO.O 

.000397 

,000??0 

.000077 

•?0 

IflO.O 

.00040* 

.P00??4 

,00007ft 

.?? 

0.0 

.0001PO 

.P0016R 

.000051 

.?4 

0.0 

.000060 

.000064 

,0000?6 

•  ?6 

0.0 

.000030 

.000040 

•0P0O16 

.?« 

o.o 

,0000?3 

,OO00?4 

.000010 

.30 

0.0 

.000015 

,0000?0 

.OOOOOR 

.3? 

0.0 

.000007 

.000016 

.00000* 

.34 

O.o 

.00000ft 

.0000)? 

.000005 

.3* 

0.0 

.000007 

.OOOOOR 

.000003 

.3H 

O.n 

.000007 

,00000ft 

,000003 

.40 

0.0 

.000007 

,00000ft 

.000003 

.4? 

0.0 

.000007 

,000008 

.000003 

.44 

O.o 

.000007 

,00000ft 

.000003 

.46 

0.0 

.000007 

.OOOOOR 

.000003 

.4ft 

0.0 

.OOOOO8 

.OOOOOR 

. OOOOOR 

.50 

0.0 

.000008 

,00000ft 

. OOOOOR 
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1HT? 

RECORD    OF    PAW    OATA 

TIME       LOAD.LPS.  OFFLtrTTON    S,TN. 

(SEC)  X  =  0  X*3.25  x=*.7*  X=P.2* 


>.00 
.0? 

.04 
.06 
.OP 
.10 
.1? 
.14 
.16 
.18 

•  ?0 
.?? 
.?4 
.26 
.2* 
.30 
.3? 

•  34 
.36 
.3* 
.40 
.4? 
.44 
.46 

•  4R 
.50 


0.0 

O.n 

100.0 

300.0 

375.0 

375.0 

375.0 

375.0 

375.0 

375.0 

375.0 

375.0 

O.o 

O.o 

0.0 

Co 

O.o 

0.0 

o.o 

0.0 
0.0 
0.0 
0.0 

o.o 

0.0 

o.o 


0.000000 
0.000000 

.000030 
.0004PB 
.000*25 
.0010*0 
.00109* 
.001110 
.001125 
.001140 
.0011*5 
.001170 
.000600 
.000210 
.000150 
.000135 
.00012* 
.000120 
.00010* 
.000090 
.0000*3 
.000075 
.00007* 
.000075 
.00007* 
.000075 


0.00 
0.00 

•  nO 
,00 
.00 
.00 
.00 

•  OP 
.00 
.00 
,00 
.00 
.00 
.00 

,00 

.00 
.00 
,00 
.00 
.00 
.00 
.00 
.00 
.00 
,00 
.00 


0000 
0000 
0040 
0400 
0*20 
0*60 
0*76 
0**4 
0*9? 
0*00 
060* 
0616 
0440 
0216 
0144 
0120 
0112 
010* 
0104 
0100 
0096 
009? 

oopp 

00*4 
00*0 
00*0 


0.000000 

o.oooooo 

.00001* 
.0001?* 
.00017ft 
.0001 92 
.00019« 
.000200 
.000201 
.00020* 
.000210 
,000?1 1 
.0000P0 
,00004* 
.000045 

.00004? 
.000040 

.ooooiR 

. 00003* 
,000034 

,oooo*»? 

.00003? 
♦00003? 
.OOOOT? 
.00003? 
.00003? 
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1HT? 

RESPONSE    FUNCTIONS    FROM    IMPLICIT    CONVOLUTION 
TIME  RFSPONSF       rUNCTlONS 

(SEC)  X=3.?5  X=k.7?  X=8,?* 


.000 

•00000040 

.00000053 

.000000?! 

.n?o 

.00000647 

,000o05?o 

.ooooni*<5 

.040 

.00001044 

.00000*46 

.00000?19 

.060 

,ooooi?48 

.00000637 

.00000?19 

•  080 

.00001187 

.00000594 

.00000?06 

.100 

.00001087 

.00000545 

.00000187 

•  1?0 

.00000996 

.00000500 

.0000017? 

•  140 

.00000915 

.0*000460 

•00000161 

•  160 

•0000084? 

.00000424 

.00000148 

•  1«0 

.00000776 

.00000391 

.000001 34 

•  ?oo 

-.00000059 

•000001?2 

-.0000005? 

.??o 

-.00000525 

-.00000141 

-. 0000007ft 

.?40 

-.00000457 

-.000001(6,4 

-.OOOOOOM 

•  ?60 

-.00000305 

-.000001  16 

-.000000?8 

•  ?80 

-.000001&? 

-.00000055 

-.00000007 

.300 

-.00000045 

.00000001 

.00000010 

.3?0 

.00000041 

.00000047 

.000000?3 

.340 

.00000110 

.000000*3 

.00000036 

.360 

,000001 76 

.000001 13 

.00000065 

.380 

•00000?26 

.00000136 

.00000054 

.400 

.00000199 

.00000130 

.00000043 

.4?0 

,00000126 

.O0OOO0Q7 

.00000031 

.440 

.00000066 

.00000065 

•000000?3 

•  460 

•00000028 

.00000046 

.00000017 

•  4fl0 

-.0000004? 

-.00000018 

-.OOOOOOOft 

.500 

-.00000094 

-.00000070 

-,000000?7 
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1HT2 

DEFLECTION*;  FROM  RFSPONSF  FUNCTIONS  AND  I  OAP 

TIME      LOAD. LPS.         D'EFLEpTTON    S.IN. 
fSEC)  X  =  0  X*3.25  Xe5.75  Xs8.?e; 


.00  0.0  0.000000  0.000000  o.oooooo 

.02  0.0  0.000000  0.000000  0.000000 

.04  100.0  .000030  .000040  .00001* 

•Oft  300.0  .0004RR  .000400  .0001?* 

•  OR  375.0  .000*2*  .000*20  .00017* 

.10  375.0  .001050  .000*60  .000192 

.1?  375.0  .001095  .000*176  .0001QP 

.14  375.0  .001110  .0005*4  .000200 

.16  375.0  .001125  .000592  .000203 

•lfl  375.0  .001140  .000*00  .00020* 

.20  375.0  .0011**  .00060*  .000210 

,2?  375.0  .001170  .000616  .000?11 

.24  0.0  .000600  .000440  .0000*0 

.26  0.0  .000210  .00021*  .00004* 

.2*  0.0  .000150  .000144  .00004* 

.30  O.O  .00013*  .000120  .000042 

.32  0.0  .00012*  .000112  .000040 

.34  0.0  .000120  .00010*  .0000"^* 

.36  0.0  .000105  .000104  .0000?* 

.3*  0.0  .000090  .000100  .000034 

.40  0.0  .0000*3  .000096  .0000*J2 

.42  0.0  .000075  .000092  .0000?? 

.44  0.0  .00007*5  .0000**  .000032 

.46  0.0  .000075  .0000*4  .000032 

.4*  0.0  .000075  .0000*0  .000032 

.50  0.0  .000075  .0000*0  .000032 
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iHn 


RECORD    OF    PAW    DATA 


TIME      L0AD»LRS.         DFFLErTION    S.IN. 
(SEC)  X  =  0  X=3.25  *  =  5. 75  X=8.?5 


0.00 

0.0 

0.000000 

0.000000 

0.000000 

.0? 

0.0 

0.000000 

0.000000 

0.000000 

.04 

3^0.0 

.O00OQ0 

.000080 

,000040 

.06 

600.  0 

.001500 

.000640 

.000224 

•  on 

750.0 

.001710 

.000960 

,000344 

.10 

7^0.0 

.0019^0 

.001 0*0 

•00038* 

.1? 

750.0 

.002040 

.001072 

.00040* 

.14 

750.0 

.002100 

.001 1 04 

.000412 

.16 

750.0 

.002115 

.001120 

.00041ft 

.1* 

750.0 

.002130 

.0011*2 

,000424 

#20 

750.0 

.002160 

.001  1(S8 

,00042* 

.22 

7*0.0 

.00217^ 

.0011*4 

,00043? 

.24 

0.0 

.001200 

.000480 

.000320 

.26 

0.0 

.000450 

.000280 

.000160 

.28 

O.O 

.000240 

.000208 

.000112 

.10 

0.0 

.000210 

.000176 

.0000*0 

.1? 

0.0 

.000180 

.0001*0 

.000076 

.14 

0.0 

.000165 

.000152 

.00007? 

.36 

0.0 

.000150 

.000144 

.000064 

•  3R 

0.0 

.000120 

.000136 

.000060 

.40 

O.O 

.000105 

.000120 

.0000** 

.4? 

0.0 

.000090 

.000112 

.0000*2 

.44 

0.0 

.000090 

.000104 

.000048 

.46 

0.0 

.000090 

.000096 

.000044 

.48 

0.0 

.000090 

.OO008R 

.000044 

.50 

0.0 

.000090 

,n0008« 

,000044 
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RESPONSE    FUNrTIONS    FROM    TmPiKTT    COMVOLllTION 
TIME  RESPONSE       FUNCTIONS 

(SEC)  Xs3.2*  X-c,7R  Xs8.?s 


.000 

.00000060 

.00000053 

.00000077 

.0?0 

.00000995 

. 000004?? 

.00000147 

.040 

•000010^* 

.00000601 

.00000715 

.060 

.000011  10 

.000O0S9P 

.00000774 

.080 

.00001060 

.00000559 

.00000715 

.100 

,00000993 

.00000«^?4 

.00000196 

.1?0 

•  00000*503 

.00000481 

.00000179 

#140 

•0O000S20 

•000004S4 

.000001*6 

.160 

.0*000757 

.00000418 

.000001*5? 

•  180 

•00000690 

.000^0387 

.000001 79 

•  POO 

-.00000024 

-.000001 17 

.000000*3 

.??o 

-.00000437 

-.00000?06 

-.0000004* 

.?40 

-.0^000436 

-.00000175 

-, 000000*4 

•  ?60 

-.00000301 

-.000001  10 

-.00000048 

•  ?80 

-.0000018? 

-.00000060 

-,000000?4 

.300 

-.00000073 

-.ooooroo* 

-.00000004 

.3?0 

.00000017 

.00000038 

.00000009 

.340 

.00000079 

.0000007* 

.00000073 

.360 

.00000138 

.000001 00 

.00000033 

.380 

.00000185 

•000o01?4 

.00000041 

.400 

.00000165 

.00000095 

.00000040 

.4?0 

.00000104 

.0000005Q 

,00000079 

.440 

.000000*9 

.oonooo?^ 

.00000070 

.460 

.00000013 

.oooooou 

.00000013 

.4R0 

-.00000037 

-.000000?3 

-.00000005 

.500 

-.0000007? 

-.0000005? 

-.00000070 
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DFFIFCTIONS  FROM  RFSPOMSF  FUNCTIONS  ANn  LOAO 

TIMF   LOAD»LPS.    D  F  F  I  E  r  T  I  0  N  S.IN. 
(SFC)      X  =  0       X*3.?5      y  =  c;.75      X  =  8.?e; 


0.00 

0.0 

0.000000 

0.000000 

0.000000 

.0? 

o.o 

o.oooooo 

0.000000 

o.oooooo 

.04 

350. o 

.000090 

. OO00R0 

.000040 

.06 

600.0 

.001500 

.000*40 

,000??4 

.09 

750.0 

.001710 

.0009*0 

.000144 

.10 

750.  0 

.001950 

.001040 

.000388 

.1? 

750.0 

.002040 

.00107? 

.000408 

.14 

750.  0 

.002100 

.001104 

.00041? 

.16 

750. o 

.002115 

.001120 

,00041* 

.18 

750.0 

.002130 

.0011=;? 

.000424 

,20 

750.0 

.  0021*0 

.0011*8 

,0004?8 

.?? 

750.0 

.00217^ 

.001184 

.0004^? 

.?4 

0.0 

.001200 

.0004^0 

. 000320 

.26 

0.0 

.000450 

,OO0?8O 

.0001*0 

•  ?fl 

0.0 

.000240 

,000?08 

.00011? 

.30 

0.0 

.000210 

.000176 

.OOOOPO 

.3? 

o.o 

.000180 

.0001*0 

.00007* 

.14 

0.0 

.00016=; 

.00015? 

.00007? 

.36 

o.o 

.0001^0 

.000144 

.0000*4 

.38 

o.n 

.000120 

.00013* 

.OOOOf 0 

.40 

0.0 

.00010^ 

.000120 

.onooc* 

.4? 

0.0 

.000090 

.00011? 

.00005? 

•  44 

O.o 

.000090 

.000104 

.000048 

.46 

o.o 

.OO0OQ0 

. orooo* 

,000044 

.48 

0.0 

.000090 

.OOOOBR 

.000044 

•  50 

0.0 

.000090 

. 000088 

.000044 

?LI1 


RFCOPf)    OF    PAW    HATA 


TIMF       LOADtLPS.  DFFLFrTTO*'    S.IN. 

(SEC)  X  =  0  X  =  3.?«^  X  =  R.7C.  x=8.?s 
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0.00 

0.0 

0.000000 

o.oooooo 

0.000000 

.0? 

O.P 

0.000000 

0. oooooo 

0.000000 

.04 

100.0 

.00001 1 

.000010 

.000003 

.06 

150.0 

.000045 

.00003? 

.00000ft 

.08 

175.  n 

.000064 

.000034 

.000008 

.10 

107.  c 

.000068 

.oooo?* 

.00001 0 

.1? 

187.* 

•00006P 

.0000^7 

.000011 

.14 

187. S 

.000071 

.000038 

.00001? 

.1ft 

187.* 

.000071 

.000040 

.000014 

.18 

187. c 

.00007S 

.00004? 

.  c;  o  o  o  i  ft 

.20 

187. ^i 

. 00007^ 

.00004? 

.000018 

.22 

O.o 

.00004^ 

. 00002*5 

.000008 

.?4 

0.0 

.000038 

.00001ft 

.OOOOOA 

.26 

0.0 

.000030 

.000010 

.000005 

.28 

O.P 

.00002? 

. oooooft 

.00000? 

.30 

O.o 

.00001^ 

.OOOOO* 

.ooooo? 

.3? 

0.0 

.00001 1 

.000003 

.ooooo? 

.34 

0.0 

.000008 

.00000? 

.00000? 

.36 

0.0 

.000004 

.ooooo? 

.00000? 

.38 

O.o 

.000004 

,00000? 

.ooooo? 

.40 

0.0 

.000004 

.ooooo? 

.ooooo? 

.42 

0.0 

.000004 

.00000? 

.00000? 

.44 

0.0 

.000004 

.00000? 

.ooooo? 

.4ft 

0.0 

.000004 

.00000? 

.ooooo? 

.48 

0.0 

.000004 

.OOOOO? 

.00000? 

.so 

0.0 

.000004 

,00000? 

.ooooo? 
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2LT1 


RESPONSE    FUNCTIONS   FROM    IMPLICIT    COMVOl  iiTTON 
TIME  RESPONSE       r    UNCTIONS 

(SEC)  X=3.25  X=^.7S  X=B,?* 


.000 

.0^000029 

.00000027 

•OOOOOOOP 

•  o?o 

,0^0001  IP 

.OQ0OO0P3 

,00000015 

.040 

.00000159 

.OOOOOOR? 

.00000019 

.060 

.00000155 

.0000007* 

.000000?3 

.ORO 

.00000139 

.00000074 

.00000OP3 

•  TOO 

•  0«0001-*3 

.00000069 

.00OO0O?6 

•  1?0 

.00000120 

.00000067 

.00000026 

•  140 

.000001  IB 

.00^00066 

.00000029 

•  160 

•00000106 

•oOoooosq 

.00000032 

.190 

.0000001H 

.OOOOOOPl 

,00000003 

.200 

•ooooooop 

-.OO00OOO8 

-.00000002 

•  220 

•oooooooi 

-.0000001 6 

-.00000003 

.240 

-.00000003 

-.0000001P 

-.00000006 

•  ?60 

-.ooooooio 

-•oOooooi 1 

-.00000005 

•  2B0 

-.00000007 

-.OOOOOOOP 

-.00000002 

.300 

-.0^000002 

-•OOOflOOO^ 

.00000001 

.320 

-.00000001 

.00000004 

.00000004 

.340 

.ooooooio 

.00000009 

.00000007 

•  3(S0 

.00000011 

•  OOOOOOI  1 

,0000000ft 

.380 

•ooooooii 

•roonoooQ 

,00000005 

•  400 

•ooooooio 

.00000006 

.00000005 

.4?0 

.00000009 

.OOOP0004 

.00000004 

.440 

.00000007 

.00000002 

.00000003 

•  460 

.00000005 

.OOOOOOOI 

.00000002 

•  4B0 

-.OOOOOOOI 

-.00000002 

-.OOOOOOOI 

•  soo 

-.00000006 

-.00000004 

-.00000003 
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PL  T  1 


DFFLFCTTONC  FROM  PFSPONSF  FUNCTIONS  ANO  1040 

TIMF       LOAD, IPS.  DFFL     ErTTOM    S,IN, 

(SFC)  X=0  X=3.?5  X*5.75  X=P.?R 


0.00 

0.0 

0.000000 

0.000000 

0.000000 

.0? 

o.o 

0.000000 

O.nOOOOO 

0.000000 

.04 

100,0 

.00001 1 

.000010 

.000003 

•  Oft 

150.0 

.000045 

.00003? 

.00000ft 

•  OR 

175.0 

.000064 

,000034 

.OOOOOR 

.10 

187. 5 

.OOOOftfl 

.000035 

.000010 

.1? 

1R7.5 

.OOOOftp 

.000037 

.000011 

.14 

187. «! 

.000071 

.00003R 

.oooon 

.1ft 

187. c 

.000071 

.000040 

.000014 

.18 

187.  e, 

.00007R. 

.00004? 

.00001ft 

•  20 

187. R 

.00007* 

.00004? 

.000018 

•  ?? 

O.n 

.000045 

.O000?9 

.OOOOOR 

.?4 

0.0 

.00003* 

.00001ft 

.oooooft 

•  ?6 

o.o 

.000030 

,00001 0 

.OOOOOS 

.?R 

0.0 

•  oooo?-* 

.oooonft 

.OOOOOT 

.30 

0.0 

.000015 

.000005 

.00000? 

.3? 

0.0 

.000011 

.000003 

.00000? 

.34 

0.0 

.000008 

.00000? 

.00000? 

.3ft 

0.0 

.000004 

.00000? 

.00000? 

.38 

0.0 

.000004 

,0^000? 

.00000? 

•  40 

0.0 

.000004 

.00000? 

.00000? 

.4? 

0.0 

.000004 

.ooooo? 

.00000? 

•  44 

O.o 

.000004 

.00000? 

.00000? 

•  4ft 

0.0 

.000004 

.00000? 

.00000? 

•  4R 

0.0 

.000004 

.OOOOO? 

.ooooo? 

•  SO 

0,0 

.000004 

.00000? 

.00000? 

?LT? 


RFCOPD  OF  PAW  DATA 


TIME"       LOADtLPS.         DFFLFrTTPN    StIN, 

(SEC)  x  =  o  x=3.?*  x=5.75  x*R.?* 
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.00 

0.0 

0.000000 

0.000000 

0.000000 

.0? 

0.0 

0.000000 

0.000000 

0.000000 

.04 

200. n 

.000023 

.ooOOOft 

.00000* 

.06 

350.0 

.000045 

.00001 0 

,00000ft 

•  OR 

375.0 

•ooooon 

.00003P 

.000010 

.10 

375.0 

.000113 

.000054 

.000016 

.1? 

375.0 

.000124 

,000070 

.oooo?i 

.14 

375.0 

.0001?fl 

.oooopo 

.0000?4 

.1ft 

375.  o 

.000131 

.000083 

.000056 

.1H 

375.0 

.000139 

.OOOOPO 

,0000?7 

.20 

375.0 

.000143 

.OOOOPft 

,0000?9 

.2? 

0.0 

.000090 

. 0000ft7 

.0000?? 

.?4 

0.0 

.0000P3 

.000054 

.000019 

.?ft 

O.P 

•  00004e; 

.00002ft 

.OOOOOP 

.2fl 

O.o 

.000030 

.00001° 

.00000* 

.30 

0.0 

.0000?ft 

.000013 

.000007 

.3? 

0.0 

.000023 

.000010 

.00000? 

.34 

0.0 

.000019 

.OOOOOP 

.00000? 

.3ft 

o.o 

.000015 

.00000^ 

.00000? 

.3P 

0.0 

.00001 1 

.00000* 

.00000? 

.40 

O.o 

.OOOOOP 

.000003 

.00000? 

.4? 

0.0 

.OOOOOP 

.00000? 

.00000? 

.44 

o.o 

.OOOOOP 

.000002 

.00000? 

•  4ft 

O.o 

.OOOOOP 

.00000? 

.ooooo? 

.4P 

0,0 

.ooooop 

.00000? 

.ooooo? 

.50 

0.0 

.OOOOOP 

.000002 

.00000? 
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2LT? 

RFSPONSF    FUNCTIONS    FROM     IMP|   TTTT    CPf»VOl.l  |TI  OK 
TIME  RESPONSF      FUNCTIONS 

(SEC)  XrT.?c  X=c.75  X=P.?* 


.000 

.00000031 

.00000008 

.00000007 

•  o?o 

.00000057 

.OOOnOOl 3 

.00000007 

.040 

•onpooi 1? 

.00000049 

.0000001 ? 

.0*0 

.00000131 

.oooooo*^ 

,0000001 q 

.080 

.00000133 

.000OO0B0 

,000000?4 

.100 

.00000125 

.00000086 

.000000?5 

•  l?n 

.000001  17 

.00000081 

.000  0  0  075 

•  140 

.000001 16 

•00000076 

.00000074 

.160 

.00000109 

.00000069 

•000000P4 

.180 

.00000031 

.0000003* 

.0000001 3 

•  200 

•00000024 

•00000018 

•00000009 

•  ??0 

-.00000018 

-.00000017 

-.00000006 

.?40 

-.00000024 

-.00000018 

-.00000007 

•  260 

-.00000013 

-.00000016 

-.00000007 

•  ?«0 

-.00000003 

-.00000010 

-.OOOOOOoe, 

.300 

.00000003 

-,00000003 

-.0000000? 

•  3?0 

•OOOOOOOQ 

.00000002 

.00000001 

.340 

.00000014 

.00000007 

.0000000? 

.360 

.00000012 

.00000008 

.00000004 

•  3P0 

.00000013 

.00000007 

. OOOOOOOE 

.400 

.00000010 

.OOOOOOOQ 

.00000004 

.4?0 

.00000006 

.00000003 

.0000000? 

•  440 

.00000004 

•OOOoOOOl 

.0000000? 

.460 

.00000004 

-.00000000 

.00000001 

.480 

-.00000002 

-.0000000? 

-.00000001 

.500 

-.00000006 

-.oonooooi 

-.0000000? 
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?LI? 

DEFLECTIONS  FROM  RESPONSE  FUNCTIONS  ANn,  LOAD 

TIME       LOADtlPS.  DEFLErTTON    S.IN. 

(SEC)  X  =  0  X=3.??  Y  =  5.7c,  X»R.?F 


0.00 

0.0 

0.000000 

0. oooooo 

o. oooooo 

.0? 

0.0 

0.000000 

o.onoooo 

0.000000 

.04 

?00.o 

.0000?3 

.00000* 

.oooon* 

.06 

350. n 

.00004S 

.00001 0 

.00000* 

.OR 

375.0 

.000090 

.00003P 

.000010 

.10 

375.0 

.0001)3 

.0000^4 

.000016 

.1? 

375.0 

.0001?4 

,oroo70 

,0000?1 

.14 

375.0 

.0001?* 

.0000*0 

.000024 

.16 

375.0 

.000131 

.OO00R3 

.0O0O?6 

•  lfl 

375.0 

.000139 

.OOOOP^ 

,0000?7 

#20 

375.0 

.000143 

.000086 

.000079 

.?? 

CO 

.000090 

.000067 

.0000?? 

•  ?4 

0.0 

.0000B3 

.0000^4 

.000019 

•  ?ft 

0.0 

.000045 

.000026 

.OOOOOP 

.28 

0.0 

.000030 

.000019 

.oooooo 

.30 

O.o 

.0000?6 

.000013 

.000001 

.3? 

O.o 

.O000P3 

.000010 

.00000? 

.34 

0.0 

.000019 

.OPOOOfl 

.00000? 

.36 

0.0 

.000015 

.000006 

.00000? 

.38 

0.0 

.000011 

.opooo* 

.00000? 

.40 

0.0 

.000O0R 

.000003 

.ooono? 

.42 

0.0 

.000008 

.00000? 

.00000? 

.44 

0.0 

.000008 

.00000? 

.00000? 

.46 

0.0 

•OOOOOP 

,00000? 

,00000? 

•  48 

o.o 

.000008 

.00000? 

.00000? 

.50 

0,0 

.000008 

.00000? 

.00000? 

?i_n 


RECORD    OF    PAW    OATA 


TIMF       LOAD. LPS.  DEFLEfTTON    S.IN. 

(SEC)  X=0  X=3.?5  Xe^.lS  *=*.?* 
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.00 

0.0 

0.000000 

0.000000 

0.000000 

.0? 

0.0 

0.000000 

o.nooooo 

0.000000 

.04 

400.0 

.00007S 

.oooo?4 

.000040 

.06 

700.0 

.000150 

.ooooPB 

.oooopo 

•  OR 

750.0 

•  00037«; 

,non?4P 

.000120 

•  10 

7R0.0 

.00043^ 

.n00304 

.0001?4 

.1? 

750.0 

.000443 

,noo3?4 

.0001?P 

.14 

750.0 

.000450 

,nn03?p 

.00013? 

.16 

7S0.O 

.000465 

,ooo3?P 

.0001 3? 

.IP 

750.0 

.0004PO 

.n0033? 

,00013ft 

•  ?C 

750.0 

.0004PP 

,00033ft 

.00013ft 

.?? 

0.0 

.000??^ 

,roo?4fl 

.OOOOPO 

.?4 

0.0 

.oooios 

.000064 

,0000?& 

.?ft 

o.n 

.00007=; 

,0000?0 

.OOOOOP 

.?R 

0.0 

.000060 

.00001? 

,000004 

.30 

O.n 

.00004F 

.ooooop 

.000004 

.3? 

0.0 

.00003P 

,000004 

.000004 

.34 

0.0 

.000030 

.000004 

,000004 

.36 

o.n 

.0000?3 

.000004 

.000004 

.3P 

0.0 

.oooom 

.000004 

.000004 

.40 

o.o 

.000015 

.000004 

.000004 

.4? 

0.0 

.ooooi*; 

.000004 

.000004 

.44 

o.o 

•  ooooi6; 

.000004 

.000004 

.4* 

0.0 

.00001^ 

.000004 

.000004 

.48 

o.r 

.000015 

.000004 

.000nn4 

.50 

0.0 

•  ooooi6; 

.000004 

.000004 
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2U3 

PESPONSF    FUNCTIONS    FROM    IMP|   TTTT    CPK'VOl  llTlON 
TIMF  RFSPONSF       FUNCTIONS 

(SEC)  xei.?';  x=«=;.75  xsR.?6; 


.000 

.00000050 

.  0  0  0  n  0  0  1  6 

.00OO0OP7 

.020 

.00000095 

,ooooooe;7 

.000000*1 

.040 

.0^000236 

.0ftftftPl*.8 

.00000073 

•  060 

.0ft000253 

.00000181 

.000000*8 

.OflO 

.00000234 

.00000176 

.000000*4 

.100 

.00000215 

.00000161 

.000000*0 

.120 

.00000203 

.OOOOOK^ 

.000000=14 

•  140 

.00000193 

•  0  0  0  ft  0  1  J  3 

•OOOOOOci 

.160 

.00000179 

.0000012? 

.00000046 

•  180 

-.00000010 

.000000=13 

.00000007 

•  200 

-.00000080 

-.00000070 

-.00000026 

•  220 

-.00000069 

-.00ftft0077 

-♦000000?7 

•  ?40 

-.0^000047 

-.00000057 

-.oooooopo 

•  260 

-.00000029 

-.00000036 

-.0000001 ? 

•  PRO 

-.00000009 

-.00000019 

-.0000000=, 

.300 

.00000007 

-.0000000? 

.00000001 

•  3?0 

.00000021 

.OOOftOOl 1 

.000000ft* 

•  340 

•00000032 

•oOnooo?? 

.0000001 0 

.360 

.00000028 

.00000026 

.0000001 0 

.3P0 

•0*000017 

•000ft00l6 

.00000006 

.400 

.00000008 

.00000007 

.00000003 

•  4?0 

.00000003 

•OOOftOOOO 

.00000001 

•  440 

-.00000001 

-.000*0003 

-.00000001 

.460 

-.00000001 

-,00OftO005 

-.00000001 

.480 

-.00000006 

-.00000006 

-.00000002 

.500 

-.00000008 

-.00000006 

-.00000003 
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?Ln 

DEFLECTIONS  FROM  RFSPONSF  FUNCTIONS  AND  Loan 

TIMF       LOAD»l  PS.  DFFI      ErTTON    S.IN. 

(SEC)  X=o  X«3.?5  Y=5.7*  x=R.?P 


0.00 

0.0 

0.000000 

o.oooooo 

0.000000 

.0? 

o.o 

0.000000 

0.000000 

0.000000 

.04 

400,  n 

.000075 

,n000?4 

.000040 

.Of 

700.0 

.000150 

.OOOOPR 

.OOOOPO 

.00 

750. 0 

.00037* 

,000?4P 

.0001 ?0 

.10 

7S0.O 

.00043* 

.000304 

.00O1P4 

.1? 

750.0 

.000443 

,0003?* 

.0001PR 

.14 

750.0 

.000450 

,0003?P 

.00013? 

.1* 

750.0 

. 0004*5 

,0003?P 

.00013? 

•  1R 

750.0 

.0004P0 

.00033? 

.000136 

.20 

750.0 

.0004PP 

.  000336 

,0001 36 

.?? 

0.0 

,ooo??* 

#O00?4R 

.OOOORO 

.24 

0.0 

.00010* 

.000064 

,0000?4 

•  ?6 

O.o 

.00007c; 

.oooo?0 

,ooooo« 

.?P 

0.0 

.000060 

.00001? 

,000004 

.30 

0.0 

.000045 

.OOOOOR 

.000004 

.3? 

O.O 

.ttO003P 

.000004 

.000004 

.34 

O.o 

.000030 

.000004 

.000004 

.36 

0.0 

.000<i?3 

.000004 

.000004 

.38 

0.0 

.000015 

.000004 

.000004 

.40 

0.0 

.000015 

.000004 

.000004 

.4? 

0.0 

,000015 

.000004 

.000004 

•  44 

0.0 

.000015 

.000004 

.000004 

.46 

0.0 

.00001* 

.000004 

.000004 

.4fl 

0.0 

.00001* 

.000004 

.000004 

.50 

0  ,n 

.000015 

.000004 

.000004 

2MT1 


RECORD    OF    PAW    DATA 


TIMF      LOAD.LPS.         OFFI     ErTIOM    S»IN. 
(SEC)  X  =  0  X*3.?5  y=5.75  X  =  fl.;?5 
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.00 

O.o 

0.000000 

0.000000 

0.000000 

.0? 

0.0 

0.000000 

0.000000 

0.000000 

.04 

100. n 

.000015 

.00001 2 

.00000* 

.06 

175.0 

.00010] 

.000072 

.000041 

.OP 

107. c 

.00013* 

.000108 

,ooooa« 

.10 

1R7.«; 

.000157 

.000120 

.0000*6 

.12 

lft7.e; 

.000167 

.000124 

.OOOOfP 

.14 

1P7.«? 

.000180 

.00012P 

.000070 

•  If 

187.  5 

.OOOI84 

.00011? 

.000070 

.lfl 

187. e: 

.00018R 

.0001 36 

.000072 

.20 

187. q 

.00019! 

.000140 

.000074 

.?? 

O.o 

.00013* 

.oroioft 

.0000*8 

.?4 

0.  n 

.000101 

.000072 

.000043 

.26 

0.0 

.000067 

.000032 

.000020 

•  2ft 

0.0 

.000030 

.000024 

.000013 

.30 

0.0 

.000022 

.00001 6 

.000010 

.32 

0.0 

.000015 

•OOOOOR 

.000006 

.34 

0.0 

.000007 

.000004 

.000003 

.36 

0.0 

.000007 

.000004 

.000003 

.3ft 

0.0 

.000007 

.000004 

.000003 

.40 

0.0 

.000007 

.000004 

.000003 

.42 

0.0 

.000007 

.000004 

.000003 

.44 

0.0 

.000007 

.000004 

.000001 

.46 

0.0 

.000007 

.000004 

.000003 

•  4R 

0.0 

.000007 

.000004 

.000003 

.50 

0.0 

.000007 

.000004 

.000003 

2MI] 


178 


RESPONSE  FUNCTIONS  FROM  IMP|  KIT  CONVOLUTION 


RESPONSF   FUNCTIONS 
Xs3»25        ^a*^^5^        **ft.25 


0*000 
.020 

•  040 

•  060 

•  080 

•  1-00 
.120 
.140 

•  1*0 

•  ISO 
.200 
*220 

•  240 

•  260 
.280 
.300 
.320 

•  340 

•  360 
.380 

^400 

•  420 

•  440 

•  4*0 
.480 
.500 


•00000040 

.00000266 

.00000331 

.00000358 

.00000348 

*0OOOO348 

.00000323 

.00000301 

.0000028] 

.00000106 

.00000031 

-  .00000029 

-.00000091 

-.00000068 

-.00000046 

-.00000029 

.00000004 

r00000032 

.00000040 

.00000039 

.00000032 

•00000020 

.00000011 

.0^000005 

-.00000008 

-.00000017 


.00000032 
.00000189 
r000O0267 


.00000273 
.00000256 

TfKmtre-24-r 

.00000228 

.00000216 

.0000020** 

.00000102 

.00000014 

-.-001^00068 

••00000056 

-.00000046 

-.00000058 

-.00000022 

.00000002 

.00000022 

.00000030 

.00000029 

.000*00^9 

•00000011 
.00000006 
r00O*00frl 
.00000007 
.00000011 


.00000013 

.000001 14 

.^0000142 

.00000149 

,00000143 

.00000135 

.000001 19 

.000001 12 

.00000105 

.000000*3 

.00000020 

-.00000006 

-.00000014 

-.00000025 

-.00000017 

-.00000012 

.00000000 

-T00  0000+t 

.00000015 

.00000016 

r00O0OO13 

.00000009 

.00000005 

r0000O003 

-.00000003 

-.00000007 
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?MTl 


DEFLECTIONS  FROM  PFSPONSF  FUNCTIONS  ANr,  !  HAD 

TIMF   LOADttPS.    OFF!  E  C  T  I  0  M  S,IN. 
(SEC)     X  =  0       X=3.25      YrS.75      X  =  8.?e; 


.00 

O.n 

0. oooooo 

o.nnoooo 

0  .000000 

.0? 

0.0 

0.000000 

O.nnonoo 

0.000000 

.04 

100.0 

.000015 

.00001 ? 

.oooons 

.0* 

175. n 

.oooioi 

.nnon7? 

.00004"? 

•  OR 

187, c 

.000135 

.nnoiOR 

.OOOO^P 

.10 

187. c 

.000157 

,nooi?0 

,000066 

.1? 

187, c 

.000167 

,onoi?A 

.000069 

.14 

187.  * 

.000180 

.0001?8 

.000070 

.If 

187. c 

.000184 

.00013? 

.000070 

.18 

187. c 

.000187 

.nnoi 36 

.00007? 

.20 

187.  * 

.000191 

.000140 

.000074 

.?? 

O.O 

.00013* 

.000108 

.0000*8 

.?4 

0.0 

.00010] 

,nn007? 

.00004  7 

•  ?6 

O.n 

.000067 

,nnoo3? 

«000O?Q 

•  ?R 

o.n 

.000030 

,nooo?4 

.000017 

•  ^0 

O.n 

•0000P3 

.nnooi6 

.000010 

.3? 

O.n 

.000015 

.nnoooR 

.000006 

.34 

O.n 

•OOOOOR 

,nnooo4 

.000007 

.36 

o.n 

.000007 

.000004 

.000003 

•  3B 

O.n 

.000007 

.000004 

.000003 

.40 

o  .n 

.000007 

,nooo04 

.000007 

.4? 

o.n 

.000007 

.000004 

.000003 

•  44 

o.n 

.000007 

.000004 

.000003 

.46 

O.n 

.000007 

.000004 

.000007 

•  4fl 

0,0 

.000007 

.nnoooR 

.000007 

.5n 

o.n 

.000007 

.000006 

.000007 

?MT? 


RFCOPD  OF  RAW  0/*TA 


TIME   LOADH.PS.    D  F_  F  I  F  r  T  T  0  M  S.TN. 
(SEC)     X  =  0       X=3,?5      YrS.TS      x=8.?c; 
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0.00 
.0? 
.04 
.06 
.08 

•  10 
.1? 
.14 
.lft 
.1R 

•  ?0 
.?? 
.?4 
.?ft 
.?R 
.^0 
.32 
.34 
.3ft 
.38 
.40 
.4? 

•  44 

•  4ft 
.4R 
.50 


0.0 

0.0 

200.0 

350.0 

375.  n 

375.0 

375.0 

37*.  n 

37S.0 

37*. o 

37*. r 

0.0 

0  ,0 

0.0 

O.o 

0.0 

0.0 

O.o 

o.o 

0.0 
0.0 
O.o 
O.o 

0.0 

o.o 

0  .0 


0.000000 

0.000000 

.00001* 

.00013* 

.ooopo-* 

.000270 
.000330 
.00034* 
.0003ft0 
,O003ft7 
.00037*; 
.000270 
.000172 
.00013* 
.OOOOftO 
.00004^ 
,000030 
.000015 
.00001* 
.00001* 
.00001* 
.00001* 

.ooooi*; 

.00001* 
.00001* 
.00001* 


O.oO 

0.00 

.  00 

.oo 
.on 

.no 
.or 

.00 

.no 
.on 
.00 
.no 
.oo 
.no 
.oo 

.00 
.00 
.00 
.00 
.00 
.00 
.00 
.00 
.00 

.on 

.00 


0000 
0000 
00OP 
0072 
0108 
017ft 
0212 
0?48 
OPftP 
0272 
0?7ft 
017ft 

010R 
007? 
003? 
0024 

onift 

000P 
OOOB 

ooofi 

0008 
0008 
0008 
OOOR 
OOOR 
OOOR 


0.000000 
0.000000 
.OOOOOP 
,00007? 
.00010P 
,0001?& 

.ooon? 

,00013ft 
,000140 
.000144 
.000144 
.000108 
,00007? 
.0000^? 
,0000?4 
,00001ft 
.OOOOOP 
.000004 
.000004 
.000004 
.000004 
.000004 
.OOOOOA 
.000004 
.000004 
.000004 
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?MT? 

PESPONSF    FUNCTIONS   FROM    TmP|  tttt    CONVOLUTION 
TIMF  RFSPONSF      rlJNfTT0N9 

'SEC)  Xs3.?*  X=^.75  X=8.2C 


.000 

.00000020 

.00000011 

.0000001  1 

.020 

.000001 78 

,000nO09c> 

•OOOOOOcc 

.040 

.00000251 

.000001 34 

.000001  "54 

.060 

.00000317 

.00000212 

.00000142 

.080 

.0^000365 

•00000239 

.000001 19 

.100 

.00000349 

.00000263 

.000001 10 

.120 

.00000334 

.000002*4 

.000001 ?2 

•  140 

•00000311 

.00000243 

.000001  i=; 

•  160 

•00000289 

.00000224 

.000001 04 

•  1«0 

.00000122 

.00000069 

.0000004* 

•  200 

-.00000003 

-.00000020 

.00000003 

•  220 

-.00000029 

-. 000ft00^3 

-.00000018 

.240 

-.00000094 

-.OOOOOOBO 

-.00000031 

•  260 

-.00000069 

-.00000059 

-. 00000024 

.280 

-.000000*7 

-.00000038 

-.OOOOOOiq 

.300 

-.00000028 

-.0000001 8 

-.0000001 1 

.320 

.00000006 

.0000000" 

.00000002 

.340 

.00000034 

.00000030 

.00000012 

.360 

.00000043 

.00000034 

.00000016 

.380 

.00000039 

,00000079 

.00000015 

.400 

.00000032 

.00000020 

.00000009 

.420 

.00000019 

.0000001 0 

.00000005 

.440 

.00000010 

.00000003 

.00000002 

.460 

.00000005 

-.00000001 

.00000000 

.480 

-.00000009 

-.00000008 

-.00000004 

.500 

-.00000017 

-.O00n00l2 

-.0000000* 
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?MT? 

DEFLECTIONS  FRO*  RESPONSE  FUMCTTONS  ANp  I  D/>0 

TIME   LOAD»LPS.    OFF!  F  r  T  I  0  W  5»IM. 
(SEC)      X  =  0       X=3.?5      Y=*.7*      X*8.?5 


.00 

0.0 

0.000000 

o .oooooo 

0.000000 

.0? 

0.0 

0.000000 

0.000000 

0.000000 

.04 

?00.0 

.000015 

.P00008 

.000008 

.0* 

350.0 

.00013* 

.A0007? 

.00007? 

.08 

375.0 

.000?0? 

.OO01P8 

.00010" 

.10 

375.0 

,000?70 

.000176 

,0001?4 

.1? 

37*.  0 

.000330 

.000?! ? 

.0001"»? 

.14 

375.0 

.000345 

,000?48 

.0001-36 

.1ft 

375.0 

.000360 

.ppppfip 

.000140 

.IS 

375. p 

.000367 

,000?7? 

.000144 

.?p 

375. r 

.00037* 

,OP0?76 

.000144 

.?? 

0.0 

.000770 

,0001 76 

.000108 

.24 

0.0 

.000173 

.00010* 

.00007? 

.?6 

O.p 

•P0013* 

.00007? 

.0000"'? 

•  ?8 

0.0 

.000060 

.00003? 

.0000?4 

..10 

0.0 

.00004* 

.0000?4 

.000016 

.3? 

0.0 

.000030 

.000016 

.OOOOOP 

.34 

0.0 

.000015 

.OOOOOP 

.000004 

.36 

0.0 

.000015 

.000008 

.000004 

.38 

0.0 

.00001* 

.000008 

.000004 

•  4n 

0.0 

.000015 

.000008 

.000004 

.4? 

0.0 

.000015 

.D00008 

.000004 

.44 

0.0 

.000015 

.OOOOOP 

.000004 

.46 

0.0 

.00001* 

.000008 

.000004 

•  4R 

0.0 

.000015 

.000008 

.000004 

.50 

0.0 

.000015 

.000008 

.000004 

?MI3 


RECORD    OF    PAU     OATA 


TIME       LOAOtlPS.  OFFl     FfTTON    S.T.N. 

(SEC)  X  =  0  X=3.?*  y=*,7*  XsP.?6; 


183 


0.00 

O.P 

0.000000 

0.000000 

0.000000 

.0? 

0,0 

0,000000 

O.onoOOO 

0.000000 

.04 

400. n 

,000090 

.opoopo 

.000040 

.06 

700.0 

.000675 

.00043? 

.000?0O 

•  OP 

750. n 

,000750 

.00047? 

,ooo??« 

.10 

750.0 

.0007P0 

.0004RO 

.000?3? 

.1? 

750.0 

.000R10 

.0004PP 

.OOO?"^ 

•  14 

750.0 

.000P17 

.000496 

.000?40 

.16 

750.0 

,O00R?5 

.000496 

,000?44 

.IP 

750,0 

.OO0P32 

.000*04 

•000?4« 

.20 

750.0 

•OO0R4O 

.00051? 

.000?5? 

.?? 

0.0 

.000600 

.O003?0 

.OOOlf 0 

.24 

P.O 

.000300 

,000?OP 

.OOOOPO 

.26 

0.0 

.000150 

.000160 

,000040 

.28 

0.0 

,000075 

.ftOOOPO 

.00003? 

.30 

0.0 

,000060 

.00004P 

.0000?4 

.3? 

0.0 

,000045 

.00003? 

.000016 

.34 

0.0 

,000030 

,0000?4 

.OOOOOP 

.36 

0.0 

,000015 

.000016 

.000004 

.38 

0.0 

.000015 

.OOOOOP 

.000004 

.40 

0,0 

.000015 

.noonoP 

.000004 

.42 

0.0 

.000015 

.OOOOOP 

.000004 

•  44 

0.0 

.000015 

.OOOOOP 

,000004 

.46 

0.0 

.000015 

.OOOOOP 

,000004 

.4P 

0.0 

.000015 

.OOOOOP 

,000004 

.50 

0.0 

.000015 

•OOOOOP 

,000004 
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2MT3 

RFSPONSF    FUNCTIONS    FROM    ImP|  TCIT    CONVOLUTION 
TIME  PFSPONSF      FUNCTIONS 

(SEC)  X«3.25  X^.^5  X«8.?5 


.000 

•oooooo^o 

.000000^3 

.00000027 

.020 

.00000444 

.00000283 

.000001^1 

.040 

.00000453 

.00000283 

.00000137 

.060 

.00000427 

.00000260 

.000001 26 

.080 

.00000404 

.00000239 

.000001 16 

•  TOO 

.00000369 

.00000220 

.000001 07 

.1?0 

.00000337 

.00000198 

.O0OO00Q9 

.140 

.00000308 

.00000184 

.00000002 

tun 

.00000282 

.00000170 

.00000085 

•  lpo 

.00000099 

.00000011 

.00000018 

•  200 

-.00000067 

-.00000020 

-.oooooo?*; 

•  220 

-.000001  16 

-.00000022 

-.00000035 

•  ?40 

-.00000112 

-.00000047 

-.U00000P5 

.260 

-.00000071 

-.00000040 

-.00000016 

.280 

-.00000036 

-.00000024 

-.00000009 

.100 

-.00000009 

-.00000007 

-.00000003 

.120 

.00000013 

.00000006 

.00000004 

.340 

•0000004O 

.00000018 

.00000012 

.360 

.00000046 

.00000019 

.00000012 

.3A0 

.00000035 

.00000015 

.00000009 

•  400 

.00000020 

.00000011 

.00000004 

.420 

.00000006 

.00000006 

.00000001 

.440 

-.00000001 

.00000001 

-.00000000 

.460 

-.00000005 

-.00000002 

-.00000001 

.480 

-.00000010 

-.00000005 

-.00000003 

.500 

-.00000013 

-.00000006 

-.00000003 
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DEFLECTIONS    FROM    RESPONSE    FUNCTIONS    AND    LOAO 

TIMF       LOAD*LPS.  DEFL    ErTION    S.IN. 

(SEC)  X=o  X=3.?5  X=5.75  X*fl.?^ 


.00 

0.0 

0.000000 

o.oooooo 

0.000000 

.02 

0.0 

0.000000 

0.000000 

0.000000 

.04 

400.0 

.000090 

.000080 

.000040 

.Oft 

700.0 

.00067*1 

.000432 

.000200 

•  OR 

750.0 

.000750 

.000472 

.000228 

.10 

750.0 

.000780 

.000480 

.000732 

.1? 

750.0 

.000810 

.000488 

.000716 

.14 

750.0 

.000817 

.000496 

.000740 

.16 

750.0 

.0008?*; 

.000496 

.000244 

.1* 

750.0 

.000832 

,oooc;o4 

.000248 

.20 

750.0 

.0^30840 

.000S12 

.000752 

.22 

0.0 

.000600 

.000320 

.000160 

.24 

0.0 

.000300 

.000208 

.000000 

.2ft 

0.0 

.O001SO 

.000160 

.000040 

.28 

0.0 

.00007*; 

.000080 

.000012 

.30 

0.0 

.000060 

.000048 

.000024 

.32 

0.0 

.000045 

.000032 

.000016 

.34 

0.0 

.000030 

.000024 

.000008 

.36 

0.0 

.00001S 

.000016 

.000004 

.38 

0.0 

.000015 

.000008 

.000004 

.40 

0.0 

.ooooi*; 

.000008 

.000004 

.42 

o.o 

.000015 

.000008 

,000004 

.44 

0.0 

.00001^ 

.OOD008 

.000004 

.4ft 

0.0 

.000015 

.000008 

.000004 

•  48 

0.0 

.oooois 

.000008 

.000004 

.50 

0.0 

.ooooi6; 

.000008 

.000004 
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?HT1 


PECOPP  OF  PAW  DATA 

TIME      LOAD. LPS.         DEFLECTION    S,IN. 
(SFC)  X=0  X=3.25  x=5.75  X=P.2« 


0.00 

0.0 

0.000000 

O.POOOOO 

0.000000 

.0? 

0.0 

0.000000 

0.000000 

0.000000 

•  04 

100.0 

.OOOOftO 

.000032 

,00001ft 

•  Oft 

175. p 

.000225 

.oooopo 

.000032 

•  OP 

187.5 

.000240 

.000160 

.0000** 

.10 

1R7.5 

.000255 

.000200 

.oooo«o 

.12 

1P7.5 

.000270 

.00020P 

.00009ft 

.14 

1R7.* 

.0002R5 

,00021ft 

.000100 

•  16 

1R7.C 

.000300 

.000224 

.000104 

•  IP 

1P7.5 

.000315 

.000232 

.00010P 

.20 

1R7.5 

.000330 

.000240 

.0001 12 

.22 

O.n 

.000150 

.POOlftO 

.oooopo 

.?4 

0.0 

.000120 

.POOORO 

.000040 

•  2ft 

0.0 

.000090 

.0000ft4 

.000032 

.?P 

0.0 

.000060 

.P0004P 

.000024 

.30 

0.0 

.000045 

.000032 

.00001ft 

.32 

0.0 

.000030 

.000024 

.000012 

.34 

0.0 

.00001  «=. 

.OPOOlft 

.OOOOOP 

.3ft 

0.0 

.000015 

.OOOOOP 

.000004 

.3P 

0.0 

.ooool*; 

.OOOOOR 

.000004 

.40 

0.0 

.000015 

.OOOOOP 

.000004 

.4? 

0.0 

.00001* 

.OOOOOP 

.000004 

•  44 

0.0 

.000015 

.OOOOOP 

.000004 

•  4ft 

0.0 

.000015 

.OOOOOP 

.000004 

.4P 

0.0 

.000015 

.OOOOOP 

.000004 

•  50 

0.0 

. 000015 

.OOOOOP 

.000004 
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?HJ1 


RESPONSE    FUNCTIONS    FROM    IMP|  T<~TT    CONVOLUTION 
TlMf  RFSPOMSF       FUNOTTONS 

(SEC)  X*3.25  y  =  c.7e;  x*8,?s 


1.000  .nroooi^O  .00000085  .00000043 

•  020  .00000585  .00000205  .OOOOOOpI 

.040  .00000569  .00000399  .000001^7 

.060  .00000552  .00000467  .000001P8 

.OflO  .00000537  .0000044?  .00000212 

.100  .00000523  .00000419  .0000020? 

•1?0  .00000511  .000O03QR  .000001Q? 

.140  .00000500  .000003HO  .000001P4 

.160  .00000490  .00000363  ,0000017ft 

.180  -.00000023  .000001??  .00000077 

•200  -.00000046  -.0OPO00P5  -.000000^0 

•??0  -.00000065  -.O0OOO0R1  -.00O00OT5 

.240  -.00000083  -.00000068  -.000000^4 

•260  -.00000061  -.00000060  -.00000031 

.280  -.00000043  -.00000033  -.00000018 

.300  -.00000027  -.00000011  -.00000008 

.320  .00000026  .00000007  .00000001 

.340  .00000072  .00000042  .00000018 

.360  .00000063  .00000051  .00000024 

.380  .0000005?  .000OO0T7  .00000019 

.400  .00000040  .00000025  .00000014 

.420  .00000028  .00000016  .00000009 

.440  .00000019  .00000008  .00000005 

.460  .00000013  .00000004  .00000002 

•4S0  -.00000011  -.OOOOOOOP  -.00000004 

.500  -.00000028  -.00000018  -.OOOOOOOQ 
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?HT1 


DEFLECTIONS  FROM  PFSPONSF  FUNCTION?  ANn  I.OAO 

TIME       LOADtLPS.  DFFl     ErTTOM    StIN. 

(SEC)  X=0  X=3.25  *=5.75  X=fl.?R 


.00 

0.0 

0.000000 

0.000000 

0.000000 

.02 

0.0 

0.000000 

0.000000 

0.000000 

.04 

100.0 

.000060 

.0000^2 

.00001* 

.06 

175.0 

.00022*; 

. oooopo 

.000032 

.08 

187.5 

.000240 

.000160 

.000056 

.10 

1R7.C 

.000255 

.000200 

.OOOOPO 

.1? 

1R7.5 

.000270 

.00020P 

•OOOOQA 

•  14 

lR7.cs 

•0002R5 

.00021* 

.000100 

.16 

1.87.5 

.000300 

.000224 

.000104 

•  1R 

lR7.e? 

.000315 

.000232 

.00010R 

•  20 

lR7.«i 

.000330 

.000240 

.000112 

.22 

o.n 

.000150 

.nOOlfrO 

.OOOOPO 

.24 

0.0 

.000120 

.OOOOPO 

.000040 

.26 

0.0 

.000090 

.n00064 

.000032 

•  2fl 

0.0 

.000060 

.00004P 

.000024 

.10 

0.0 

.000045 

.000032 

.000016 

.32 

0.0 

.000030 

.000024 

.000012 

.34 

0.0 

.000015 

.noOOl* 

.OOOOOR 

.36 

0.0 

.oooois 

.oooooR 

.000004 

•  3fl 

0.0 

.000015 

.OOOOOP 

.000004 

.40 

0.0 

.000015 

.OOOOOP 

.000004 

.42 

0.0 

.000015 

.OOOOOR 

.000004 

.44 

0.0 

.000015 

.OOOOOR 

.000004 

.46 

0.0 

.oooom 

.OOOOOR 

.000004 

•  4R 

0.0 

.000015 

. OOOOOR 

.000004 

.50 

0.0 

.ooooi«=> 

.OOOOOR 

.000004 
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2HT2 

RFCORD  OF  PAW  DATA 

TIMF   LOADtLPS.    DFFI  E  r  T  T  0  N  S.IN. 
(SEC)      X  =  0       X=3.25      y=5.75      X=P.?5 


0.00 

0 

.0 

0.000000 

0.000000 

0.000000 

.02 

0 

»0 

0,000000 

0.000000 

0. 000000 

.04 

200 

•  0 

.000120 

.0000*4 

.000040 

•  Of 

350 

,n 

•0004P0 

.000144 

.000136 

•  OP 

375 

.0 

.000600 

.000272 

.00016P 

.10 

375 

1  p 

.ooo6m 

.O003ci2 

.0001*3? 

.1? 

375 

.0 

.000630 

.000432 

.00020P 

.14 

375 

ifl 

.000645 

.00044P 

.00021? 

.1* 

375, 

>0 

.000660 

.000464 

.00021* 

.IP 

375, 

,n 

.000675 

.000472 

,0002?0 

.20 

375, 

.0 

.0006QO 

•0004PO 

.000224 

.22 

0 

.0 

.000330 

.000352 

.00017* 

.24 

0 

.0 

.000210 

.000272 

.000116 

.26 

0 

.0 

.0001PO 

.000192 

.00009* 

.2P 

0, 

.0 

.000150 

.000112 

.0000*4 

.30 

0, 

,n 

.000120 

.OOOOPO 

.00004P 

.32 

0, 

.0 

.000090 

,000064 

.00003? 

.34 

o, 

,0 

.000075 

.00004P 

,0000?P 

.36 

0, 

n 

.000060 

.000040 

,0000?4 

•  3P. 

0, 

.0 

.000045 

.000032 

.000020 

•  40 

0, 

.n 

.000030 

,000024 

.000016 

.42 

0, 

,0 

.ooooi*; 

.00001* 

.000012 

•  44 

0, 

.0 

.000015 

.OOOOOP 

•OOOOOP 

.46 

0, 

.n 

.000015 

.OOOOOP 

.000004 

•  4P 

0, 

»0 

.00001^ 

, OOOOOP 

.000004 

.50 

0, 

,0 

.000015 

.OOOOOP 

.000004 
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?HI? 

RESPONSE    FUNCTIONS    FROM    T^R|  T<~TT    CONVOLUTION 
TIME  RFSPONSF       FUNCTIONS 

(SEC)  X=3.25  X  =  <=.75  Xsfl.?" 


.000 

.00000160 

.oOonooR^ 

.000000*3 

.020 

.00000625 

,00000lfl4 

,00000176 

.040 

•0O000726 

.00000337 

.00000202 

•  060 

.00000674 

.0000041 1 

.00000214 

.080 

.00000626 

.00000477 

.00000214 

•  100 

•00000583 

.00000451 

.00000108 

•  120 

•00000545 

•00000427 

•00000183 

•  140 

•00000510 

•00onO3O5 

.00000170 

•  160 

•00000479 

•00000366 

•OOOOOICQ 

•  1*0 

-•000000-*3 

.00000167 

•000000P4 

•  200 

-.00000131 

.00000061 

.00000039 

•  220 

-.00000086 

-.00000019 

.00000002 

•  240 

-.00000050 

-.OOOOOOR^ 

-.00000020 

•  260 

-.0000002? 

-.00000070 

-.00000018 

•  280 

-.00000001 

-.OOOOOOOO 

-.00000017 

.300 

.00000033 

-.00000014 

-.00000003 

•  3?0 

.00000061 

.00000016 

•OOOOOOOO 

.340 

.00000083 

.00000041 

.00000019 

•  360 

.0000005? 

.00000043 

.00000020 

.380 

.00000013 

.00000034 

.0000001 7 

.400 

.00000003 

.00000018 

.00000010 

.420 

-.0000000? 

.00000008 

.00000002 

.440 

-.00000004 

-.oooooooo 

-.oooooooo 

.460 

-.00000004 

-.00000004 

-.00000002 

.480 

-.0000001 0 

-.0000001 0 

-.00000005 

.500 

-.00000013 

-.00000013 

-.00000006 
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?HT? 

DFFLECTIONS  FRCM  PFSPONSF  FUNCTIONS  ANP  LDAn 

TIMF      LOAD, |  PS.         DEFLEcTTOM    S.IN, 
(SEC)  X  =  0  X  =  3.?5  V=5.75  X=P.?* 


.00 

0.0 

0.000000 

0.000000 

0.000000 

.02 

0.0 

0.000000 

O.ooOOOO 

0.000000 

.04 

200.  0 

.000120 

. 000064 

.000040 

.06 

350.  n 

.0004P0 

, 000144 

.000136 

•  OR 

375.0 

.000600 

.00027? 

.00016P 

.10 

375.0 

.000615 

.000352 

.000192 

.1? 

375.0 

.000630 

.00043? 

.00020« 

.14 

375.0 

.000645 

.OO044P 

.00021? 

.16 

375.0 

.000660 

.000464 

.000216 

.IP 

375.0 

.000675 

.00047? 

.000220 

.20 

375.0 

.000690 

.0004PO 

.000224 

.22 

0.0 

.000330 

.00015? 

.000176 

.24 

0.0 

.000210 

,000?7? 

.00013* 

.26 

O.n 

.0001PO 

.000192 

.00009* 

.28 

0.0 

.000150 

.000112 

.0000*4 

.30 

0.0 

.000120 

.oOOOflO 

.00004R 

.32 

O.o 

.000090 

.0000*4 

.00003? 

.34 

0.0 

.000075 

,ro004« 

•0000?R 

.36 

0.0 

.000060 

.000040 

.0000?4 

.3fi 

0.0 

.000045 

,000032 

.000020 

.40 

0.0 

.000030 

.000024 

.000016 

.42 

0.0 

.000015 

.000016 

.000012 

.44 

0,0 

.000015 

,«00O0P 

.OOOOOP 

•  46 

0.0 

.000015 

.OOOOOP 

.000004 

.4P 

0.  0 

.000015 

.OOOOOP 

.000004 

.50 

0.0 

.000015 

.OOOOOP 

.000004 
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?HP 


RECORD  OF  PAW  DATA 


TIME      LOAD. LPS.         DEFLEpTTON    S.IN, 
(SEC)  X  =  p  X=3.25  Xr5.75  X  =  fi.?e; 


.00  0,0  0.000000  O.nOOOOO  0.000000 

0?  0.0  0.000000  0.000000  0.000000 

04  400.0  .000375  .fl00224  .000064 
06  700.0  .000825  .000736  .000240 

05  750.0  .000863  .000848  .0003?? 
10  750.0  .000975  .00091?  .000400 
1?  750.0  .001050  .OP09?R  .000416 
14  750.0  .OOlOPfl  . 000944  .000424 
16  750.0  .0011?5  ,000960  .0004?? 
IP  750.0  .001163  .000976  .000440 
20  750.0  .0/01200  .0009^2  .00044** 
2?  0.0  .000525  .000640  .000320 
24  0.0  .000450  .00044P  .000240 
26  O.o  .000375  .000256  .000160 

>28  0.0  .000300  .000160  .000096 

30  0.0  .000225  .000128  .OOOOPO 

32  0.0  ,^)0018P  .000096  .000064 

34  0.0  .000150  .000064  .00004R 

i36  0.0  .000113  .000048  . 000012 

>3B  0.0  .000075  .000032  .000024 

,40  0,0  ,000075  .000016  .000016 

42  0.0  ,00007ci  .000016  .OOOOOP 

44  0.0  ,000075  .000016  .OOOOOP 

,46  0.0  .000075  .000016  .OOOOOP 

48  0.0  .000075  .000016  .OOOOOP 

,50  0.0  .000075  .000016  .OOOOOP 
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RESPONSE  FUNCTIONS  FROM  TmP|  TCIT  CONVOLUTION 
TIME        PESPONSF   FUNCTIONS 
(SEC)        X*3.25         X=<=.75         X  =  fi.?S 


0*000 
.020 

•  040 

•  060 
.080 
.100 
.120 

•  140 
.160 
.IPO 
.200 
,220 
.240 
.260 
.280 
.300 
.320 
.340 
.360 
.380 

•  400 
.420 
.440 
.460 
.4*0 
.500 


.00000250 

.00000527 

.00000501 

.00000526 

.00000523 

.00000496 

.00000471 

.00000449 

.00000429 

-.00000039 

-.00000036 

■.00000032 

-.00000026 

-.00000021 

.00000005 

.00000027 

.00000045 

,00000058 

.00000048 

.00000040 

.00000033 

.00000027 

.00000022 

.00000020 

■• 00000004 

•.00000024 


.00000149 

,0O0n0477 

.0000050ft 

,O0OnO4Q8 

.00000459 

.00000423 

.0000039? 

.00000363 

.00000337 

.00000083 

-.00000007 

..OOonOOP4 

-.0000009] 

-.00000057 

-.00000030 
-.00000009 

.00000018 

.000  00039 

.00000033 
.00000029 

•0000001P 
•000ft0007 

.00000000 
.00000003 
.00000009 
.OOooOOl 1 


.00000043 
.000001** 

.00000216 
.00000227 
.00000215 
.00000109 
.000001P4 
.00000171 
.000001*9 

. 000000*2 

.00000017 

-.00000016 

-.00000035 

-.00000021 

-.00000010 

-.00000001 

.00000006 

,0000001ft 

,00000015 

,00000010 

,00000007 

,00000003 

,00000001 

-.00000000 

-.00000003 

-.00000005 
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?HI3 

DEFLECTIONS  FROM  RESPONSF  FUNCTION?  AND  LOAD 

TIME      LOADtLPS.         DFFLErTIOM    StIN. 
(SEC)  X=0  X*3.25  Xe5.75  X*P.2«< 


•  00 

0.0 

0.000000 

0.000000 

0.000000 

.^2 

0.0 

0.000000 

0,000000 

0,000000 

.04 

400.0 

.000375 

.000224 

.000064 

•  06 

700.0 

•  00082«S 

.000736 

.000240 

.08 

750.0 

.000*63 

.00084* 

.00035? 

.10 

750.0 

.000975 

.00091? 

.000400 

.1? 

750.0 

.001050 

.000928 

.000416 

•  14 

750.0 

.001088 

,000944 

.000424 

.16 

750.0 

.001125 

.000960 

.00043? 

.18 

750.0 

.001163 

,000976 

.000440 

.20 

750.0 

.001200 

.000992 

,000448 

•  22 

0.0 

.000525 

•  000640 

,0003?0 

.24 

0.0 

.000450 

.000448 

,000?40 

.26 

0.0 

.000375 

.000256 

.000160 

.28 

0.0 

.000300 

.000160 

.000096 

.30 

0.0 

.000225 

.000128 

.0000*0 

.32 

O.p 

.00^187 

.000096 

,0000f4 

.34 

0.0 

.000150 

.000064 

.000048 

.36 

0.0 

.000112 

.000048 

.00003? 

.38 

0.0 

.000075 

.00003? 

.0000?4 

.40. 

0.0 

.000075 

.000016 

.000016 

•  42* 
.44 

0.0 
0*^ 

.000075 
.000075 

.000016 

AAAA1  ^ 

.000008 

flAAAAQ 

•  0000 1 o 

. UuUOOn 

.46 

0.0 

.00007*; 

.000016 

.000008 

•  48 

0.0 

.000075 

.000016 

.000008 

•  50 

0.0 

.000075 

.000016 

.000008 

APPENDIX  B 
STATIC  LOAD  PROGRAM 
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PROGRAM  STALOD( INPUT, OUTPUT, TAPE5= INPUT, TAPE6=OUTPUT) 

C 

r      ****************************************************** 

C 

C      COMPUTER  PROGRAM  TO  PREDICT  STATIC  LOAD  RESPONSE 

C 

r  ****************************************************** 

C 

C      THE  PROGRAM  COMPUTES  THE  STEADY  STATE  RESPONSE.  THE 

C      INPUTS  ARE  THE  LOAD  MAGNITUDE  AND  THE  PARAMETERS 

C      DETERMINED  FROM  IMPULSE  TESTS. 

C      K=  CONVERSION  FACTOR.  K=  10  IN  THIS  STUDY. 

C      TN=  TIME  AT  WHICH  IMPULSE  INPUT  FIRST  BECOMES  NONZERO 

C      TN=  .04  IN  THIS  STUDY. 

C      NS=  NUMBER  OF  CASES  ANALYZED 

C      FO=  FORCE  MAGNITUDE  IN  POUNDS 

C      STPRSP=STEP  LOAD  RESPONSE  TO  BE  CALCULATED 

C      ALPHA,  BETA,  AND  GAMMA  ARE  PAVEMENT  PARAMETERS 

C      THE  OTHER  VARIABLES  ARE  SELFDESCRI PTIVE  AND  REFER  TO 

C      THE  STATIC  LOAD  EQUATION  IN  THE  TEXT 

C 

C 

DIMENSION  LABLE(4) 

INTEGER  LABLE 

REAL  K 

READ(5,1)  NS,K,TN 

DO  99  1=1, NS 

READ(5,2)  (LABLE(J)  ,J  =  1,3)  ,FO 

READ(5,3)  ALPHA, BETA, GAMMA 

DEN=BETA**2+GAMMA**2 

ROTDEN=SQRT(DEN) 

TANPH I = GAMMA/ BETA 

PHI=ATAN(TANPHI) 

ARGSIN=GAMMA*TN+PHI 

STPRSP=FO*ALPHA* (K*GAMMA/DEN- (K-l . ) /ROTDEN*EXP (- BETA* 
1TN)*SIN(ARGSIN)) 
C      MULTIPLY  RESULT  BY  10000  FOR  PRESENTATION 

STPRSP=STPRSP*10  000. 
C      THE  FOLLOWING  IF  STATEMENTS  ALLOW   PRINTING  CERTAIN 
C      GROUPS  OF  RESULTS  CN  ONE  PAGE 

IF(I.EQ.  1)  GO  TO  20 

IF(I.EQ.  7)  GO  TO  20 

IF(I.EQ.16)  GO  TO  20 

IF(I.EQ.25)  GO  TO  2  0 

IF(I.EQ.34)  GO  TO  20 

IF(I.EQ.43)  GO  TO  20 


196 


WRITE(6,5)  (LABLE(M) ,M=1,3) ,F0 , ALPHA , BETA, GAMMA , STPRSP 

GO  TO  99 

20  WRITE(6,4) 

WRITE (6,5)  (LABLE(M) ,M=1,3) , FO , ALPHA, BETA , GAMMA , STPRSP 

99  CONTINUE 
p      ****************************************************** 

C 

C      FORMAT  STATEMENTS 

C 

Q         ****************************************************** 

c 

1  FORMAT(I3,2F10.5) 

2  FORMAT (3A6,2X,F10. 2) 

3  FORMAT(E18.8,2F18.8) 

4  FORMAT (1H1,///9X, 'SERIES' ,8X, 'LOAD*  ,11X, 'P  A  R  A  M  E  * 
1,'T  E  R  S',7X,'STEP  RESP '// 23X ,' LBS .', 7X ,' ALPHA ', 7X , 

2 'BETA' ,6X, 'GAMMA' , 2X,'IN.*1/10000') 

5  FORMAT(//3A6,F10.2,E13.3,3F10.3) 
C 

STOP 
END 


APPENDIX  C 
REPEATED  LOAD  PROGRAM 


c 
c 
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PROGRAM  REPLOD  (INPUT , OUTPUT ,TAPE5= INPUT ,TAPE6=OUTPUT) 
C 

Q  ****************************************************** 

c 

C  COMPUTER  PROGRAM  TO  PREDICT  REPEATED  LOAD  RESPONSE 

C 

r  ****************************************************** 

C 

C  THE  PROGRAM  COMPUTES  THE  TIME-DEPENDENT  RESPONSE  OF 

C  THE  SYSTEM.  THE  INPUTS  ARE  THE  RATE  (OR  DURATION)  AND 

C  THE  AMPLITUDE  OF  THE  LOADING  FUNCTION,  AND  THE  PARAMETERS 

C  DETERMINED  FROM  IMPULSE  TESTS. 

C 

C  K=  CONVERSION  FACTOR.  K=  10  IN  THIS  STUDY. 

C  TN=  TIME  AT  WHICH  IMPULSE  INPUT  FIRST  BECOMES  NONZERO. 

C  TN=  .04  IN  THIS  STUDY. 

C  FO=  AMPLITUDE  OF  REPEATED  LOADING  FUNCTION. 

C  P=  PERIOD  OF  LOADING. 

C  ALPHA,  BETA  AND  GAMMA  ARE  PAVEMENT  PARAMETERS. 

C  NS=  NUMBER  OF  CASES  ANALYZED. 

C  NP=  NUMBER  OF  PREDICTIONS  DESIRED  FOR  EACH  CASE. 

C  N=  NUMBER  OF  LOAD  APPLICATIONS. 

C  TPEAK=  TIME  WHEN  THE  LOAD  IS  AT  PEAK  POSITION  AFTER  N 

C  LOAD  APPLICATIONS. 

C  RPEAK=  RESPONSE  WHEN  THE  LOAD  IS  AT  PEAK  POSITION 

C  AFTER  N  LOAD  APPLICATIONS. 

C  TOFPK=  TIME  WHEN  LOAD  IS  AT  OFF  PEAK  POSITION  AFTER  N 

C  LOAD  APPLICATIONS. 

C  ROFPK=  RESPONSE  WHEN  LOAD  IS  AT  OFF  PEAK  POSITION 

C  AFTER  N  LOAD  APPLICATIONS. 

C  THE  OTHER  VARIABLES  ARE  SELFDESCRI PTI VE  AND  REFER  TO 

C  THE  REPEATED  LOAD  EQUATIONS  IN  THE  TEXT. 


DIMENSION  LABLE(ll) 

INTEGER  LABLE 

REAL  K 

READ(5,1)  NS,NP,P,K,TN 

DO  200  M=l ,NS 

WRITE  (6,2) 

READ  (5,3)  (LABLE(J) ,J=1,11) ,FO 

READ (5, 4)  ALPHA,  BETA,  GAMMA 

WRITE(6,5)  (LABLE(J)  ,J=1,11) 

WRITE(6,6)  FO 

WRITE (6, 7)  ALPHA, BETA, GAMMA 
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WRITE(6,8)  P 

WRITE (6, 9) 

DEN=BETA**2+GAMMA**2 

RTDEN=SQRT(DEN) 

TNPH I = GAMMA/ BETA 

PHI=ATAN(TNPHI) 

TWOPI=2.*3. 14159265 

OMEGA=TWOPI/P 

DEN1=BETA**2+ (GAMMA-OMEGA) **2 

RTDEN1=SQRT(DEN1) 

TNPHI1-C GAMMA- OMEGA) /BETA 

PHI1=ATAN(TNPHI1) 

DEN2=BETA**2+(GAMMA+OMEGA)**2 

RTDEN2=SQRT(DEN2) 

TNPHI2=(GAMMA+OMEGA)/BETA 

PHI2=ATAN(TNPHI2) 

N=l 

DO  50  1  =  1, NP 

TPEAK=(2*N-l)*P/2. 

TOFPK=N*P 

RPEAK=  FO*ALPHA/ 2 . * (K*GAMMA/DEN+K/ 2 . * ( (GAMMA- OMEGA) /DEN1 
l+CGAMMA+OMEGA)/DEN2)-(K-l.)/RTDEN*EXP(-BETA*TN)*SIN( 
2GAMMA*TN+PHI)-(K-1.)/2.*(EXP(-BETA*TN)/TRDEN1*SIN(( 

3GAMMA-0MEGA)*TN+PHI1)+EXP(-BETA*TN)/RTDEN2*SIN((GAMMA+ 
40MEGA)*TN+PHI2))-EXP(-BETA*TPEAK)/RTDEN*SIN(GAMMA*TPEAK 
5  +  PHI ) -  5*  (EXP (- BETA*TPEAK) /RTDEN1*SIN ( (GAMMA-OMEGA) * 
6TPEAK+PHI1)+EXP(-BETA*TPEAK)/RTDEN2*SIN((GAMMA+0MEGAJ* 

7TPEAK+PHI 2 ) ) ) 

ROFPK=  FO*ALPHA/ 2 . * ( K* GAMMA/ DEN  - K/ 2 . * ( ( GAMMA -  OMEGA) /DEN 1 
l+CGAMMA+OMEGA)/DEN2)-CK-l.)/RTDEN*EXP(-BETA*TN)*SIN( 

2GAMMA*TN+PHI ) + (K-l . )/2 . *EXP (- BETA*TN) /RTDEN1*S IN ( ( 
3GAMMA-0MEGA)*TN+PHI1)+EXP(-BETA*TN)/RTDEN2*SIN((GAMMA+ 

40MEGA) *TN+PHI 2) ) -EXP ( -BETA*TOFPK) /RTDEN*SIN (GAMMA*TOFPK 
5  +  PHD+  5*  (EXP  (-BETA*TOFPK)/RTDENl*S  IN  ((GAMMA- OMEGA)* 
6TOFPK+PHIl)+EXP(-BETA*TOFPK)/RTDEN2*SIN((GAMMA+OMEGA)* 

7TOFPK+PHI2))) 
WRITE (6, 10)  TPEAK,RPEAK,TOFPK,ROFPK 

50  N=N+1 
200  CONTINUE 

Q  ****************************************************** 

c 

C      FORMAT  STATEMENTS 

q  ****************************************************** 

C 

1  FORMAT(2I3,3F10.5) 

2  FORMAT (1H1) 

3  FORMAT(HA6,F10.1) 

4  FORMAT (El  8. 8, 2F1 8. 8) 

5  FORMAT (11A6) 
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6  FORMAT C///6X, 'THE  LOAD  MAGNITUDE  =*,F10.1,'   LBS  .  '  J 

7  FORMAT (7//6X, 'THE  PARAMETERS  ARE »  ,//8X , ■ ALPHA= '  ,E12 . 3  , 
1'    BETA='  ,F8.3,  '    GAMMA= '  , F8 . 3) 

8  FORMAT (//6X, 'THE  PERIOD  =',F10.5,'   SEC.*) 

9  FORMAT(///8X, 'TPEAK' ,4X, 'PEAK  RESPONSE ', 5X , ' TOFPK *, 5X , 
l'OFFPK  RESPONSE') 

10  FORMAT (//5X,E9. 2, El 4. 3, El  4. 2, El  5. 3) 
C 
q  ****************************************************** 

C 

STOP 
END 
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